Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,111,527	G→A	100%	C287Y (TGT→TAT)	ydiD →	short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,111,527	0	G	A	76.3%	63.8 / 18.3	38	C287Y (TGT→TAT)	ydiD	short chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (8/1);  new base A (0/29);  total (8/30)
Fisher's exact test for biased strand distribution p-value = 1.84e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACC  >  minE/1111486‑1111558                                          |                                cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcagcagc          >  1:1806352/1‑65 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:30995/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:1313223/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:1346415/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:679228/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:566123/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt        >  1:526132/1‑67 (MQ=255) cgcTCGCAGCGTGTTGTTAGATATTTTAACTCCTGATGCGTGTCTCGCgct                        >  1:184187/1‑51 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1234013/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:2093137/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:2025215/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:213810/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:2309922/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:346221/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:443517/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:462421/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:569891/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:65676/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:818460/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:849016/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:918556/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:978350/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:2088754/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:2069413/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:106114/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1929755/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1905820/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1881278/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1775782/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1701889/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1507063/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1457170/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1446174/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1258597/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1230224/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1189112/66‑1 (MQ=255)        aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1154680/66‑1 (MQ=255)         gCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCAcc  <  1:1362311/65‑1 (MQ=255)                                          |                                CGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACC  >  minE/1111486‑1111558

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