Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,593,534	G→T	I113I (ATC→ATA)	cysA ←	sulfate/thiosulfate transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,593,534	0	G	T	100.0%	127.9 / NA	40	I113I (ATC→ATA)	cysA	sulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (0/40);  total (0/40)

CACTTTCGCTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCA  >  minE/1593522‑1593582             |                                                 cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGTCAGCACCGTCAGGCCa  <  1:2683622/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:526953/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:3265042/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:3272395/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:3280672/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:331871/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:364770/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:42153/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:485620/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:497391/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:3098301/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:543503/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:586025/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:688798/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:712963/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:750087/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:750678/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:785792/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:858550/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:884948/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:3083481/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1561841/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1576558/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1648537/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1701374/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1816173/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1952581/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1969843/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2116346/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2307642/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2401656/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2446470/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2687381/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:288218/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2906389/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2998536/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1413702/61‑1 (MQ=255) cACTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGAAGCACCGTCAGGCCa  <  1:853061/61‑1 (MQ=255)  aCTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:2096658/60‑1 (MQ=255)  aCTTTCGCTTTTATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCa  <  1:1828711/60‑1 (MQ=255)             |                                                 CACTTTCGCTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCA  >  minE/1593522‑1593582

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