Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	375,067	A→G	intergenic (+186/‑392)	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,067	0	A	G	93.8%	26.0 / ‑3.1	16	intergenic (+186/‑392)	mhpE/mhpT	4‑hydroxy‑2‑oxovalerate aldolase/3‑hydroxyphenylpropionate/3‑ hydroxycinnamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (4/11);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCG  >  NC_000913/374969‑375074                                                                                                   |        gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTg       <  1:500055/101‑1 (MQ=25) gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTg       >  1:500520/1‑101 (MQ=38)  cGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTgg      >  1:500521/1‑101 (MQ=25)  cGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTgg      <  2:500068/101‑1 (MQ=25)   gTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgc     <  1:500063/101‑1 (MQ=18)   gTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgc     <  2:500120/101‑1 (MQ=18)   gTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgc     >  2:500523/1‑101 (MQ=25)    tGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcg    <  1:500087/101‑1 (MQ=17)    tGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcg    <  2:500206/101‑1 (MQ=17)     gagaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcgc   <  1:499776/101‑1 (MQ=14)     gagaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcgc   <  1:500100/101‑1 (MQ=14)     gagaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcgc   <  2:500043/101‑1 (MQ=14)     gagaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGCCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcgc   <  2:499932/101‑1 (MQ=9)      agaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGcgcg  <  1:500020/101‑1 (MQ=11)      agaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGcgcg  >  1:500525/1‑101 (MQ=11)                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGcgcg  >  2:500463/1‑43 (MQ=35)                                                                                                   |        GCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCG  >  NC_000913/374969‑375074

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