Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,757,061	C→A	100%	A251S (GCA→TCA)	fmt ←	10‑formyltetrahydrofolate:L‑methionyl‑ tRNA(fMet)N‑formyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,757,061	0	C	A	77.8%	10.6 / 2.7	9	A251S (GCA→TCA)	fmt	10‑formyltetrahydrofolate:L‑methionyl‑ tRNA(fMet)N‑formyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (0/7);  total (2/7)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGCTTTCCAGACTTTAACCGGCTGTCCTTCAATTTCCAGCCAGCTCATTGGCCATGGAT  >  minE/2757023‑2757117                                       |                                                         tCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:1332630/67‑1 (MQ=255) tCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:1480955/67‑1 (MQ=255) tCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:1758445/67‑1 (MQ=255) tCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:1968275/67‑1 (MQ=255) tCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:576953/67‑1 (MQ=255)  ccTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:64808/66‑1 (MQ=255)   cTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGATTTCCAGACTTTAACCGGCTGTCCTTCa                              <  1:342426/65‑1 (MQ=255)                              tGACCGATGCTTTCCAGACTTTAACCGGCTGTCCTTCAATTTCCAGCCAGCTCATTGGCCATGGAt  >  1:1641529/1‑66 (MQ=255)                              tGACCGATGCTTTCCAGACTTTAACCGGCTGTCCTTCAATTTCCAGCCAGCTCATTGGCCATGGAt  >  1:521139/1‑66 (MQ=255)                                       |                                                         TCCTGGTGCAGCGTTGGTTGCCGTATCAATGACCGATGCTTTCCAGACTTTAACCGGCTGTCCTTCAATTTCCAGCCAGCTCATTGGCCATGGAT  >  minE/2757023‑2757117

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