Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,532,228	A→T	100%	D166E (GAT→GAA)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,532,228	0	A	T	95.3%	43.4 / ‑4.1	21	D166E (GAT→GAA)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (0/20);  total (1/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCG  >  minE/1532206‑1532288                       |                                                             cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2528438/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:461837/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:384659/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:3433037/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:3195152/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2996289/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2994602/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2938176/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:277767/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2646980/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:1098693/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2307573/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:1796947/69‑1 (MQ=255) cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:1785016/69‑1 (MQ=255)  gACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa                <  1:2470535/68‑1 (MQ=255)          aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt        <  1:1585249/68‑1 (MQ=255)          aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt        <  1:1511539/68‑1 (MQ=255)          aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt        <  1:3226501/68‑1 (MQ=255)          aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt        <  1:1375292/68‑1 (MQ=255)          aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt        <  1:1358441/68‑1 (MQ=255)                   tGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCg  >  1:1286996/1‑65 (MQ=255)                       |                                                             CGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCG  >  minE/1532206‑1532288

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