Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,115 A→C 16.0% R713R (CGT→CGG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1150AC16.0% 160.4 / 31.9 81R713R (CGT→CGG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (39/29);  new base C (2/11);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCC  >  minE/608045‑608184
                                                                      |                                                                     
cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                                       >  1:2963865/1‑71 (MQ=255)
    gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                   <  1:3744660/71‑1 (MQ=255)
      cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                                                    <  1:2729242/68‑2 (MQ=9)
       agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg                                                                     <  1:263141/66‑1 (MQ=9)
       aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                >  1:396873/1‑71 (MQ=255)
       aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                <  1:385533/71‑1 (MQ=255)
        gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                               <  1:803100/68‑1 (MQ=9)
         aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                               >  1:3380607/1‑70 (MQ=255)
           acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                            >  1:1149639/1‑71 (MQ=255)
            cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                            <  1:251882/70‑1 (MQ=255)
               ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                                        <  1:1812077/71‑2 (MQ=255)
               ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                         >  1:2612834/1‑70 (MQ=255)
                     tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                                      <  1:160098/51‑1 (MQ=255)
                     tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                  <  1:928150/71‑1 (MQ=255)
                      cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                   >  1:254837/1‑69 (MQ=255)
                       aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                >  1:1010417/1‑71 (MQ=255)
                       aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                <  1:3241491/71‑1 (MQ=255)
                       aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                >  1:3155460/1‑71 (MQ=255)
                        gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                              >  1:2916952/1‑56 (MQ=255)
                        gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                 <  1:1912623/69‑1 (MQ=255)
                         cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                              >  1:40621/1‑71 (MQ=255)
                         cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                              <  1:2538351/71‑1 (MQ=255)
                         cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                              <  1:3583132/71‑1 (MQ=255)
                           taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                            >  1:1934025/3‑70 (MQ=255)
                           gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                            >  1:3690568/1‑71 (MQ=255)
                             ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                                      >  1:2967272/1‑43 (MQ=255)
                              gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                         <  1:1929460/71‑1 (MQ=255)
                              gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                         >  1:616967/1‑71 (MQ=12)
                                 atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                                          <  1:1328372/51‑1 (MQ=255)
                                 atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                      >  1:564393/1‑71 (MQ=255)
                                  tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                     <  1:2149654/71‑1 (MQ=255)
                                   ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                    >  1:3251582/1‑71 (MQ=255)
                                   ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                    >  1:534188/1‑71 (MQ=255)
                                   ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                    >  1:855116/1‑71 (MQ=255)
                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                    <  1:3278536/70‑1 (MQ=255)
                                        cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                               <  1:868636/71‑1 (MQ=255)
                                        cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                               <  1:1278827/71‑1 (MQ=255)
                                         gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                               <  1:558239/38‑1 (MQ=255)
                                         gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                              <  1:1627687/71‑1 (MQ=255)
                                           gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                             >  1:1553947/1‑70 (MQ=255)
                                            ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                <  1:237535/50‑1 (MQ=255)
                                            ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                <  1:1953432/50‑1 (MQ=255)
                                            ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                           >  1:1659397/1‑71 (MQ=255)
                                              ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                          >  1:1977977/1‑70 (MQ=255)
                                               tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                        >  1:3689004/1‑71 (MQ=255)
                                                ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc                        <  1:3419300/70‑1 (MQ=255)
                                                ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                        <  1:1213134/70‑1 (MQ=255)
                                                 cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                                               <  1:904983/46‑1 (MQ=255)
                                                 cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt                      >  1:3171065/1‑71 (MQ=255)
                                                  aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                  >  1:1568133/1‑58 (MQ=255)
                                                  aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                  >  1:1577473/1‑58 (MQ=255)
                                                  aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                     >  1:350886/1‑71 (MQ=255)
                                                   tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                    >  1:136118/1‑71 (MQ=255)
                                                    cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                >  1:1356558/1‑42 (MQ=255)
                                                     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                             >  1:2517574/1‑44 (MQ=255)
                                                     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                             >  1:3904503/1‑44 (MQ=255)
                                                     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                             >  1:1872647/1‑44 (MQ=255)
                                                       gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                         >  1:3597212/1‑46 (MQ=255)
                                                       gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                >  1:2158204/1‑71 (MQ=255)
                                                        tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                               <  1:2830291/55‑1 (MQ=255)
                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                           >  1:3094586/1‑41 (MQ=255)
                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                     <  1:37952/63‑1 (MQ=255)
                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:1775066/70‑1 (MQ=255)
                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              >  1:1551415/1‑70 (MQ=255)
                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:1361846/70‑1 (MQ=255)
                                                          gAGCAGCAATTCACGGTTCATCACGTTAACGTTCAg                                                >  1:1227772/1‑36 (MQ=255)
                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg              <  1:3562883/69‑1 (MQ=255)
                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa             <  1:1099635/70‑1 (MQ=255)
                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa             <  1:1307701/70‑1 (MQ=255)
                                                           agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                    <  1:1502184/47‑1 (MQ=255)
                                                           agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa            <  1:2845961/71‑1 (MQ=12)
                                                           agcagcATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                         >  1:2861022/1‑42 (MQ=255)
                                                            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa           <  1:2209198/71‑2 (MQ=255)
                                                             cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt          <  1:410139/71‑1 (MQ=255)
                                                               gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa        >  1:686513/1‑71 (MQ=255)
                                                               gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa        >  1:690240/1‑71 (MQ=255)
                                                                 aTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc      <  1:3187811/71‑1 (MQ=255)
                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                <  1:3859320/60‑1 (MQ=255)
                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca     >  1:653763/1‑71 (MQ=255)
                                                                   ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat    >  1:1077683/1‑71 (MQ=255)
                                                                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                          <  1:2052940/48‑1 (MQ=255)
                                                                     cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                >  1:923323/1‑57 (MQ=255)
                                                                      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc  <  1:597418/70‑1 (MQ=255)
                                                                      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc  <  1:1859894/70‑1 (MQ=255)
                                                                      |                                                                     
CGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCC  >  minE/608045‑608184

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