Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,720	G→A	12.5%	I178I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,720	0	G	A	12.5%	191.1 / 15.5	72	I178I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (3/6);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTC  >  minE/609650‑609788                                                                       |                                                                     tGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                                      <  1:2263731/71‑1 (MQ=255)  gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                                     >  1:2835699/1‑71 (MQ=255)    aCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                                     <  1:1435398/69‑1 (MQ=255)        tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCTCAACGCGACGGTAGTCACCGATGATacgg                                                                <  1:3307852/70‑1 (MQ=255)        tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacgg                                                                <  1:2629354/70‑1 (MQ=255)        tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                                 <  1:404804/69‑1 (MQ=255)           cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                            >  1:2075395/1‑71 (MQ=255)           cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                            >  1:3306986/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)                   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:1682036/50‑1 (MQ=255)                           aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                                 <  1:618009/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:2225540/1‑63 (MQ=255)                                gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                               >  1:3429532/1‑63 (MQ=255)                                gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                        <  1:1835672/70‑1 (MQ=255)                                gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                       >  1:2999190/1‑71 (MQ=255)                                gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                       >  1:1429495/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:651234/71‑1 (MQ=255)                                   gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                    >  1:772263/1‑71 (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)                                                            tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg           <  1:3363619/71‑1 (MQ=255)                                                              gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg           >  1:1254477/1‑69 (MQ=255)                                                               tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                     <  1:2711657/58‑1 (MQ=255)                                                               tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                  <  1:3535560/61‑1 (MQ=255)                                                               tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg        <  1:2118094/71‑1 (MQ=255)                                                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCa                                 >  1:2631320/1‑44 (MQ=255)                                                                 aCCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAgc                               <  1:3941157/46‑2 (MQ=255)                                                                  ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt     <  1:1036060/71‑1 (MQ=255)                                                                  ccGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGt                                   <  1:2146723/41‑1 (MQ=255)                                                                   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg        >  1:2603756/1‑67 (MQ=255)                                                                   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg        >  1:762822/1‑67 (MQ=255)                                                                   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg    >  1:1591625/1‑71 (MQ=255)                                                                   cGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAtt                     <  1:1316290/54‑1 (MQ=255)                                                                    gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt   >  1:3618208/1‑71 (MQ=255)                                                                    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt   >  1:2490979/1‑71 (MQ=255)                                                                    gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat   <  1:3795941/71‑3 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                  >  1:2591318/1‑39 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                  >  1:1744646/1‑39 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc          >  1:563614/1‑63 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc          >  1:1713932/1‑63 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc  >  1:8738/1‑71 (MQ=255)                                                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag   >  1:1777735/1‑68 (MQ=255)                                                                       |                                                                     TGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTC  >  minE/609650‑609788

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