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

GGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGG  >  NC_000913/941466‑941962                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           ggTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAc                                                                                                                                                                                                                                                         >  1:315385/1‑251 (MQ=255) ggTACGCGGGGCGGACCCACGACCCGCTCATGGCCGCCGGGAAAGTCCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAAGTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGGTGGGAGTGCTGTTTGGTAATATCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAc                                                                                                                                                                                                                                                         <  2:369115/251‑1 (MQ=255)  gTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                                                                                                                        <  1:55374/251‑1 (MQ=255)    aCGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                                                                                                                                       <  1:35553/250‑1 (MQ=255)                   aTGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGGGGCGGGGGGGGGGATTACTATCTTGAACaagg                                                                                                                                                                                                                                       >  1:53525/1‑248 (MQ=255)                        ccGCGCCTGGCCGCAGGGAAAGACCCGGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAGTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                                                                                                                                  <  2:103588/251‑1 (MQ=255)                                ggCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAAGAACCCTGGCGAAACGTGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                                                                                                                            <  1:376628/249‑1 (MQ=255)                                                      cgcgcgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                                                                                                                     <  2:99649/250‑1 (MQ=255)                                                         gcgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                                                                                                                                 <  1:216851/251‑1 (MQ=255)                                                               gatgaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACGCCGCCGCGCAAATTGCGGAAGGTTTGAACTACACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATTAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCTATCCGCGCt                                                                                                                                                                                           <  2:449559/251‑1 (MQ=255)                                                                          gcggcgGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                                                                                                                                <  1:379220/251‑1 (MQ=255)                                                                               ggcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGc                                                                                                                                                                          <  1:221650/251‑1 (MQ=255)                                                                                gcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGcc                                                                                                                                                                         <  1:395053/251‑1 (MQ=255)                                                                                                cGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGt                                                                                                                                                         <  1:328459/251‑1 (MQ=255)                                                                                                    gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                                                                                                                                  <  2:188300/207‑1 (MQ=255)                                                                                                           cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                                                                                                                              <  1:137154/251‑1 (MQ=255)                                                                                                               cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                                                                                                                          <  1:253521/251‑1 (MQ=255)                                                                                                                      aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTCTATATTTAACAGGGACGCCAGAAAACTAAAGCCCCCCAGATCATCCTCGATCCGCGCTATAACGACACC‑GGGGCCGGGGGCGGAGATGAGTGGAGCCCTCTTCGTCCgggga                                                                                                                                    >  1:430451/1‑248 (MQ=255)                                                                                                                                      acgataCCTACGGCGGCTGGGCAGAGGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGGGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGGCCGGGAACAGATGCCGCACTGGtt                                                                                                                   <  2:172798/248‑1 (MQ=255)                                                                                                                                        gagtccTACGTCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAGGAGTTGCGGTGGGGTGACTTACTATCTTGAACAGGCACGGCAGAAATCTAAGGCCCGCATGATCATCATCGATCCGCGCTATACCGACTCC‑GTTGCCGGGCGCGAAGTGGAGTGGATCCCTATTCGTCCGGGAACAGAGGCCGCACTGGTTa                                                                                                                  <  2:504958/246‑1 (MQ=255)                                                                                                                                                    ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTg                                                                                                                                                                           <  1:406103/181‑1 (MQ=255)                                                                                                                                                    ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAATATTAAGCCCGCCAGATTCTCATCCATCCCGGCTCAACCG‑CACCCGGTGCCGGGGGCGAAGGTGAATGGGTCCCCCTTCGTCCGGGGACACATGCCGGCCTGGTTTACGGGCTgggggt                                                                                                     >  2:237826/1‑247 (MQ=255)                                                                                                                                                        tggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTGCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATTATCATCATCGATCCGCGCTATGCCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGAGCCCTATTCGTCCGGGAACGGAGGCCGCACTGGTTAACGGTCTTGCGTACgg                                                                                                  <  2:480632/249‑2 (MQ=255)                                                                                                                                                          gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAAGGCGAATGAGTGGCGGTGGGGGGTCTGACTATCTTGAACAGGCACGCCAGAAATATAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAGACC‑GGTGCCGGGCGCGAAGATGAGTGGATCGCTATTCGTCCGGGAACAGATGCCGCACTGGTGAACTGGCTTGCGTTCGTCaa                                                                                               <  2:118169/251‑2 (MQ=255)                                                                                                                                                                ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCa                                                                                          <  1:431157/250‑1 (MQ=255)                                                                                                                                                                      aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATTCCCCCATGATCATCCTCGATCCGCGCTCTAAAAACACC‑GGTGCCGGGGGGGAAGATGAGTGGATACCTATTCGTCCGGGGACAGATGCCGCACTGGGTAAAGGTCTGGGGTaaggcctggtcaccgaaa                                                                                    >  1:53917/1‑233 (MQ=255)                                                                                                                                                                         ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCt                                                                                <  1:67283/251‑1 (MQ=255)                                                                                                                                                                          cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCCCATGATCCTCATCGATTCGCGCTATAACGACACC‑GGGGCCGGGCGCGGAGATGAGGGGGTCCCTCTTCGTCCGGGAACAGATGCCGCACAGGTTAACGGGCTGGCGTACGGCCTGATTCCCGAAAACCtg                                                                               >  1:21062/1‑251 (MQ=255)                                                                                                                                                                             cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGACCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCtggtg                                                                            <  2:332124/251‑1 (MQ=255)                                                                                                                                                                                     gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTAACTTTAAAAAGCCCCCCAGAAAACTAAAGCCCCCATGAACATCATCCATCCCCCCTCTACCGACACC‑CGGGCCGGGCGCGGAGAAGAGGGGGACCCTATTCGTCCGGGGACAGAGGCCGCCCCGGGTAACGGTCGGGCGGACGTCCTGATCACTGAAAACCCGGGGGATccagc                                                                    >  2:138882/1‑247 (MQ=255)                                                                                                                                                                                      aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGCAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGGGCCGGGGGCGAAGATGAGGGGAATCCTATTCGTCCGGGGAAAGAAGGCGCACTGGTTAACGGTCTGGCGGACCTCATGATCTCCGAAAAACTGGGGGATCaagc                                                                    >  1:466009/1‑247 (MQ=255)                                                                                                                                                                                         aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTc                                                                <  1:49407/251‑1 (MQ=255)                                                                                                                                                                                           aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCt                                                              <  1:144841/251‑1 (MQ=255)                                                                                                                                                                                              aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCg                                                            <  1:199085/250‑1 (MQ=255)                                                                                                                                                                                              aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAAAAGGCACGCCAGAAAACTAATGGCCGCCTGAACATCATCGATCCGCGCTATACCGACACC‑GGGGCCGGGGGCGGAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCCCTGGGTTACGGGCTGGGGTACGTCATGATCACTGAAAAACTGGGGGGTCAGGCaattctcgg                                                           >  1:850/1‑243 (MQ=255)                                                                                                                                                                                               aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGAt                                                          <  1:267778/251‑1 (MQ=255)                                                                                                                                                                                                   ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCAATAAat                                                      <  2:233909/251‑1 (MQ=255)                                                                                                                                                                                                         gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCg                                                <  1:10420/251‑1 (MQ=255)                                                                                                                                                                                                              ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTgg                                            <  1:401213/250‑1 (MQ=255)                                                                                                                                                                                                                     taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGAACCGCGCTATACCGACACC‑CGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCACTGGTTTACGGTCTGGCGTACGTCCTGGACACTGAAAAACTGGTGGGTCAGGCATTCCTCCGTAAATATTGCGTTGGGTACGATg                                    >  1:30346/1‑251 (MQ=255)                                                                                                                                                                                                                                      gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgccagcc                   <  1:73101/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgccagcca                  <  2:66119/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgccag                     <  1:257397/248‑1 (MQ=255)                                                                                                                                                                                                                                               gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCg          <  1:405168/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ggggTGACTTGCTATATTGGACAGGCACGCCAGAAATCTAATTCCCGCATGATCATCATCGATCCGCGCTATGCCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTTTGGCGTACGTCATGATCACTGAAAACCTTGTGGATCAGGCTTTCCGCGATAAATATTGCGTTGGCTACGATGAGAAATCCCTGCCAGCCAGTGCGCCGaaaaa     <  2:437004/251‑1 (MQ=255)                                                                                                                                                                                                                                                       gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATgg  <  1:152017/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           GGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACC‑GGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGG  >  NC_000913/941466‑941962

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