Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,295 A→C 11.9% A653A (GCT→GCG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2950AC11.9% 159.8 / 20.6 76A653A (GCT→GCG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base C (3/6);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

CCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCC  >  minE/608227‑608363
                                                                    |                                                                    
cccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                                      <  1:1176120/69‑1 (MQ=255)
  cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                                      >  1:1048478/1‑67 (MQ=255)
   aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                                 >  1:7101/1‑71 (MQ=255)
    gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                                <  1:3339879/71‑1 (MQ=255)
     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               <  1:269442/71‑1 (MQ=255)
     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               <  1:1722830/71‑1 (MQ=255)
     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTATAGc                                                                    >  1:2105193/1‑66 (MQ=255)
       cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               >  1:3867587/1‑69 (MQ=255)
       cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               >  1:3834838/1‑69 (MQ=255)
        gTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAg                                                            >  1:369184/1‑71 (MQ=255)
          tCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                           >  1:361587/1‑70 (MQ=255)
           cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                           >  1:3140270/1‑69 (MQ=255)
            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                                  >  1:1670850/1‑61 (MQ=255)
            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:694873/71‑1 (MQ=255)
            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:2044591/71‑1 (MQ=255)
            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:3333443/71‑1 (MQ=255)
             gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:3827985/70‑1 (MQ=255)
              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         >  1:1982606/1‑68 (MQ=255)
              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         >  1:3352267/1‑68 (MQ=255)
              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         >  1:1412528/1‑68 (MQ=255)
               aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                       >  1:2992609/1‑67 (MQ=255)
               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                     <  1:324297/71‑1 (MQ=255)
               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       >  1:3302705/1‑69 (MQ=255)
               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       >  1:3283272/1‑69 (MQ=255)
                  attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                  <  1:2742794/68‑1 (MQ=255)
                  attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                  <  1:267180/68‑1 (MQ=255)
                  aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                  <  1:2445083/71‑1 (MQ=255)
                   tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               >  1:1607299/1‑57 (MQ=255)
                   tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                  >  1:428323/1‑70 (MQ=255)
                    agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                                      >  1:2221826/1‑65 (MQ=255)
                      agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                  <  1:2212037/67‑1 (MQ=255)
                      agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACt                                              >  1:2131970/1‑71 (MQ=255)
                        aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                             >  1:523529/1‑54 (MQ=255)
                        aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                             >  1:1123783/1‑70 (MQ=255)
                            tgtACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGac                                                             <  1:2837786/50‑2 (MQ=255)
                             gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                       >  1:953217/1‑71 (MQ=255)
                              tacgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCttt                                      <  1:1753908/68‑1 (MQ=255)
                                ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                     <  1:2401660/70‑1 (MQ=255)
                                    aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                <  1:1654394/71‑1 (MQ=255)
                                      aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                           <  1:2029132/42‑1 (MQ=255)
                                      aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                >  1:1331851/1‑69 (MQ=255)
                                        cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGa                            >  1:720036/1‑71 (MQ=255)
                                         aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                  <  1:3942664/48‑1 (MQ=255)
                                         aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACATTTCTGGTCACGAc                           >  1:3260619/1‑71 (MQ=255)
                                          tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAc                           >  1:1051499/1‑70 (MQ=255)
                                            tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                        >  1:1980820/1‑71 (MQ=255)
                                              tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGTACCTTTCTGGTCACGACCGTg                       >  1:87529/1‑70 (MQ=255)
                                              tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                          <  1:2037412/67‑1 (MQ=255)
                                                 cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                   <  1:3250035/71‑1 (MQ=255)
                                                 cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTAACGACCGTGc                      <  1:913851/68‑1 (MQ=255)
                                                 cGTAAGAAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                   <  1:3164314/71‑1 (MQ=255)
                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                     <  1:3559758/51‑1 (MQ=255)
                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                 >  1:1019022/1‑71 (MQ=255)
                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAATg                                             <  1:355408/43‑1 (MQ=255)
                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                            >  1:2251034/1‑42 (MQ=255)
                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCATCTTTCTGGTCACGACCGTGCATCGGGt               <  1:3696793/71‑1 (MQ=255)
                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt               >  1:850934/1‑71 (MQ=255)
                                                     aGCAAACGGCAGTTTAGCAACGGAAGGCAGCGCGGCTACTGCACCTTTCTGGTCACGACCGTGCATCggg                <  1:3344432/70‑1 (MQ=255)
                                                      gtgaaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGtt              <  1:2517230/68‑1 (MQ=255)
                                                      gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                 <  1:1104843/52‑1 (MQ=255)
                                                      gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                 >  1:1086037/1‑68 (MQ=255)
                                                      gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt              <  1:3646045/71‑1 (MQ=255)
                                                         aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAc                                               >  1:1135848/1‑35 (MQ=255)
                                                         aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                           >  1:1794823/1‑39 (MQ=255)
                                                         aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                   <  1:1799456/63‑1 (MQ=255)
                                                          aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCa          <  1:552843/71‑1 (MQ=255)
                                                           cGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcgc         >  1:2072388/1‑69 (MQ=255)
                                                            ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc        <  1:1770690/71‑1 (MQ=255)
                                                               aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg     <  1:1867519/71‑1 (MQ=255)
                                                               aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg     <  1:539325/71‑1 (MQ=255)
                                                                gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCCTGCATCGGGTTAGCGCCCggc    <  1:3263540/71‑2 (MQ=255)
                                                                 tttAGCATCGGAAGTCTGAGAGGCTACTGCACCTTTc                                     <  1:2397951/37‑1 (MQ=255)
                                                                 tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                          <  1:2470539/48‑1 (MQ=255)
                                                                 tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc   >  1:3031956/1‑71 (MQ=255)
                                                                   tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:2954541/1‑70 (MQ=255)
                                                                   tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:3643701/1‑70 (MQ=255)
                                                                   tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:878034/1‑70 (MQ=255)
                                                                   tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc   >  1:2297856/1‑69 (MQ=255)
                                                                    aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc        >  1:323299/1‑63 (MQ=255)
                                                                    |                                                                    
CCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCC  >  minE/608227‑608363

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