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
*minE609,7200GA12.5% 191.1 / 15.5 72I178I (ATC→ATTpflBpyruvate 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: