Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,329 T→C 5.0% Q642R (CAG→CGG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,3290TC 5.0% 233.4 / 3.0 79Q642R (CAG→CGG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (46/29);  new base C (2/2);  total (48/31)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTG  >  minE/608264‑608397
                                                                 |                                                                    
tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                                   <  1:433874/69‑1 (MQ=255)
 aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                                >  1:2276328/1‑71 (MQ=255)
 aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                                >  1:1997276/1‑71 (MQ=255)
 aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                                >  1:3711575/1‑71 (MQ=255)
     tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                                      <  1:931244/61‑1 (MQ=255)
     tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                                      <  1:57827/61‑1 (MQ=255)
     tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                                    <  1:2276715/63‑1 (MQ=255)
     tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                                            >  1:2824860/1‑71 (MQ=255)
       tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                                          >  1:86378/1‑71 (MQ=255)
         tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         >  1:2274216/1‑70 (MQ=255)
         tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         >  1:162429/1‑70 (MQ=255)
         tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         >  1:2871084/1‑70 (MQ=255)
         tAGCGTAAGCAAACCGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         >  1:3170732/1‑70 (MQ=255)
             gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                                    <  1:3083121/71‑1 (MQ=255)
             gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                                    >  1:3206440/1‑71 (MQ=255)
              tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                                    <  1:1532937/54‑1 (MQ=255)
                aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                                 >  1:1485182/1‑71 (MQ=255)
                aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                                 >  1:3070467/1‑71 (MQ=255)
                aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                                 >  1:1218867/1‑71 (MQ=255)
                    aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                                <  1:51290/68‑1 (MQ=255)
                     aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                                                  >  1:1588358/1‑49 (MQ=255)
                      cGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAc                                           >  1:2351881/1‑71 (MQ=255)
                       ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                                          <  1:1757453/71‑1 (MQ=255)
                         cAGTTTCGCCACCGAGGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCg                                        >  1:244307/1‑71 (MQ=255)
                           gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCggc                                      >  1:442481/1‑70 (MQ=255)
                            tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                                              <  1:3418642/62‑1 (MQ=255)
                            tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                                              <  1:3423367/62‑1 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                                   >  1:2472862/1‑55 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                    >  1:2686875/1‑70 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                    >  1:2777947/1‑70 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                    >  1:249141/1‑70 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                    >  1:1308062/1‑70 (MQ=255)
                              tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                    >  1:816383/1‑70 (MQ=255)
                               aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                                <  1:3357585/57‑1 (MQ=255)
                                gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa                                  <  1:2898274/70‑1 (MQ=255)
                                 cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCGCGACCGTGCATCgg                                                   <  1:2989599/52‑1 (MQ=255)
                                  aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                                   >  1:1076369/1‑51 (MQ=255)
                                  aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg                               >  1:3648859/1‑71 (MQ=255)
                                  aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg                               <  1:1243177/71‑1 (MQ=255)
                                    cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACTAcc                                                            >  1:3876436/1‑40 (MQ=255)
                                      gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                             <  1:18356/69‑1 (MQ=255)
                                       aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg                          >  1:3561611/1‑71 (MQ=255)
                                         gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                                 <  1:2751199/46‑1 (MQ=255)
                                         gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGcc                        <  1:3416744/71‑1 (MQ=255)
                                          tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         <  1:2555676/37‑1 (MQ=255)
                                          tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                         <  1:735118/37‑1 (MQ=255)
                                          tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                                              <  1:2509218/48‑1 (MQ=255)
                                           cAGAGAGGCTACTGCACCTTTCTGGTCACGAACGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                      >  1:867294/1‑71 (MQ=255)
                                            agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:3127725/1‑59 (MQ=255)
                                            agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                      <  1:2643037/70‑1 (MQ=255)
                                            agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                      <  1:3095710/70‑1 (MQ=255)
                                            agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc                     <  1:1417372/71‑1 (MQ=255)
                                               cagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg                  <  1:2164731/70‑1 (MQ=255)
                                                agGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg                  <  1:1962260/70‑1 (MQ=255)
                                                  gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                                              <  1:813800/40‑1 (MQ=255)
                                                      cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg           >  1:2325586/1‑71 (MQ=255)
                                                      cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg           <  1:360988/71‑1 (MQ=255)
                                                       tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt          >  1:1337665/1‑71 (MQ=255)
                                                       tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt          >  1:4043511/1‑71 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:817330/1‑47 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:827428/1‑47 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:40307/1‑47 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                      <  1:2366974/58‑1 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc         >  1:2613453/1‑71 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc         >  1:2520849/1‑71 (MQ=255)
                                                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc         >  1:2494558/1‑71 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:1682477/1‑46 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:1488334/1‑46 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                 >  1:247056/1‑46 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt        <  1:3660757/71‑1 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt        >  1:960620/1‑71 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt        <  1:350474/71‑1 (MQ=255)
                                                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt        >  1:3574887/1‑71 (MQ=255)
                                                          aCCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                                          >  1:2477791/1‑36 (MQ=255)
                                                           ccTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgg      <  1:2790816/71‑1 (MQ=255)
                                                            cTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCg                                   <  1:166460/41‑1 (MQ=255)
                                                              ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCAcgacga              >  1:105284/1‑60 (MQ=255)
                                                              ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt   >  1:4116194/1‑71 (MQ=255)
                                                               tCTGGTGACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCa                    >  1:5650/1‑53 (MQ=255)
                                                               tCTGGTGACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCa                    >  1:1188/1‑53 (MQ=255)
                                                                cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGTCGCGCCAGCACGACGACCGTCTGGGgtg  <  1:1871833/70‑1 (MQ=255)
                                                                 |                                                                    
TACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTG  >  minE/608264‑608397

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