Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	111,195	0	T	G	52.9%	4.3 / 11.1	17	V51G (GTG→GGG)	mutT	8‑oxo‑dGTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (9/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTG‑AACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGA  >  NC_000913/110947‑111439                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGc                                                                                                                                                                                                                                                     <  2:205043/251‑1 (MQ=255)       aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTG‑AAc                                                                                                                                                                                                                                              <  1:117853/251‑1 (MQ=255)               cAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGGGGGGCGGG‑AACTTCAGga                                                                                                                                                                                                                                       >  1:144021/1‑250 (MQ=255)                                      aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACATGCGGTGGTGCGTG‑AACTTCAGGAAGAAGTCGGGATTACCCCCCAACa                                                                                                                                                                                                               <  1:19263/251‑1 (MQ=255)                                                      gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGCAACAGCTGCAAATTGCGGTAGGCATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGCGGGGGGGG‑GATATCAGGGAGAAGGCGGGGTTACCCCCCAACAttttttgtcttttta                                                                                                                                                                                                >  1:160069/1‑240 (MQ=255)                                                                          cgcCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGGGGGGCGTG‑TACTTTAGGGAGAAGACGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGGATATGAATTccc                                                                                                                                                                           >  1:37854/1‑251 (MQ=255)                                                                                  ttATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTG‑AACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAgg                                                                                                                                                                    <  1:137775/250‑1 (MQ=255)                                                                                                cAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCCGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGTG‑AACTTTTAGAGGAAGACGGGGGTACCCCCCCACCTTTTTTTAACTTTGAAAAACACGGATATGAATTCCCGTACAGGCATATTAAACTGTgg                                                                                                                                                     >  1:246008/1‑251 (MQ=255)                                                                                                cAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGTG‑AACTTCCGGAAGAAGTCGGGATTACCCCCCCACATTTTTCTCTATTTGAAAAACTGGGATATTAATTCCCGGACAGGCATATAAAACTGTgg                                                                                                                                                     >  1:167779/1‑251 (MQ=255)                                                                                                           cTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGGGAAACGCCGGAACAGGCGGTGGTGCGTG‑AACTTCCGGAAGAAGTCGGGATTACCCCCCAAAATTTTTCTCTTTTTTGAAAACTGGAATAAGAAATCCCCGACCGGCATATAacaactttgtttttggttgt                                                                                                                                          >  2:110205/1‑234 (MQ=255)                                                                                                                              ttattCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGAGTG‑AAATTCAGGAAAAAGTCGGGGGTACCCCCCAAAATTTTTCGCTATTTTAAAAAATGGGATATTAATTCCCGGGCACGCATATAAAACTTGGGTTTTGGGTgggtggaagatggggagggggg                                                                                                                       >  1:221035/1‑231 (MQ=255)                                                                                                                                                           ataACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGGG‑AAATTTAGGGAGAAATCGGGGTTTCCCCCCAACATTTTTTGCCTTTTGAAAAACTGGGAAATGAATTCCCGGACAGGCCTATAACACTGTGGTTTTTGCTGGGCGGACGCTGGGGAGGGGGGGCGTGGGGTGAAGAAGGGCAACCCggggg                                                                                          >  2:239185/1‑248 (MQ=255)                                                                                                                                                            taACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTTAAATGGGTGAAACGCCGGAAAAGGCGGTGGTGCGGG‑AACTTCAAGAAGAAGTCGGGGTTACCCCCCCACATTTTTCGCTCTTTTAAAAACCGGAGAATGAATTTCCCGCCAGGCATTTAACACCGTGGTGTTGGCTGGTCGAAAGCTGGGAAGGGGGGCCCAGGGGTAAAGAAGGGGAAACCgggt                                                                                           >  2:70361/1‑247 (MQ=255)                                                                                                                                                              aCGCGTCGCGCAGCAGATGCGCACATGGCGAATACACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGTGGAATTTC‑GGGAGGAGACGGGGTTTCCCCCCAACCTTTTTCGCTATTTTAAAAAATGGGATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTTGGTGGGTGAACGCTGGGAAGGGGAGGCGGGGGGGAAAGAAGGGCAACCCGGTGAGTg                                                                                       >  2:108000/1‑251 (MQ=255)                                                                                                                                                                           cagATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGGGTG‑AACTTCAGGAAGAAATCGGGGTTACCCCCCCACATTTTTCGCTATTTGAAAAAATGGAATATGAATTCCCGGACAAGCCTATAACACCGTGGTTTTGGCTGGTCGAACGATGGGAAGGGGGGCCGGGGGGGAAAGAAGGGCAACCCGGGGAGAGGGTGTCggcgggg                                                                          >  1:11613/1‑245 (MQ=255)                                                                                                                                                                                                               gcggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTG‑AACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTTGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGaa                                      <  2:244596/251‑1 (MQ=255)                                                                                                                                                                                                                                                  gCGGTGGTCCTTA‑AACTTCATGAAGAATTCGGGATTACCCCCAAACATTTTTCGCTATTTGAAAAACGGGAAATTTAATCCCCGGACAGGCATTTACCATTTTGGTTTTGGCTGGTCAACCGTTGGGAAGGGGAGCC‑TGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGa  <  1:236714/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      GTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTG‑AACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGA  >  NC_000913/110947‑111439

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