Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,708	G→A	5.5%	Y182Y (TAC→TAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,708	0	G	A	5.5%	226.1 / 2.3	72	Y182Y (TAC→TAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (39/29);  new base A (2/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAA  >  minE/609638‑609777                                                                       |                                                                      ccTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                       <  1:62663/71‑1 (MQ=255)  cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                       <  1:3103040/70‑1 (MQ=255)  cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGt                                                                      <  1:2508675/71‑1 (MQ=255)  cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGt                                                                      <  1:2508047/71‑1 (MQ=255)     cAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                   >  1:2879555/1‑71 (MQ=255)     cAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                   >  1:2624355/1‑71 (MQ=255)          aaGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCAcc                                                               >  1:1064537/1‑70 (MQ=255)             tGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                           <  1:2263731/71‑1 (MQ=255)              gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTc                                                                  <  1:3950375/63‑1 (MQ=255)              gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                          >  1:2835699/1‑71 (MQ=255)                aCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                          <  1:1435398/69‑1 (MQ=255)                 cTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCa                                                                 >  1:2789970/1‑61 (MQ=255)                  tgtgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                       <  1:3229497/54‑1 (MQ=255)                    tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCTCAACGCGACGGTAGTCACCGATGATacgg                                                     <  1:3307852/70‑1 (MQ=255)                    tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacgg                                                     <  1:2629354/70‑1 (MQ=255)                    tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                      <  1:404804/69‑1 (MQ=255)                       cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAg                                                                    >  1:2959841/1‑52 (MQ=255)                       cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                   >  1:3442974/1‑53 (MQ=255)                       cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                 >  1:3306986/1‑71 (MQ=255)                       cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                 >  1:2075395/1‑71 (MQ=255)                       cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                 >  1:3202839/1‑71 (MQ=255)                        aGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacg                                                >  1:1162447/1‑71 (MQ=255)                        aGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacg                                                >  1:1433351/1‑71 (MQ=255)                         gTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacgg                                               >  1:3906498/1‑71 (MQ=255)                         gTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacgg                                               >  1:3550676/1‑71 (MQ=255)                             gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggcc                                             <  1:1633531/69‑1 (MQ=255)                               cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCa                                                                 <  1:3920702/47‑1 (MQ=255)                               cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACAGCCACGGCCatat                                         >  1:3121457/1‑71 (MQ=255)                                tttCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                           >  1:1992020/1‑52 (MQ=255)                                  tcatcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                      >  1:2098277/1‑71 (MQ=255)                                   catcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                     >  1:3428461/1‑71 (MQ=255)                                     tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                   <  1:2305776/71‑1 (MQ=255)                                       aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                      <  1:618009/50‑1 (MQ=255)                                       aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                      <  1:1682036/50‑1 (MQ=255)                                       aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                 >  1:3538684/1‑71 (MQ=255)                                       aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                 >  1:1033102/1‑71 (MQ=255)                                        ggTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                                                   >  1:2380777/1‑52 (MQ=255)                                        ggTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                 >  1:1976427/1‑70 (MQ=255)                                         gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCa                               <  1:2268395/71‑1 (MQ=255)                                           aGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                             >  1:878512/1‑71 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCata                                          >  1:3657247/1‑57 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                    >  1:3429532/1‑63 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                    >  1:2225540/1‑63 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                             <  1:1835672/70‑1 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                            >  1:1429495/1‑71 (MQ=255)                                            gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                            >  1:2999190/1‑71 (MQ=255)                                             tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATac                                                       <  1:1141364/43‑1 (MQ=255)                                             tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                            <  1:914329/54‑1 (MQ=255)                                             tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCata                                          >  1:961492/1‑56 (MQ=255)                                             ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                  >  1:208337/3‑64 (MQ=255)                                             ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                  >  1:941580/3‑64 (MQ=255)                                               gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                         >  1:772263/1‑71 (MQ=255)                                               gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                         <  1:651234/71‑1 (MQ=255)                                               gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                         >  1:1975028/1‑71 (MQ=255)                                               gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCg                          <  1:3947610/70‑1 (MQ=255)                                                 tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggcc                                             <  1:2008917/49‑1 (MQ=255)                                                 tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                            >  1:1150230/1‑66 (MQ=255)                                                 aaCCGTACAGCGCACCGCGACGGTAGTCACCGATGATacggccacggcc                                             <  1:3550818/48‑1 (MQ=255)                                                  aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                         >  1:3559219/1‑68 (MQ=255)                                                  aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                        >  1:2029271/1‑69 (MQ=255)                                                   ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGAATCTGGCAGACCGGTCa                      <  1:2138738/70‑1 (MQ=255)                                                    cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGa                    <  1:614699/71‑1 (MQ=255)                                                     gTACAGCGCAACGCGACGGTAGTCACCGATGATacggcca                                                  >  1:3048815/1‑40 (MQ=255)                                                     gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                   <  1:2199452/71‑1 (MQ=255)                                                        cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                   <  1:1353960/68‑1 (MQ=255)                                                         aGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGcc               <  1:3794751/71‑1 (MQ=255)                                                            gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa            <  1:2344681/68‑1 (MQ=255)                                                            gcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                   <  1:1224188/48‑1 (MQ=255)                                                               aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACAcc               <  1:3513508/65‑1 (MQ=255)                                                               aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAttt         <  1:3367157/71‑1 (MQ=255)                                                                cgcgACGGTAGTCACCGATGATACGGCCACGGCCata                                          >  1:1296742/1‑37 (MQ=255)                                                                cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa        >  1:1437964/1‑71 (MQ=255)                                                                cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa        <  1:3746778/71‑1 (MQ=255)                                                                    aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg     >  1:3381941/1‑70 (MQ=255)                                                                     cGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc    >  1:2347202/1‑70 (MQ=255)                                                                      ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa  >  1:277157/1‑71 (MQ=255)                                                                       |                                                                      CCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAA  >  minE/609638‑609777

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