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%	115.7 / NA	32	Q167R (CAG→CGG)	hlyE	hemolysin E
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/27);  total (5/27)

GCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAG  >  minE/782067‑782164                                           |                                                        gCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAAt                              >  1:1465998/1‑70 (MQ=255) gCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAAt                              >  1:3753136/1‑70 (MQ=255) gCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAAt                              >  1:2239953/1‑70 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGcttctt                                        <  1:4237177/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:2519833/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:3447673/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:3731796/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:3917956/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:4059085/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:2791560/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:2705425/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:435306/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:2054904/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:69064/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:1768687/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:1665912/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:905892/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:111224/46‑1 (MQ=255)               ttccttccTGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:1295665/46‑1 (MQ=255)                tccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCttttt                                        <  1:956004/45‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:3496755/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:3240445/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:3200195/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:3107039/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:3982205/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:2775373/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:623868/69‑1 (MQ=255)                 ccttcctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGc               <  1:1911415/69‑1 (MQ=255)                     cctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCt          <  1:2863818/70‑1 (MQ=255)                     cctGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCt          <  1:2540150/70‑1 (MQ=255)                            ttATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCcag  >  1:1578938/1‑71 (MQ=255)                            ttATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCcag  >  1:1349800/1‑71 (MQ=255)                                           |                                                        GCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACTGGCTATCTAACGCCAG  >  minE/782067‑782164

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