Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,386	C→T	59.4%	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	59.4%	14.5 / 33.1	32	S288F (TCC→TTC)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base T (9/10);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGA  >  minE/327317‑327452                                                                      |                                                                   tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGATGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTcc                                                                   >  1:324227/1‑71 (MQ=255)         ggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCt                                                                 <  1:3391055/65‑1 (MQ=255)         ggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGc                                                           >  1:1525595/1‑71 (MQ=255)            tttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcg                                                         >  1:2089103/1‑70 (MQ=255)              tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcggcg                                                      >  1:1359148/1‑71 (MQ=255)              tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcgg                                                        <  1:3932999/69‑1 (MQ=255)                    ggATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaa                                                <  1:1229002/71‑1 (MQ=255)                      aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaaa                                              <  1:2799900/71‑1 (MQ=255)                      aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaaaaaa                                              <  1:950066/71‑1 (MQ=255)                      aTTGTCTGGAGACGCTGGAAGAGATT‑CCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAat                                             <  1:2692643/71‑1 (MQ=255)                         gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAata                                            <  1:53508/70‑1 (MQ=255)                            tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                         <  1:595411/70‑1 (MQ=255)                              gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTtcc                                                                  <  1:3595214/43‑1 (MQ=255)                              gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAata                                      >  1:3509793/1‑71 (MQ=255)                               agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                     <  1:2615197/71‑1 (MQ=255)                                   gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTc                                 <  1:1349961/71‑1 (MQ=255)                                   gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTc                                 <  1:1629149/71‑1 (MQ=255)                                     tGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCg                               >  1:3365178/1‑71 (MQ=255)                                       gAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc                             >  1:2109379/1‑71 (MQ=255)                                        aaGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc                             >  1:3360174/1‑70 (MQ=255)                                         agagATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAat                                       >  1:1071327/1‑59 (MQ=255)                                           agaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCg                               <  1:962987/65‑1 (MQ=255)                                            gaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTa                        <  1:430985/71‑1 (MQ=255)                                             aTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGaa                                        <  1:2747939/54‑1 (MQ=255)                                               tGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatata                                    >  1:314510/1‑56 (MQ=255)                                                 ccGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                    >  1:1268177/1‑70 (MQ=255)                                                     gCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGc               >  1:3176086/1‑71 (MQ=255)                                                      cAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc              <  1:999103/71‑1 (MQ=255)                                                       aaaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc              <  1:1596167/70‑1 (MQ=255)                                                        aaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCg             >  1:160135/1‑70 (MQ=255)                                                          aCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                    >  1:478958/1‑61 (MQ=255)                                                            cGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACa        <  1:132399/71‑1 (MQ=255)                                                              tGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACata      <  1:3430348/71‑1 (MQ=255)                                                                aGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGaa          >  1:519257/1‑65 (MQ=255)                                                                   tctcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGa  >  1:3843990/1‑70 (MQ=255)                                                                      |                                                                   TGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGA  >  minE/327317‑327452

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