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%	22.9 / NA	8	intergenic (‑437/+54)	gudP/ygdL	predicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (2/6);  total (2/6)
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%	23.2 / NA	8	intergenic (‑437/+54)	gudP/ygdL	predicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (2/6);  total (2/6)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCG  >  minE/1891951‑1892062                                                                      ||                                            cTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCag                                              <  1:234570/70‑2 (MQ=255)               aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                               >  1:28382/1‑71 (MQ=255)                           tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                              <  1:209117/60‑1 (MQ=255)                                aCCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt               <  1:34176/70‑1 (MQ=255)                                  cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                   <  1:328245/48‑1 (MQ=255)                                    cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg           >  1:44123/1‑70 (MQ=255)                                             cATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg  <  1:173470/70‑1 (MQ=255)                                              aTTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                <  1:69556/39‑1 (MQ=37)                                               tttACCGCTCGCTGAACATATCAGCTAAGCACAGa                                   <  1:313393/35‑1 (MQ=25)                                                                      ||                                            CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCG  >  minE/1891951‑1892062

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: ‑