Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,892,019:1	+GC	100%	intergenic (‑437/+54)	gudP ← / ← ygdL	predicted D‑glucarate transporter/conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,892,019	1	.	G	100.0%	16.1 / NA	7	intergenic (‑437/+54)	gudP/ygdL	predicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (3/4);  total (3/4)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*	minE	1,892,019	2	.	C	100.0%	16.0 / NA	7	intergenic (‑437/+54)	gudP/ygdL	predicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (3/4);  total (3/4)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTC  >  minE/1891959‑1892055                                                              ||                                     tGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagga                                     >  1:4545361/1‑61 (MQ=255)                       cACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt        <  1:1057122/71‑1 (MQ=255)                           gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg    >  1:2106524/1‑71 (MQ=255)                               aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtc  <  1:2170658/69‑1 (MQ=255)                                                  ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAg             <  1:2399963/39‑1 (MQ=25)                                                  ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAg             <  1:4890514/39‑1 (MQ=25)                                                    gAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc           >  1:391254/1‑39 (MQ=25)                                                    gAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc           >  1:4546749/1‑39 (MQ=25)                                                              ||                                     TGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTC  >  minE/1891959‑1892055

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