Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,979	0	G	A	72.5%	59.0 / 21.5	40	V71M (GTG→ATG)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (2/9);  new base A (8/21);  total (10/30)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCC  >  NC_000913/339681‑340260                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                          cATCCAGCTTTTGCCTGCGTTGAGAGAGCTAGTTTTGGCTGGGTCGAATGTTACGGAGTCTGGCCTACATGTTCATGCTGGGTGCGGCGTGAACGTCTTATCCTGCCTACGAATGGTGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCATTAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgc                                                                                                                                                                                                                                                                                          <  2:719326/300‑2 (MQ=14)                    tCGACAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGATCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGcc                                                                                                                                                                                                                                                                      <  1:242384/300‑1 (MQ=255)                         agagctaTTCGTGGCTGGGGCGAACGTTACGGCGTCTGGCCTACATGTTCATGCNGGATGCGGCGTGAACGCCTTATCTGGCGTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATa                                                                                                                                                                                                                                                                 <  2:625840/294‑1 (MQ=14)                              cattttttattGGTGCGAACGTTACGGCGTCTGACTTGCATGTTCATGCCGGATTCGGCGTGAACGCCTTATCCGGCCTATGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTTACCATTGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGt                                                                                                                                                                                                                                                            <  2:453338/290‑1 (MQ=14)                                   ttggctagtgcgagcgttgctttgTCTGACCTACATGTTCAGGCCGGATGCGGCTTGAACACCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGAGAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAGTGGTGCGAGAAT‑TGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                                                       <  2:163703/277‑1 (MQ=11)                                   ttGGCTGGGGCGGACGTTATGGCGTCTGACCTACCTGTTTGTGCCGGATGCGGCATGAACGCCTTATTCGGCCTACGAGTGGCGCAAGAGGCTGTAGGCCTGATAAGCGTAGCGCATCATGCATTTGTTACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAACCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTCATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAGGCGTAgcgc                                                                                                                                                                                                                                                       <  2:177788/300‑1 (MQ=11)                                              gAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGCTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCAGCGTGAACGCCTTATCCGGCATACGCGTGGCGAGAGAAT‑TGTAGGCCTGATAAGCGTAGAGCATCAggtagtc                                                                                                                                                                                                                                            >  2:638883/1‑295 (MQ=14)                                                   ttACGGCATCTGACCTACATGTTCATGCTGGAGGCGTCGTGAACGCCTTAGCCGGCCTACGAACGGCGCAAGAATCTGGAGGCCCGATAAGCGTAGCGCATCAGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGGATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcg                                                                                                                                                                                                                                       <  2:109706/300‑9 (MQ=11)                                                    tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCTCAAGAATCTGTAGGCGTGATAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                                                      <  2:509667/300‑10 (MQ=11)                                                    tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCNTACGAATGGCGCAAGAATCTGTAGGCNTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                                                      <  2:569264/300‑10 (MQ=11)                                                    tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                                                      <  1:22707/300‑10 (MQ=255)                                                         cGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagcctggcgttggtc                                                                                                                                                                                                                                 >  1:139529/1‑286 (MQ=11)                                                                                                                                                           gacctaaatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCGGGCATTTGTCGCCATTGCCGGATGCGGTGTGAACGTCTTATCCGGCCTCCGAATGGCTAGAGAATTGTAAGCTTGATACGCGCCGCGCCCCAGGCAGTCTGGCGTTGGTCATagc                                                                                                                                >  2:19548/15‑298 (MQ=11)                                                                                                                                                            acctgcgtgttcatGCCGGATGCGTCGTGAACACCTTATCCTGCCTGCGAATGCNGCAAGAATCTTTATGCCTGATAAGCGTTACGCATCGGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGAATGCGGCGTGAACGCCTTACCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                                               <  2:635013/287‑1 (MQ=11)                                                                                                                                                            acctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                                               <  1:165262/287‑1 (MQ=11)                                                                                                                                                            accgacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                                               <  1:321387/287‑1 (MQ=11)                                                                                                                                                             cctgcatgttcatGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCTTAGCACATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGGTGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                                                                              <  2:44894/288‑1 (MQ=11)                                                                                                                                                                          gcctgATGCTGCGTGAACGCCTTATCCGGCCTACGAACGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCGGATAAGCGTGGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                                                                 <  1:635894/296‑1 (MQ=14)                                                                                                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATGAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                                                                 <  1:798161/300‑1 (MQ=255)                                                                                                                                                                            cGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                                                               <  1:264644/300‑1 (MQ=255)                                                                                                                                                                            cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGCGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                                                               <  1:266323/300‑1 (MQ=255)                                                                                                                                                                               attcggcgTGAACGCCCTATCCGGCTTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                                            <  2:434287/297‑1 (MQ=14)                                                                                                                                                                               aTGCGGCGTGAACGCCTTATCTTGCCTACGAATTGCGCAAGCATCTGTAGGCCTAATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTACCGGATGCGGCGTGAACGCCTTATCCGGCCTACGATTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                                            <  2:465861/300‑1 (MQ=14)                                                                                                                                                                               aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAATAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                                            <  1:560656/300‑1 (MQ=255)                                                                                                                                                                                   ggcgTGAACACCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCAACCTACAAATGGCGCAAGAATCTGTAGGCCTAATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGACATACGAATGGCGCGAGGGTTGTAGGCCTGATAAGCGCAGCGCACCAGACAGTCTGACGTCGGTCATAACCACCTCACCCTCTACAGCGGatat                                                                                                        >  2:303982/1‑297 (MQ=11)                                                                                                                                                                                      gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaac                                                                                                     <  1:453483/300‑1 (MQ=255)                                                                                                                                                                                           cGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                                                                                                >  1:151564/1‑300 (MQ=255)                                                                                                                                                                                                gattcGGCCTACGAACGGCGCAAGAAGCTGTAGGCCTGATAAGCGTAGCGAAGCAGGCATTTGGCAACAGTGCCGGAGGCGGCGGGAACGCCTTATCCGGCCTACGAATNGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCTGACGAGTGGTGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAag                                                                                           <  1:600291/296‑1 (MQ=17)                                                                                                                                                                                                      gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTCGGCCTGCTAAGCGTAGCGCATCAGGCCGTCTGGCGTTGGTCATAACCCCATCACCCCCTGTAGCGGACATAACAACATGAAAGAGCTtg                                                                                     >  2:237954/1‑300 (MQ=17)                                                                                                                                                                                                                 gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGANTTGTCGGCCTGATACGCTTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGAACATAACAACATGAAAGAGCGTGTGGTCGTTGcc                                                                          >  2:727541/1‑300 (MQ=17)                                                                                                                                                                                                                  cacAAGNATCTGTATGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGNCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCa                                                                         <  2:695636/298‑1 (MQ=255)                                                                                                                                                                                                                   gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCAt                                                                        >  1:734371/1‑300 (MQ=255)                                                                                                                                                                                                                          tctgtggacttgatgaGCGTGGCGCATCCTGCATTTGGCACTCGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGc                                                                 <  2:358214/290‑1 (MQ=18)                                                                                                                                                                                                                             gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAAc                                                              <  2:266863/300‑1 (MQ=255)                                                                                                                                                                                                                                          gCGTAGCGCATCAGGCATTTGTAACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCCGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAGTGGCGCGAGAGTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATCTGGCGTTGGTCACCACCCAATCCACCTCTGTAGAGGAAACAACAACATGAAAGAGCCCGTGGTCGTTGACATTGGTGGCGACAGCATCATCAAAg                                                 >  2:11973/1‑300 (MQ=14)                                                                                                                                                                                                                                            gTAGCGCATCAGGCGTTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGNCCTACGAGTGGCGCGAGAATTGTAGGNCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGAt                                               <  2:573573/300‑1 (MQ=255)                                                                                                                                                                                                                                                      aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCNGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCNTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgccagcc                                     <  2:659969/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    attGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCAt                       <  1:437799/298‑1 (MQ=255)                                                                                                                                                                                                                                                                         cGGATGCTGCGTGAACGCCGTATCCGGCCGACGAATGGCGCAAGAATCTGTAGGACTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCNCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCaggc                  <  2:187920/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                     gaaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGGGaa      <  1:696734/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                         gCCTTATCCGGCCTACGAATGGCGCAAGAATNTGTAGGCCTGATAAGCGTAGCGCACCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGAGTGGCGCGGGGGTNGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCACAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGAGGTCAAAGcc  >  2:581694/1‑300 (MQ=25)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                          CATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCC  >  NC_000913/339681‑340260

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