Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,386	C→T	66.7%	S288F (TCC→TTC)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	327,386	0	C	T	66.7%	13.8 / 8.9	24	S288F (TCC→TTC)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (10/6);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG  >  minE/327337‑327442                                                  |                                                         gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaa                                        >  1:1864116/1‑68 (MQ=255) gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGa                                         >  1:1219268/1‑67 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:2795573/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:462241/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:3397521/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:3389383/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:3340371/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:3029049/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:226273/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:1569480/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:1555911/1‑68 (MQ=255)   ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                      >  1:1337251/1‑68 (MQ=255)             aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:908913/69‑1 (MQ=255)             aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:2769788/69‑1 (MQ=255)             aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:1003171/69‑1 (MQ=255)             aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:2888558/69‑1 (MQ=255)             aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:524140/69‑1 (MQ=255)             aCGCTGCAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                           <  1:1678927/69‑1 (MQ=255)                      gagaTTGCCGAGCAAAACCGTGAGGTCTCCCTCgg                                                    >  1:2858176/1‑35 (MQ=255)                      gagaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc                   >  1:80764/1‑68 (MQ=255)                        gaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcgc                 >  1:3313461/1‑68 (MQ=255)                               gAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc          <  1:2415645/68‑1 (MQ=255)                                      aCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg  <  1:3088028/69‑1 (MQ=255)                                      aCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg  <  1:238497/69‑1 (MQ=255)                                                ctcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                  <  1:1370512/43‑1 (MQ=255)                                                ctcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                  <  1:820704/43‑1 (MQ=255)                                                 tcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                  <  1:2512692/42‑1 (MQ=255)                                                  |                                                         GATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG  >  minE/327337‑327442

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