Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	232,635	T→C	55.6%	N43D (AAC→GAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	232,635	0	T	C	55.6%	22.8 / 75.3	63	N43D (AAC→GAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base C (19/16);  total (35/28)
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

TCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAA  >  minE/232565‑232701                                                                       |                                                                   tCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGtt                                                                    >  1:1212220/1‑71 (MQ=255)  cAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGtt                                                                    <  1:327285/70‑1 (MQ=255)     tttCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCa                                                                <  1:2987566/71‑1 (MQ=255)        cgcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAAt                                                              <  1:2812139/70‑1 (MQ=255)           gTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATttct                                                          >  1:2628671/1‑71 (MQ=255)            tGCCAGATGTTTCTCTGGAATGCGCATCACTCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttctt                                                         <  1:1651130/71‑1 (MQ=255)               cAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttctt                                                      >  1:329818/1‑71 (MQ=255)                aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATttcttctt                                                      >  1:439843/1‑70 (MQ=255)                   tGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaa                                                   >  1:1643437/1‑70 (MQ=255)                      ttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                 >  1:3499514/1‑69 (MQ=255)                      ttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaaca                                                 >  1:8063/1‑69 (MQ=255)                       tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                 >  1:2361702/1‑68 (MQ=255)                           tGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGa                                             <  1:1682554/68‑1 (MQ=255)                           tGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaaca                                                 >  1:358527/1‑64 (MQ=255)                           tGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGAt                                            <  1:466541/69‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATc                                           <  1:461391/68‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAAt                                                              <  1:1775801/49‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAAt                                                              <  1:3294026/49‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAAt                                                              <  1:3474733/49‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATc                                           <  1:2587540/68‑1 (MQ=255)                             gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCgg                                         >  1:1109367/1‑70 (MQ=255)                              aaTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCa                                       <  1:3596361/71‑1 (MQ=255)                              aaTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCa                                       <  1:2003279/71‑1 (MQ=255)                                tGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCaa                                      >  1:3205452/1‑70 (MQ=255)                                tGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCaa                                      >  1:3563028/1‑70 (MQ=255)                                tGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaac                                                  >  1:2264943/1‑58 (MQ=255)                                  cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGtt                                                                    >  1:991623/1‑38 (MQ=255)                                  cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGtt                                                                    >  1:3512305/1‑38 (MQ=255)                                  cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGAt                                            >  1:1241731/1‑62 (MQ=255)                                  cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGAt                                            >  1:1735105/1‑62 (MQ=255)                                   gcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                  <  1:3770999/71‑1 (MQ=255)                                   gcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                  <  1:2035025/71‑1 (MQ=255)                                     aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAAAAAAGATCGGCCACCCCAg                                >  1:1029329/1‑71 (MQ=255)                                       cACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt                              >  1:1208061/1‑71 (MQ=255)                                       cACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt                              >  1:1566417/1‑71 (MQ=255)                                       cACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaacaa                                                >  1:2260529/1‑53 (MQ=255)                                        aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCGATTTCTTCTTCAACAAAGATCGGCAACACCAGGTc                             >  1:273143/1‑71 (MQ=255)                                         cGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCg                            <  1:1551413/71‑1 (MQ=255)                                         cGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTGTTCAACAAAGATCGGCAaca                                    >  1:1096269/1‑63 (MQ=255)                                          gCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCg                                          >  1:1043892/1‑56 (MQ=255)                                           ccTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                 <  1:923597/48‑1 (MQ=255)                                           ccTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGtt                          <  1:2111603/71‑1 (MQ=255)                                               gCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                      >  1:1219260/1‑71 (MQ=255)                                               gCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                      >  1:781103/1‑71 (MQ=255)                                               gCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                      >  1:776579/1‑71 (MQ=255)                                                 aTGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt                    <  1:1499125/71‑1 (MQ=255)                                                  tGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa                   <  1:1892248/71‑1 (MQ=255)                                                  tGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaa                                                   >  1:1272841/1‑39 (MQ=255)                                                      ttCAAGGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg               >  1:850670/1‑71 (MQ=255)                                                       tCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtgt              >  1:1788619/1‑71 (MQ=255)                                                        cAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg               <  1:1461459/69‑1 (MQ=255)                                                        cAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgtt             >  1:2287923/1‑71 (MQ=255)                                                         aaCGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgttg            >  1:781488/1‑71 (MQ=255)                                                         aaCGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgtt             <  1:3778594/70‑1 (MQ=255)                                                           cGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa                   >  1:2117246/1‑62 (MQ=255)                                                            ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGAt                                            >  1:2262845/1‑36 (MQ=255)                                                            ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGt                           <  1:659177/53‑1 (MQ=255)                                                            ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                      >  1:1351376/1‑58 (MQ=255)                                                            ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct         <  1:2882777/71‑1 (MQ=255)                                                             gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtctc        <  1:299191/71‑1 (MQ=255)                                                               tttGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAgg                               <  1:1435913/46‑1 (MQ=255)                                                                 tGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct         >  1:2409214/1‑66 (MQ=255)                                                                   tAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtctc        >  1:1352717/1‑65 (MQ=255)                                                                   tAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtctc        >  1:1058798/1‑65 (MQ=255)                                                                   tAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAaca                                    <  1:1109756/37‑1 (MQ=255)                                                                   tAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCaaa  <  1:1846564/71‑1 (MQ=255)                                                                   tAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCaaa  <  1:2158986/71‑1 (MQ=255)                                                                       |                                                                   TCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAA  >  minE/232565‑232701

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