Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	782,109	T→C	Q167R (CAG→CGG)	hlyE ←	hemolysin E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	782,109	0	T	C	100.0%	121.6 / NA	34	Q167R (CAG→CGG)	hlyE	hemolysin E
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (1/33);  total (1/33)

ACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACT  >  minE/782055‑782148                                                       |                                        acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:134197/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:855061/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:706162/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3203627/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3162799/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:305061/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3049445/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2878878/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2793347/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2751445/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2661640/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2630271/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2610430/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2608326/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2551816/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2381529/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1156249/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1534309/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1564377/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1608109/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1625095/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1647366/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1671332/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1920146/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1921523/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1942587/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2147074/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2295328/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2299392/71‑1 (MQ=255) acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTCTGACCGGAAATAGCTGCtttt                         <  1:3140874/71‑1 (MQ=255)  cacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:580782/70‑1 (MQ=255)           cGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2686987/61‑1 (MQ=255)               aCCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2949256/57‑1 (MQ=255)                        tGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACt  >  1:886534/1‑71 (MQ=255)                                                       |                                        ACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACT  >  minE/782055‑782148

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