Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	34.9%	70.6 / 11.7	43	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/27);  new base G (14/1);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 2.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAA  >  NC_000913/941466‑941958                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               ggTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGAc                                                                                                                                                                                                                                                             <  1:14576/251‑1 (MQ=255)      gCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGTTGACTTACt                                                                                                                                                                                                                                                        <  1:453250/251‑1 (MQ=255)            cggcACCATTACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTGTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑g                                                                                                                                                                                                                                                 <  1:545462/251‑1 (MQ=255)              gcACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑Gaa                                                                                                                                                                                                                                               <  1:550166/251‑1 (MQ=255)                          gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTGGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTT‑GATCAGGCACG‑CCAg                                                                                                                                                                                                                                  <  2:543618/251‑1 (MQ=255)                             ccTGGCCGCAGGGAGATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAAACCTTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGGCACAGCCCGTTGGATATCGAAAAACA‑TAAGCTGGTAGTGGTGGTTGGTAATAACCCTGGGGAAAAGCGCATGAGGGGCGGGTGGGGGGACTTataatttggtatgggacgccccaa                                                                                                                                                                                                                                   >  2:91148/1‑227 (MQ=255)                                 gccgccGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTa                                                                                                                                                                                                                           <  2:131275/251‑1 (MQ=255)                                            aTACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATg                                                                                                                                                                                                                <  1:232377/251‑1 (MQ=255)                                              aCCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGa                                                                                                                                                                                                               <  1:50553/250‑1 (MQ=255)                                               cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGatc                                                                                                                                                                                                             <  2:177376/251‑1 (MQ=255)                                                   ggTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACTGCGGCTGGGCAGATGGCAACAGCCCGTCGTATATCGAAAA‑CAGTAAGCTGGTAGTGCTGGTTGGTAATAACCCTGTCGAAACGCGAATGAGTTGCGT‑GGGGGTGTCTTGCTAGCTG‑GAATGGGCTCG‑CCAGAAATCTAATGCCCGCATGatcatc                                                                                                                                                                                                          <  1:533387/250‑1 (MQ=255)                                                   ggTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGatcatca                                                                                                                                                                                                         <  2:197262/251‑1 (MQ=255)                                                           ggCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                                                                                                                                 <  1:37550/251‑1 (MQ=255)                                                           ggCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCTATTTT‑TAAAAAGCACG‑CCCGAAATCTAATTCCCGCACGATCATTCTCGATccc                                                                                                                                                                                                 >  1:50918/1‑250 (MQ=255)                                                               gatgaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                                                                                                                               <  1:404985/249‑1 (MQ=255)                                                                  gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                                                                                                                          <  2:429797/251‑1 (MQ=255)                                                                          gcggcgGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                                                                                                                                    <  1:335270/249‑1 (MQ=255)                                                                             gcggcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTATCTT‑TAAAAAGCACG‑CCAGAAATCTTATGCCCGCATGATTATCATCGATTCCCGGCATACCGACACC‑ggg                                                                                                                                                                              >  1:90638/1‑250 (MQ=255)                                                                                gcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGcc                                                                                                                                                                           <  1:200733/251‑1 (MQ=255)                                                                                                cGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACGATACCGACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCAACACC‑GGTGCCGGGCGCGAAGATGAGt                                                                                                                                                           <  2:70828/251‑1 (MQ=255)                                                                                                   cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGCTATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGAGGCGG‑TGGGGGGGGTTACCATCTT‑TAACAGGGACGACCAGAAA‑ATTATGCCCCCATTAACATAACCTATCCGCGCCATACCGACACC‑CGGGCCGGGGGGGAAGATGGGTggg                                                                                                                                                        >  1:76313/1‑250 (MQ=255)                                                                                                           cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                                                                                                                                  <  2:195424/249‑1 (MQ=255)                                                                                                                 cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTT‑GAAAAAGCACC‑CCAGAAATTTAATGCCCGCCTGATCCTCATCGAACCCCGCTATAACGACACC‑CGGGCCGGGGGCGGAGATGAGGGGATCCCTCTTCGTCC‑g                                                                                                                                         >  1:154190/1‑251 (MQ=255)                                                                                                                         gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑g                                                                                                                                 <  2:581831/250‑1 (MQ=255)                                                                                                                             aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTT‑TAAAAAGCCCG‑CCAGGAAACTAATTCCCGCCTGAACATCCTCGAACCCCGCTATACCGACACC‑GGTGCCCGGCGCGGAGATGAGGGGATCCCTATTCGTCC‑GGGAACA‑AATGcc                                                                                                                            >  1:577716/1‑251 (MQ=255)                                                                                                                              aGGTTTGAACTATACCTGCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTTTCGGGCGCGAGGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCg                                                                                                                           <  1:34427/251‑1 (MQ=255)                                                                                                                                       cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTT‑GAACAGGCACG‑CCCGAAAACTAATGCCCGCCTTATCATCATCGATCCCCGCTATACCG‑CACCCGGGGCCGGGGGCGGAGAAGAGTGGATCCCCATTTGTCC‑GGGAAAA‑AATGCCCCAccgggtt                                                                                                                  >  1:495038/1‑245 (MQ=255)                                                                                                                                               aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTg                                                                                                       <  2:210432/251‑1 (MQ=255)                                                                                                                                                cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTTTAAAAAGGC‑CG‑CCCGAAAACTAAAGCCCGCATTATCACCATCGATACGCCCTACACCGACACC‑CGGGCCCGGGGCGAAGAAGAGGGGATCCCCATTCGTCCCGGGAA‑A‑AATGCCGCCCCT‑GG‑TTaagggctgg                                                                                                       >  1:41499/1‑244 (MQ=255)                                                                                                                                                                          cGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTTCTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑CGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑GTTACGGGCTGGCGGACGTCCTGATCACTGAAAACCtg                                                                            >  1:44558/1‑251 (MQ=255)                                                                                                                                                                           gTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCtgg                                                                           <  2:228710/251‑1 (MQ=255)                                                                                                                                                                           gTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTT‑TAACAGGCACG‑CCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGGGCGGAGATGAGTGGGTCCCCTTTCGTCC‑GGGGACA‑GATGCCGCAC‑T‑GG‑TTAAAGGGCTGGCGTACGTCATGATCACTGAAAACCtgg                                                                           >  1:417434/1‑251 (MQ=255)                                                                                                                                                                                             gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TTGGGGGGCTTCCCTTCTT‑TAAACGGCACG‑CCCGAAATCTTATTCCCGCACGATCATCCTCGAATCGCGCCACAACGACACC‑GGTGCCGGGCGCGGAGATGAGGGGATCCCCATTTGTCC‑GGGGACACAATGC‑GCAC‑TCGG‑TTAA‑GGGCTGGCGGACGTCATTATCACCGCAAAACTGGGGGATCAGGCATTtctc                                                          >  1:74663/1‑251 (MQ=255)                                                                                                                                                                                                 gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑a                                                     <  1:575217/251‑1 (MQ=255)                                                                                                                                                                                                  cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCTATCTT‑TAACAAGCACG‑CCCGAAAACTAATGCCCGCATGATCTTCCTCGATCCGCCCTATACCGACACC‑CGGGCCGGGCGCGGAGATGAGGGGGGCCCTTTTCGTCC‑GGGGACA‑GATGCCGCAC‑G‑GG‑GTTACGGTCTGGCGTACGGCAAGATCACCTAAAACCTGGGGGGTCAGGGATTCCTCGAT‑a                                                     >  2:96612/1‑250 (MQ=255)                                                                                                                                                                                                            tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTATCTT‑TAACAGGCACG‑ACCGAAATCTAATGCCCGCCAGATCATCATCGATTCGCCCCATACCGAAACC‑GGGGCCGGGGGCGGAGAAGAAGGGGGCCCTATTCGTCC‑GGGGAAA‑GATGACGCAC‑A‑GGCGTAA‑GGGCTGGCGTACGCCCTGATCACTGGACACCAGGGGGGTCAGGCATTTCCCGATTAAATAT‑GCGtt                                          >  2:435190/1‑251 (MQ=255)                                                                                                                                                                                                             gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTg                                         <  1:218982/251‑1 (MQ=255)                                                                                                                                                                                                                 ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTAATATCTT‑GAACAGGCACG‑CCCGAAATCTAATGCCCCCATGATCATCATAGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATACCTATTCGTCC‑GGGAACA‑GAAGCCGCCC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGCAAACCTGGTGGGTCAGGCATTCCTCGAT‑AAATATTTCGTTGGCt                                      >  1:483726/1‑250 (MQ=255)                                                                                                                                                                                                                       aaCCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAGCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCAGCGATTCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTTGATGCCTATTCGTCG‑GGGAACA‑TGTGGCGGAC‑G‑GG‑TTAACGGTCGGGCGTACGTCATGATCACGTAAGACCTGGTGGAGGAGGCATTCCTGGAT‑AAGTACTGCGTTGGTTACGATga                               <  2:260692/251‑1 (MQ=255)                                                                                                                                                                                                                                         aTGAGTGGCGG‑TGGGGTTACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTGGCGACGATGAGAAAACCCTGccagccag             <  2:254075/251‑1 (MQ=255)                                                                                                                                                                                                                                          tGAGTGGCGG‑GGGGGTGGATTACTATTTT‑TAAAAGGAAAC‑CCAGAAATATAATGTCCCCATTATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGGGGTGCCCTATTCTTCC‑GGGAAAA‑GAAGCCGCCC‑C‑GG‑TTAAAGGTCTGGCTTCCGGCATTATCACTGAAAACCTGGTGGAGCAGGGATTTCTCGAT‑AAAAATTGAGTTGGCTACGATGAGGAAACCCCGCCAACCAGt            >  2:429161/1‑251 (MQ=255)                                                                                                                                                                                                                                               ggCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGcc       <  1:203091/251‑1 (MQ=255)                                                                                                                                                                                                                                                      ggggTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGaaaa  <  2:445059/250‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               GGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAA‑CAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTT‑GAACAGGCACG‑CCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACA‑GATGCCGCAC‑T‑GG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAT‑AAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAA  >  NC_000913/941466‑941958

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