Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,274 A→G 8.0% A660A (GCT→GCC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2740AG 8.0% 208.0 / 7.8 75A660A (GCT→GCC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base G (2/4);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

TCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTG  >  minE/608204‑608342
                                                                      |                                                                    
tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa                                                                      <  1:3187426/71‑1 (MQ=255)
tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa                                                                      <  1:3325699/71‑1 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGc                                                                    >  1:2391648/1‑71 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGc                                                                    >  1:829027/1‑71 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGc                                                                    >  1:2641297/1‑71 (MQ=255)
    aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   >  1:2193365/1‑70 (MQ=255)
    aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   >  1:1431312/1‑70 (MQ=255)
     cGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:1569136/70‑1 (MQ=255)
        aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                               >  1:3709296/1‑70 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:3872484/66‑1 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:2255876/71‑1 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:1712302/71‑1 (MQ=255)
         cTTCGGCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:3877025/66‑1 (MQ=255)
            cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                          >  1:1025265/1‑71 (MQ=255)
            cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                          >  1:3407785/1‑71 (MQ=255)
             gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                          <  1:3206845/70‑1 (MQ=255)
                gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                        >  1:2550577/1‑69 (MQ=255)
                 tcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACAATCTTTAGCGTAAGCAAACGGCAg                                                     <  1:1804639/71‑1 (MQ=255)
                  cTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                                    >  1:1149312/1‑71 (MQ=255)
                     tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                  <  1:3813684/70‑1 (MQ=255)
                     tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                  <  1:821148/70‑1 (MQ=255)
                       cccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                 <  1:1176120/69‑1 (MQ=255)
                        ccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:3218834/51‑1 (MQ=255)
                         cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:1561477/50‑1 (MQ=255)
                         cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:2795541/50‑1 (MQ=255)
                         cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                 >  1:1048478/1‑67 (MQ=255)
                          tgTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:2371406/48‑1 (MQ=255)
                          aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                            >  1:7101/1‑71 (MQ=255)
                           gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                           <  1:3339879/71‑1 (MQ=255)
                            tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt                                                                  <  1:3541302/47‑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)
                                  cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   <  1:2998825/40‑1 (MQ=255)
                                  cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                                         >  1:101027/1‑50 (MQ=255)
                                  cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                      >  1:3140270/1‑69 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                             >  1:1670850/1‑61 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:3333443/71‑1 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:2044591/71‑1 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:694873/71‑1 (MQ=255)
                                    gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:3827985/70‑1 (MQ=255)
                                     aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                    >  1:1412528/1‑68 (MQ=255)
                                     aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                    >  1:3352267/1‑68 (MQ=255)
                                     aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                    >  1:1982606/1‑68 (MQ=255)
                                     aaCGATAGAGAAGGTGGAGGAGATACCATCTTTAGCGTa                                                                 <  1:251144/39‑1 (MQ=255)
                                      aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                  >  1:2992609/1‑67 (MQ=255)
                                      aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                <  1:324297/71‑1 (MQ=255)
                                      aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                  >  1:3283272/1‑69 (MQ=255)
                                      aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                  >  1:3302705/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)
                                                                      |                                                                    
TCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTG  >  minE/608204‑608342

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