Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,760 G→A 27.9% I498I (ATC→ATT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,7600GA27.9% 174.2 / 77.7 111I498I (ATC→ATTpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (44/36);  new base A (16/15);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 125 total aligned reads displayed.

TGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608690‑608829
                                                                      |                                                                     
tGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                                                       <  1:3700027/71‑1 (MQ=255)
  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCTTGCATGTAGTGga                                                                      <  1:180471/70‑1 (MQ=255)
  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:3045120/70‑1 (MQ=255)
  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:115108/70‑1 (MQ=255)
   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                                                    <  1:1051507/71‑1 (MQ=255)
   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                                                     >  1:1712621/1‑70 (MQ=255)
   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                                                     >  1:2498510/1‑70 (MQ=255)
    gATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   <  1:1510895/71‑1 (MQ=255)
      tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                                                 <  1:2111386/71‑1 (MQ=255)
       aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:2848019/1‑71 (MQ=255)
        acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                                                               <  1:2451649/68‑1 (MQ=255)
        aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                <  1:284503/70‑1 (MQ=255)
          gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                             <  1:3765323/71‑1 (MQ=255)
          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                                                               >  1:1375221/2‑69 (MQ=255)
            cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagc                                                           >  1:1683013/1‑70 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:583284/71‑1 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:2616968/71‑1 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:3659170/71‑1 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:2448796/1‑71 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:3690028/1‑71 (MQ=255)
             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt                                                           >  1:1746528/1‑70 (MQ=255)
               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                               <  1:1271516/64‑1 (MQ=255)
               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                               <  1:923735/64‑1 (MQ=255)
               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:3745091/1‑69 (MQ=255)
               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                              >  1:2678580/1‑65 (MQ=255)
                gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                             <  1:827209/65‑1 (MQ=255)
                atCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgc                                                       <  1:1986661/70‑3 (MQ=255)
                atCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgc                                                       <  1:2285052/70‑3 (MQ=255)
                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                                                      >  1:1361555/1‑71 (MQ=255)
                  cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      >  1:3892463/1‑67 (MQ=255)
                  cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      >  1:3818077/1‑67 (MQ=255)
                   gTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                                                            <  1:3051711/63‑1 (MQ=255)
                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                     <  1:3364346/70‑1 (MQ=255)
                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:3172857/1‑71 (MQ=255)
                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:524081/1‑71 (MQ=255)
                    tGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                     <  1:1989274/69‑1 (MQ=255)
                     gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGAt                                                   <  1:1474046/70‑1 (MQ=255)
                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgat                                                                  >  1:3949500/1‑54 (MQ=255)
                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  >  1:1270476/1‑70 (MQ=255)
                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                                                 <  1:1276183/71‑1 (MQ=255)
                         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                              >  1:1974860/1‑71 (MQ=255)
                           ccATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                            <  1:2932575/71‑1 (MQ=255)
                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:3283633/71‑1 (MQ=255)
                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:1634670/71‑1 (MQ=255)
                             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                          <  1:2296146/55‑3 (MQ=255)
                             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac                                          >  1:1489148/1‑68 (MQ=255)
                              tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTAc                                               <  1:1963599/65‑1 (MQ=255)
                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                                                       <  1:7523/41‑1 (MQ=255)
                               cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:3129243/41‑1 (MQ=255)
                               cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag                                        >  1:258021/1‑71 (MQ=255)
                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                            >  1:925830/1‑50 (MQ=255)
                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       >  1:638018/1‑71 (MQ=255)
                                 gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   <  1:2817946/42‑1 (MQ=255)
                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:213058/1‑44 (MQ=255)
                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:3791007/1‑44 (MQ=255)
                                   gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                                                 <  1:326538/42‑1 (MQ=255)
                                    agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:3524790/1‑71 (MQ=255)
                                      gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:3910658/1‑69 (MQ=255)
                                      gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:1012830/1‑69 (MQ=255)
                                        ttCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                        <  1:1460678/46‑1 (MQ=255)
                                         tCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                         >  1:2267359/1‑60 (MQ=255)
                                          cGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGt                               <  1:972323/69‑1 (MQ=255)
                                           gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCATCGCGCTGATGTACTGCACCGCCAGCCAGTCCa                            >  1:3113024/1‑71 (MQ=255)
                                            tAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCAt                           <  1:2457256/71‑1 (MQ=255)
                                            tAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCAt                           <  1:3509691/71‑1 (MQ=255)
                                             aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCAt                           >  1:3492343/1‑70 (MQ=255)
                                              gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          <  1:321491/70‑1 (MQ=255)
                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:563680/50‑1 (MQ=255)
                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:129318/50‑1 (MQ=255)
                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          >  1:1187658/1‑69 (MQ=255)
                                                 gTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      >  1:691787/1‑71 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                            >  1:3491816/1‑48 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      <  1:3668348/70‑1 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      >  1:1843002/1‑70 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTg                     <  1:2744294/71‑1 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                      >  1:3348768/1‑70 (MQ=255)
                                                  tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                      >  1:723764/1‑70 (MQ=255)
                                                    cTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                              >  1:1219141/1‑44 (MQ=25)
                                                     ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  >  1:2524114/1‑71 (MQ=255)
                                                      tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                   >  1:2902988/1‑69 (MQ=255)
                                                       gTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                  <  1:2918680/69‑1 (MQ=255)
                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGt                                                 <  1:1544398/37‑1 (MQ=255)
                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGt                                                 <  1:1612665/37‑1 (MQ=255)
                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:2311670/43‑1 (MQ=255)
                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                <  1:2677928/70‑1 (MQ=255)
                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                       >  1:1007084/1‑63 (MQ=255)
                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcgag               >  1:2502151/1‑68 (MQ=255)
                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcga                <  1:726165/70‑3 (MQ=255)
                                                             cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   <  1:2457636/46‑1 (MQ=255)
                                                             cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt          >  1:1615469/1‑71 (MQ=255)
                                                              aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc         >  1:787339/1‑71 (MQ=255)
                                                              aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc         <  1:1118954/71‑1 (MQ=255)
                                                              aTGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                       >  1:1027596/1‑57 (MQ=255)
                                                                gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa       >  1:2264940/1‑71 (MQ=255)
                                                                gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGAGATCCATGCGCTcc        >  1:2305669/1‑70 (MQ=255)
                                                                 tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      <  1:2640792/71‑1 (MQ=255)
                                                                 tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      <  1:2956175/71‑1 (MQ=255)
                                                                   gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      >  1:3067653/1‑69 (MQ=255)
                                                                    tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc   <  1:2996861/71‑1 (MQ=255)
                                                                     ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt  >  1:1342334/1‑71 (MQ=255)
                                                                      |                                                                     
TGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608690‑608829

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: