Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	375,134	C→T	intergenic (+253/‑325)	mhpE → / → mhpT	4‑hydroxy‑2‑oxovalerate aldolase/3‑hydroxyphenylpropionate/3‑ hydroxycinnamate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	375,134	0	C	T	85.7%	10.7 / ‑5.0	7	intergenic (+253/‑325)	mhpE/mhpT	4‑hydroxy‑2‑oxovalerate aldolase/3‑hydroxyphenylpropionate/3‑ hydroxycinnamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (2/4);  total (2/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCC  >  NC_000913/374938‑375358                                                                                                                                                                                                     |                                                                                                                                                                                                                                 cggcgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGTCATCTGCGCACGACGGCCGGATGCGGCGTGAACGCCTTATCCGGGCGACGGATGGCGCGGGAATTTGTAGTCCTGATAAGACCCGTTTGCGTCGCATCAGGCATCTGCGCACGACTGCCTTATGCGGCGTGAACGTCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc                                                                                                                                                                            <  1:162185/251‑1 (MQ=14)                                     agaATTTGGAGGTCTGATAAGACGCGTTAGCGTCGCAACAGGCATCTGCGCACGACGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTTCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgc                                                                                                                                        <  1:37858/251‑1 (MQ=255)                                       aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgcgc                                                                                                                                      <  1:22072/251‑1 (MQ=255)                                               ggTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGc                                                                                                                              >  1:52634/1‑251 (MQ=255)                                                                tAGCGTCGCGTCAGGCATCTGCGCACGACTGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGTCGTAAAc                                                                                                              <  1:130393/250‑1 (MQ=255)                                                                                                                                                                 ncgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAGCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCACACGACTGCCGGATGCGAtca             >  1:101421/2‑248 (MQ=255)                                                                                                                                                                           cGCATCAGGCATCTGCGCACGACTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTcc  <  1:65581/251‑1 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                                                                 CGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCC  >  NC_000913/374938‑375358

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