Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,072,731	0	T	G	55.6%	13.9 / 20.7	18	T45P (ACC→CCC)	rutB	ureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (10/0);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 6.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04

TTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAG  >  NC_000913/1072586‑1072879                                                                                                                                                  |                                                                                                                                                     ttCGATTTATGAAAATCCGGTGAGCCGGGTCCGCCAGCCGCGACATACTGTTCATCCCAGCCATTTTGAACCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTgag                                                                                                                                                  <  2:350549/150‑1 (MQ=255)  tCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGAc                                                                                                                                                <  1:54071/151‑1 (MQ=255)     aTTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGTTGGCAATGACCGGGGGGGTGGGTGAGACAt                                                                                                                                              >  2:279284/1‑150 (MQ=255)     aTTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGGTTGTATGTTGGCAATGACCGGGCGAGTGGGTGGGACAt                                                                                                                                              >  1:299477/1‑150 (MQ=255)                                          gACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGATTTAATGGTTGCAACTACCCGGCGAGGGGGTGAGACATCAAAACCGGCGAGATCTAAGGAGCCGCCTGGCGGGGCa                                                                                                        >  1:3918/1‑151 (MQ=255)                                           aCATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCCATGACCGGGCGGGGGGGTGAGACATCAAACCCGGCGGGGTCTAAGGAGCCGCCTGGCGTGGCAt                                                                                                       >  2:245051/1‑151 (MQ=255)                                             aTACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGGAATTACCGGGCGGGGGGGGGAGGCACCAAAACCGGCGCGATCTAAGGAGGCGCCTGGCGGGGGATaa                                                                                                     >  2:138161/1‑151 (MQ=255)                                             aTACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCCGGCGCGGGGGTGAGACATCAAACCCGGCGGGATCTAAGTAGCCGCCTGGCGTGGCATaa                                                                                                     >  2:335543/1‑151 (MQ=255)                                                cTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGGAATGACCGGGCGAGTGGTTGGGACATCAAACCCCGCGAGATCTAAGTAGCCGCCCGGCGTGGCGTAAGCg                                                                                                  >  1:29272/1‑151 (MQ=255)                                                cTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGc                                                                                                   <  1:111174/150‑1 (MQ=255)                                                cTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCCCGCGCTGCGGTCACGGCGGTTTGAACAATGGCAATTACCCGGGGGGTGGGTGGGAGAACAAAACCCGGGGGATCTAAGGACCCGGCTGGGGGGGCATAAgg                                                                                                   >  2:111194/1‑149 (MQ=255)                                                                 tttGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGGGGGGGTGAGAAACCACAACCGGGGGGGTCTAAGTAGCCGCCCGGGGGGGCATAAGCGTTTTTCATATacaccc                                                                                  >  2:116935/1‑149 (MQ=255)                                                                 tttGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGGTTGGATGTTGGGAATGACCGGGGGGGGGGGTGGGGCATCAAACCCGGCGGGGGTTAAGAAGCCGCCCGGGGGGGGATAAGGGTTTTGCAATTccacca                                                                                  >  2:7470/1‑150 (MQ=255)                                                                                      tCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGGTGGCAAAGACCGGGCGCGTGGGTGGGACATCAAAACCCGCGAGAGTTAAGTAGCCGCCCGGCGGGGGATAAACGGTTTGTATATCCACCACGATCAGGGCACCTTGCTGCgg                                                            >  1:320758/1‑151 (MQ=255)                                                                                                     cTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATgg                                              <  1:23061/150‑1 (MQ=255)                                                                                                                              ggCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtcgt                     >  1:8836/1‑150 (MQ=255)                                                                                                                                         ggCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATCATGCCAc          <  1:369758/150‑1 (MQ=255)                                                                                                                                                ggTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGCCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAg  >  1:409821/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAG  >  NC_000913/1072586‑1072879

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