Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_013198	1,873,078	Δ1 bp	coding (727/795 nt)	LGG_01843	hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_013198	1,873,078	0	T→.	100.0%	296.8	78	coding (727/795 nt)	LGG_01843	hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/56):  ref base (0/10):  total (12/66)

TGTGCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTAGGCGGTGATGGCC.GGGCTTTGGCCATTGCGGCTGAGGTCCTTACACGCAGACTTCTGCGCCGGCGAACGCGTT.T.ATGAGCGTGAACCGACGCGCTTAGAAACCGG  >  NC_013198/1872982‑1873163                                                                                                     |                                                                                        tgtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCa                                                                                    <  23:297551/101‑1 tgtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCa                                                                                    <  1:481407/101‑1 tgtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCa                                                                                    <  8:297551/101‑1 tgtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCa                                                                                    <  16:481407/101‑1  ttgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  6:508298/2‑101  ttgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  21:508298/2‑101  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  28:909226/1‑101  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   <  8:772647/101‑1  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   <  6:836210/101‑1  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  8:494095/1‑101  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  13:909226/1‑101  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   <  21:836210/101‑1  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   <  23:772647/101‑1  gtgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAt                                                                                   >  23:494095/1‑101   tgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAtt                                                                                  <  19:157602/101‑1   tgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAtt                                                                                  <  4:157602/101‑1   tgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAtt                                                                                  <  20:326251/101‑1   tgCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCAtt                                                                                  <  5:326251/101‑1    gCAGGTTTCTGGGCTGGGGAACGCGTTATG...AGCGTGAACCTGCGCGCTTAGGGATCGGAGTGGANNAGGCCNNNGGCG.NGATGNCCTGGGCTT.GGNCNTTg                                                                                 >  27:119529/1‑101    gCAGGTTTCTGGGCTGGGGAACGCGTTATG...AGCGTGAACCTGCGCGCTTAGGGATCGGAGTGGANNAGGCCNNNGGCG.NGATGNCCTGGGCTT.GGNCNTTg                                                                                 >  12:119529/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTg                                                                                 >  22:969714/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTg                                                                                 >  7:969714/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTg                                                                                 >  18:455080/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTg                                                                                 >  3:455080/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGGGTAAGTGGGCCTGGGC.GTGATGGCCCGGGCTT.GGCCATTg                                                                                 >  5:1202439/1‑101    gCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGGGTAAGTGGGCCTGGGC.GTGATGGCCCGGGCTT.GGCCATTg                                                                                 >  20:1202439/1‑101     cAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGc                                                                                <  8:886406/101‑1     cAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGc                                                                                <  23:886406/101‑1     cAGGTTTCTGCGCTGGTGAACGCGTTATGGTGAGCGTTAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAt                                                                                   <  19:731293/101‑1     cAGGTTTCTGCGCTGGTGAACGCGTTATGGTGAGCGTTAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAt                                                                                   <  4:731293/101‑1      aGGTTTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAtt                                                                                  <  28:849064/101‑1      aGGTTTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAtt                                                                                  <  13:849064/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  12:815406/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  27:815406/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  21:45761/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  16:440519/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  1:979191/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  16:979191/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  1:440519/101‑1      aGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCg                                                                               <  6:45761/101‑1      aGGTCTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAtt                                                                                  <  4:1132925/101‑1      aGGTCTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCGGAGCATAAGTGGCCTTGGGCG.TGATGACC.GGGCTTTGGTCAtt                                                                                  <  19:1132925/101‑1       ggTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCga                                                                              >  19:165098/1‑100       ggTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCga                                                                              >  4:165098/1‑100                gCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGGCTTGGGCGG.GATGGCCTGGGCTT.GGCCATTGCGACCAAGGTcc                                                                     >  1:89383/1‑101                gCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGGCTTGGGCGG.GATGGCCTGGGCTT.GGCCATTGCGACCAAGGTcc                                                                     >  16:89383/1‑101                gCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTcc                                                                     <  6:845145/101‑1                gCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTcc                                                                     <  21:845145/101‑1                  tGGCGGACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  29:200532/101‑1                  tGGCGGACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  14:200532/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  28:787787/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  1:904841/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  4:1027104/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  4:76013/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  1:146457/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  19:76013/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  19:1027104/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  16:904841/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  16:146457/101‑1                  tGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCtt                                                                   <  13:787787/101‑1                   ggCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTa                                                                  <  24:266188/101‑1                   ggCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTa                                                                  <  9:266188/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  10:929062/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  28:838392/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  12:652615/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  27:652615/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  25:929062/101‑1                    gCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTac                                                                 <  13:838392/101‑1                     cGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTaca                                                                <  22:654494/101‑1                     cGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTaca                                                                <  7:654494/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  3:1204742/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  4:920235/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  23:40384/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  19:920235/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  8:40384/101‑1                      gAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacac                                                               <  18:1204742/101‑1                       aaCGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacact                                                              <  11:771627/101‑2                       aaCGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacact                                                              <  8:801214/101‑2                       aaCGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacact                                                              <  23:801214/101‑2                       aaCGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTGGGCG.TGATGGCCTGGGCTT.GGCCATTGCGACCAAGGTCCTTacact                                                              <  26:771627/101‑2                             ttATG...AGCGTGAACCGACGCGCTTAGAAACCGGAGCATAAGTGGCCTTGAGC.GTGATGGCC.CGCTCTTGGCCATTGCGCTTAAGCTCCTTACACGCAGACt                                                        <  19:1229565/101‑1                             ttATG...AGCGTGAACCGACGCGCTTAGAAACCGGAGCATAAGTGGCCTTGAGC.GTGATGGCC.CGCTCTTGGCCATTGCGCTTAAGCTCCTTACACGCAGACt                                                        <  4:1229565/101‑1                                                                                       gATGACC.GGGCTTTGGTCATTGCGCCCAAGGTCCTTATGTGCAGGTTTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCgg  <  13:724536/101‑1                                                                                       gATGACC.GGGCTTTGGTCATTGCGCCCAAGGTCCTTATGTGCAGGTTTCTGGGCTGGTGAACGCGTTATGGTGAGCGTGAACCAGCCCGGTTAGAAACCgg  <  28:724536/101‑1                                                                                                     |                                                                                        TGTGCAGGTTTCTGGGCTGGCGAACGCGTTATG...AGCGTGAACCGGCGCGCTTAGGGATCGGAGTGTAAGTGGCCTTAGGCGGTGATGGCC.GGGCTTTGGCCATTGCGGCTGAGGTCCTTACACGCAGACTTCTGCGCCGGCGAACGCGTT.T.ATGAGCGTGAACCGACGCGCTTAGAAACCGG  >  NC_013198/1872982‑1873163 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG