Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,471	0	T	G	54.5%	31.2 / 36.1	22	G331G (GGT→GGG)	ghxP	guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base T (0/10);  new base G (12/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.55e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05

GGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGA  >  NC_000913/4279336‑4279613                                                                                                                                        |                                                                                                                                               ggATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAAGGCGGGGGGGGGGGAAAACggggtttcg                                                                                                                                 >  1:393190/1‑144 (MQ=255)  gATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAGGGCGGGGGGCGGGGAAACCGGTTTGACgg                                                                                                                                >  1:140208/1‑151 (MQ=255)        gATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGGGGGGGGGAAAAACGGGTTGACGGCTATCa                                                                                                                          >  2:83542/1‑151 (MQ=255)        gATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGCGGGAAAACCGGgttttcgggtct                                                                                                                            >  1:126305/1‑140 (MQ=255)                 atcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTCGACGGCTATCACCGTTGGCg                                                                                                                 <  2:128089/151‑1 (MQ=255)                  tcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGAAAGGCGGCGGGGGGGTAAAACGGGTTGAGGGCTATTACCGTTGGCgg                                                                                                                >  2:81435/1‑150 (MQ=255)                         cGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTtc                                                                                                         <  2:310040/151‑1 (MQ=255)                             gggAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGGGGCGGGGGGGGAAAAAGGGTTGTAGGCTCTCCCCGTTGGCGTGCTGTtcctc                                                                                                      >  1:218432/1‑150 (MQ=255)                              ggAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGGGGGGAAAACGGGTTGACGGCTATCACCCGTGGCGGGGTGTttctc                                                                                                      >  1:145795/1‑149 (MQ=255)                                                  tCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGc                                                                                <  1:468027/151‑1 (MQ=255)                                                            gcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGGGGGGAAAACGGTTTTTCGGCTATCACCGGTGGCGGGCTGTTTCTCCTGAGTTTGTTTCTCTCTCCGCCCTCTTTcc                                                                       >  1:251441/1‑150 (MQ=255)                                                                      tGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCgggg                                                            <  2:136065/151‑1 (MQ=255)                                                                       ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGGGGGAAAACCGGGTTTACGGCTATCACCGTTGGCGTGCTGTTTCTCCTGGTTCTGTTCCTCTCTCCGCTCTCTTTCCTCGTTCCggggg                                                           >  2:186070/1‑150 (MQ=255)                                                                       ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGGGGGCGGGAAAAACGGGTTGGCGGCTATCACCGGTTGGGGGCTGTTTCTCCTGATTTTGTTTCTCTCTCCGCCCTCTTTCCTCGTTTCggggg                                                           >  2:46843/1‑150 (MQ=255)                                                                            ggtgggtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCt                                                                      <  1:442626/135‑1 (MQ=255)                                                                            ggtgggtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGGAAAAACGGGTTTACGGCTATCAACGTTGGCGGGGTGTTTCTCCTGATTCTGTTTCTCTCTCCGCTCTCTTTCCt                                                                      >  2:442636/1‑135 (MQ=255)                                                                                                  gTATACATCGAGTCTGCGGCGGGGAAGGCGGCGGGGGGGAAAAACCGGTTGAAGGGTTTTACCGTTGGCGGGCTTTTCCTCCTGGTTCTGTTCCTCTCTCCCCTCTCTTACCCCTGTCCGGGGGATGCAACGGCGCCGGCGCCGATTTACg                                >  2:525686/1‑151 (MQ=255)                                                                                                   tataCATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGGTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGt                               <  1:204047/151‑1 (MQ=255)                                                                                                       cATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTc                                               <  2:338289/131‑1 (MQ=255)                                                                                                                  cggcggGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatgc                 <  1:127466/150‑1 (MQ=255)                                                                                                                       ggTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGc            <  1:472548/150‑1 (MQ=255)                                                                                                                                 cgggcggTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGa  <  1:30455/150‑1 (MQ=255)                                                                                                                                        |                                                                                                                                               GGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGA  >  NC_000913/4279336‑4279613

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