Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,930,047	0	G	T	28.1%	91.8 / 11.3	57	G158C (GGC→TGC)	sdaB	L‑serine deaminase II
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base T (16/0);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 4.63e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCC  >  NC_000913/2929805‑2930285                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ccGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTAGTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑Ccagcag                                                                                                                                                                                                                                                <  2:215407/251‑1 (MQ=255)      gATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGCGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGAGGC‑TTTATC‑GTTGATGGAGAGCATTTT‑TG‑CCAGCAGG‑A‑t                                                                                                                                                                                                                                           >  2:74013/1‑249 (MQ=255)                 cATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCAGTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTg                                                                                                                                                                                                                             <  2:204718/251‑1 (MQ=255)                  aTCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGAGGC‑TTTTTC‑GTTTATGTAGGGACTTTT‑GT‑GCAACAGG‑A‑TAGCGCCCCg‑ggtgg                                                                                                                                                                                                                            >  2:356857/1‑247 (MQ=255)                              aTTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATcc                                                                                                                                                                                                              <  1:321828/251‑1 (MQ=255)                                       gtgAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑Gtt                                                                                                                                                                                                     <  1:505989/250‑1 (MQ=255)                                          aaTACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGACGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGG‑TTTTTC‑GTTTGTGAAGAGGATTTT‑GG‑CCAAGAGG‑AGTA‑AGGACCG‑GTTGGAG‑TTTCTT‑ATCCGTACA‑gattgac                                                                                                                                                                                                 >  2:549613/1‑245 (MQ=255)                                                cATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GGTGATTGAGAGAATTTTTGG‑CC‑GCAGGGA‑T‑GCGCACCG‑GGTGAAG‑TTGCTT‑AACCggaacggtacggcaccg                                                                                                                                                                                             >  2:45491/1‑234 (MQ=255)                                                   ggTCGCCTGATGCTGGCACACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTTTC‑GTTGATGAAGAGCATTTT‑TG‑CCAGCCGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATTCGTACA‑GATCAGCAGCCGGTCt                                                                                                                                                                                        >  1:93540/1‑251 (MQ=255)                                                           gATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGGTGAAAAGAATTTT‑GG‑CCCGCAGG‑A‑TAGCGCACCG‑GTTGGAG‑GTCCTT‑ATCCGTACA‑GGTCAGCAGCCGATCTGCaaaaa                                                                                                                                                                                 >  1:243407/1‑250 (MQ=255)                                                             tGCTGGCAACCGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGACAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGG‑TTTATC‑GTTGTTGAAGAGCATTTT‑GGGCCGGC‑GG‑A‑TAGCGCACCG‑GGTGGAG‑TTCCTT‑TTTCGTACA‑Gtttaccagaccggtcgcaacaaaat                                                                                                                                                                              >  2:80819/1‑228 (MQ=255)                                                                tGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑GTTTTT‑GTTGTTGAAGAGCATTTT‑TG‑CCAGCAGG‑A‑GAGCGCACCC‑GTTGAAG‑TTTCTT‑ATCCGTACC‑GTTCAGCAGCCGATCTTCAAAAAAATTGt                                                                                                                                                                           >  1:76260/1‑251 (MQ=255)                                                                 ggCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAAAGCATTTT‑TG‑CCCACCGG‑A‑TAGCGCACCG‑GTTGAAG‑TTTCTT‑TTCCGGAAA‑GTTCAGCAGCCGATCTGCAAAAACAATGTc                                                                                                                                                                          >  1:236767/1‑251 (MQ=255)                                                                  gCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCa                                                                                                                                                                         <  1:172650/251‑1 (MQ=255)                                                                    aaaCGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATCCGTT‑ATGAAGAGACTTTT‑GG‑CCCCCAGG‑A‑GAGCGCACCCCGTTGA‑G‑TTACCT‑ATCCGT‑CACGTTCCGCAGCCCATCTGCCCAAACAttttacg                                                                                                                                                                        >  2:461665/1‑245 (MQ=255)                                                                    aaaCGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGGGGC‑TTTATC‑GATGATGAAGAGCAATTT‑TG‑CCAGCAGG‑A‑TAGCGCACCG‑GGTGAAG‑TTCCTT‑ATCCGTACA‑GGTCAGCCGCCGGTCTGCAAAAAAAATGTCA‑gg                                                                                                                                                                      >  1:344240/1‑251 (MQ=255)                                                                     accGGTCAGCATGAAGTTGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGa                                                                                                                                                                     <  1:450755/249‑1 (MQ=255)                                                                     aaCGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGTTGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCC‑GGTGTAG‑TTCCTT‑ATCCGTACA‑GTTTAGCAGCCGATCTGCCAAAAAATTTTCA‑ggg                                                                                                                                                                     >  1:262434/1‑250 (MQ=255)                                                                       cGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑TGTGATGAAGAGCATTTT‑TG‑GCAGCAGG‑G‑TAACGCACCG‑GGTGAAG‑GTCTTT‑ATACGGACA‑TGTCTGCAGCCGATCTTCAAAAACATTGTCA‑GGaa                                                                                                                                                                    >  2:448786/1‑250 (MQ=255)                                                                         gTCAGCATGAAGTGGAGTTCCCGGTTGATCCGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATA‑GTTGATGAAGAGAATTTT‑TG‑GCAGCAGG‑G‑TAGGGGACCG‑GTTGGAG‑TTCCTTTACCCG‑ACA‑GGTCAGGAGCCGATCTGCAAAAACATTGTTA‑GGAAAcc                                                                                                                                                                 >  1:505150/1‑251 (MQ=255)                                                                           cAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCgg                                                                                                                                                               <  2:93540/251‑1 (MQ=255)                                                                                        aGTTCCCGGTTGAGCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTTCTACTCTATTGGCGGTGGC‑TTTATC‑TTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑AACCGTACA‑TTTCAGCAGCCGATCTGCAAAAACATTTTCA‑GGAAACCGGGCTGTCActctct                                                                                                                                                  <  2:536169/251‑1 (MQ=255)                                                                                                       aGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATT‑GTTTATGAAGAGCATTTT‑TT‑CCAGCAGG‑G‑TAGCGCACCG‑GGTGAAA‑TTCCTT‑TTCCGTACA‑GTTCAGCCGCCGATCTGCAAAAACATTGTCA‑GGAAAACACGCCGTCACTCCCTGGCCtgatggttga                                                                                                                                    >  1:326701/1‑245 (MQ=255)                                                                                                          gCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATACAGTCGTTTACAGCCAGACTTACTACTATATTGGCGGAGGCCTTTTTC‑‑TTGGTAACAAGGATTTT‑GG‑CCAGCCGG‑A‑TAGCGCACCG‑GGTGAAGGTTCCT‑‑ATCCCTACA‑GTTCCACAGCCCATCTGCAACAACATTTTTA‑GGAAACAGGGCTGTCACCCTCTGGCCTGATGATGaaaaa                                                                                                                                 >  1:275647/1‑250 (MQ=255)                                                                                                                 cTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGAGGC‑TTTATC‑GTTGGTGGAGAGCATTTT‑GG‑CCCGCCGG‑A‑TAGCGCACCG‑GTTTGAG‑ATCCTT‑ATCCGTAAA‑GATCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAAACGGGCTGTCACTCTCTGGCCCGATGATGGAAAACGAGCTgg                                                                                                                         >  1:405305/1‑251 (MQ=255)                                                                                                                   tCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGAGGC‑GTTATC‑GTTGATGAAGAGCATTTT‑GG‑GCCGCAGG‑G‑TAGAGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGGTCTGCAAAAACATTGTCA‑GGGAACCGGGCTGTCACTCTCTGGCCTGATAATGAAAAACGAGCTGgcg                                                                                                                       >  1:198237/1‑251 (MQ=255)                                                                                                                       cGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGTCCTGATGATGAAAAACGAGCTGGCGCTGc                                                                                                                   <  1:389022/251‑1 (MQ=255)                                                                                                                        gCCGACACCCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTGGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGca                                                                                                                  <  2:426918/251‑1 (MQ=255)                                                                                                                           gACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAg                                                                                                               <  1:493353/251‑1 (MQ=255)                                                                                                                                           aTGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGaa                                                                                              <  2:304087/251‑1 (MQ=255)                                                                                                                                                aaCGGTATGCGCTTTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACGACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAATCCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGATCT‑GGAACAGCa                                                                                         <  1:370548/251‑1 (MQ=255)                                                                                                                                                  cGGTATGCGCATTACCGCTCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCAcc                                                                                       <  1:249146/251‑1 (MQ=255)                                                                                                                                                  cGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTTTC‑GGTGATGAAACGCAATTT‑TG‑CCCGCAGG‑G‑TAACGCCCCG‑GGTGAAG‑TTCCTT‑TTCCCTACA‑ATTCAGCAGCCGAGCTGCCAAAACATTTTCA‑GGAAAACGGGCGGTCACCCTCTGGGCCGATGATGAAAAAAGAACAGGGGCTGCCCAGCAAAAAAAAGCT‑GGaaacgcacc                                                                                       >  2:501881/1‑244 (MQ=255)                                                                                                                                                   ggTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTTTC‑GTTGATGAAGAGCATTTT‑TG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAACCGATCTGCAAAAACATTGTCA‑GGAAAACGGGCTGTCAATCTCTGGCCTGATGATGAAAAAAGAGCTGGCGCTGCACAGCAAAGAAGAGGT‑GGAAACGCACCt                                                                                      >  1:126405/1‑251 (MQ=255)                                                                                                                                                      aTGCGCATTACCGCGTTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTGCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTgg                                                                                    <  2:348950/250‑1 (MQ=255)                                                                                                                                                         cgcATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCCTTTT‑GG‑CCAGCAGG‑G‑TAGCGCCCCG‑GGTGAAG‑TTCCTT‑ATCCGGACA‑GTTCAGCAGCCGATCTGCAAAAAAATTGTCA‑GGAAACCGGGCTGTCACTCTCTTGCCTGATGAAGAAAAAAGAGCTGGCGCTGCACAGCAAAGAAAAGCT‑GGAACCGCACCTGGCGa                                                                                 >  1:439347/1‑250 (MQ=255)                                                                                                                                                           cATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGG‑TTTATT‑GTTGATGAAGAGCATTTT‑TG‑CCCGCAGG‑A‑TAGAGCCCCG‑GTTGAAG‑TTCCTT‑TTCCGTACC‑GGTCAGCAGCCCATCTGCAAAAAACTTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGCTGAAAAACGAGCTGGCGCTGCACCGCAAAGAAAAGCT‑GGAACAGACCCTGGCGAACg                                                                              >  1:370824/1‑251 (MQ=255)                                                                                                                                                                      tGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGA‑TTTATT‑GTTGGTGGAGAGCATTTT‑TG‑CCAGCAGG‑G‑TAGCGCACCG‑GGTGAAG‑TTCCTT‑TTCCGTACA‑GTTCAGACGCCGATCTGCAAAAACCTTTTCA‑GGAAACCCGGCTGTCACTCTCTGGCCTGGTGATGGAAAAAGAGGTGGCGCTGCACAACAAAAAAGAGGT‑GGGACAGCACCCGGCGAACGTCTTGGAAGTc                                                                   >  2:349190/1‑251 (MQ=255)                                                                                                                                                                      tGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGCGGC‑TTTTTC‑GTTGATGGAGAGCATTTT‑GG‑CCAGCAGG‑G‑TAGCGCCCCG‑GCTGAAG‑TTCCTT‑ATCCGTCCA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAAACGGGCAGTCACTCTCTGGCCTGATGGTGAAAAACGAGCCGGGGCTGCACAGCAAAGAAGAGGT‑GGAACAGCACCCGGCGGACGTCTggggaggc                                                                   >  1:222841/1‑246 (MQ=255)                                                                                                                                                                                 aaaGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTAGC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCgcggcg                                                         <  2:5984/250‑1 (MQ=255)                                                                                                                                                                                   aGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGG‑TTTATC‑GTTGATGAAGAGCATTTT‑TG‑CCCGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑GTTCTT‑ATCCGGACA‑GTTTAGCAGCCGGTCTTCAAAAAAATTTTCA‑GGAAACCGGGCTGTCCATCTCTGGCCTTGTGGTGAAAAACAAGCTGGCGCCGGACAGCAAAAAAGAGCTTGGAA‑AGCAACCGGCGAACGGCTGGGaaagccagcgcggcgggg                                                      >  2:403969/1‑235 (MQ=255)                                                                                                                                                                                     tCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGG‑TTTTTC‑GTTGATGAAGAGCATTTT‑GG‑CCCGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑TTCCCTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTTGCCTGATGATGAAAAACGAGCTTGCGCTGCACAGCACAGAAGAGCT‑GGAACACCACCTGGCGAACCTCTGGGGAGTCATGCGCCGCGGTAtt                                                    >  1:161707/1‑251 (MQ=255)                                                                                                                                                                                     tCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑CTTATC‑GTTGCTTACTAGCATTTT‑TG‑CCAGCAGG‑A‑TAGGGCAACG‑GTTGGAG‑TTTCTT‑TTCCGTACA‑GGTCAGCAGCCGCTCTGCAAAAACATTGTCAAGGAAAC‑GGGCTGTCACTCTCTCGCCTGATGATGAAAAACGAGCTGGCGCTGCACCACAAAGAAGAGCT‑GGGAAAGCACCCGGCGAACGCCTGGGAAGGCAAGCGCGGCGGTAtt                                                    >  2:528949/1‑251 (MQ=255)                                                                                                                                                                                      cGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCTCCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTg                                                   <  1:216823/251‑1 (MQ=255)                                                                                                                                                                                                cAGACTTACTACTCTATTGGCGGTGGC‑TTTTTC‑GTTGATGAAGAGCATTTT‑TG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATTCGTACA‑GTTCAGCAGCCGATCTGCCAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGGCTGATGATGAAAAACGAGCTGGCGCTGCACAGCCAAGAAGAGCT‑GGAACAGCACCTGGGGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCgggat                                         >  1:38135/1‑248 (MQ=255)                                                                                                                                                                                                 ataCTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTAtt                                        <  2:403389/249‑1 (MQ=255)                                                                                                                                                                                                      gactactCTATTGGCGGTGGC‑TTTATC‑GTTGAGGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCAc                                   <  1:467432/250‑1 (MQ=255)                                                                                                                                                                                                       actactCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCAcc                                  <  2:439347/251‑1 (MQ=255)                                                                                                                                                                                                       actactCTATTGGCGGAGGC‑TTTATC‑GTTGATGAAGAGAATTTT‑TG‑CCCGAAGG‑A‑TAGCGGACCG‑GTTTGAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAAAATTGTCA‑GGGAACCGGGGTGTCCCTCTCTGGCCTTATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCCTGCGCGGGGGTATTGAGCGCGGGATTTTCAcc                                  >  2:106353/1‑251 (MQ=255)                                                                                                                                                                                                          actCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGTGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGaa                               <  1:389723/251‑1 (MQ=255)                                                                                                                                                                                                                    cGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATgcg                                                              <  1:118150/210‑1 (MQ=255)                                                                                                                                                                                                                    cGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATgcg                                                              >  2:118150/1‑210 (MQ=255)                                                                                                                                                                                                                        ggC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTgg                 >  2:78142/1‑251 (MQ=255)                                                                                                                                                                                                                          c‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGAGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGc                >  1:419315/1‑250 (MQ=255)                                                                                                                                                                                                                                agcggtgatgaAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTGGAAG‑GTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCGGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGGAGGGCT‑GGAACAGCACCTAGCGAACGGCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTTTTGCCTGGCAAACt           <  2:71224/246‑1 (MQ=255)                                                                                                                                                                                                                                      gatgaAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGt     >  1:291254/1‑251 (MQ=255)                                                                                                                                                                                                                                         gaAGAGCATTTT‑GG‑CCCGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTcc  >  1:408262/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        CCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTCAGCATGAAGTGGAGTTCCCGGTTGATCAGTGCATGAACTTCCACGCCGACAACCTTTCTCTGCATGAAAACGGTATGCGCATTACCGCGCTGGCGGGCGATAAAGTCGTTTACAGCCAGACTTACTACTCTATTGGCGGTGGC‑TTTATC‑GTTGATGAAGAGCATTTT‑GG‑CCAGCAGG‑A‑TAGCGCACCG‑GTTGAAG‑TTCCTT‑ATCCGTACA‑GTTCAGCAGCCGATCTGCAAAAACATTGTCA‑GGAAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCT‑GGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCC  >  NC_000913/2929805‑2930285

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑