Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,355 A→G 5.9% G633G (GGT→GGC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,3550AG 5.9% 188.2 / 3.5 68G633G (GGT→GGC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (2/2);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAC  >  minE/608286‑608423
                                                                     |                                                                    
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)
        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                                              >  1:1308062/1‑70 (MQ=255)
        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                                              >  1:816383/1‑70 (MQ=255)
        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                                              >  1:2777947/1‑70 (MQ=255)
        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                                              >  1:2686875/1‑70 (MQ=255)
        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                                                              >  1:249141/1‑70 (MQ=255)
          gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa                                                            <  1:2898274/70‑1 (MQ=255)
            aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg                                                         <  1:1243177/71‑1 (MQ=255)
            aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg                                                         >  1:3648859/1‑71 (MQ=255)
                gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                                                       <  1:18356/69‑1 (MQ=255)
                 aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg                                                    >  1:3561611/1‑71 (MQ=255)
                   gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGcc                                                  <  1:3416744/71‑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)
                                cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg                                     <  1:360988/71‑1 (MQ=255)
                                cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg                                     >  1:2325586/1‑71 (MQ=255)
                                 tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt                                    >  1:4043511/1‑71 (MQ=255)
                                 tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt                                    >  1:1337665/1‑71 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:827428/1‑47 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:40307/1‑47 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:817330/1‑47 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                                                <  1:2366974/58‑1 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc                                   >  1:2494558/1‑71 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc                                   >  1:2520849/1‑71 (MQ=255)
                                  gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc                                   >  1:2613453/1‑71 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:247056/1‑46 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:1682477/1‑46 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                           >  1:1488334/1‑46 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                  >  1:960620/1‑71 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                  <  1:350474/71‑1 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                  >  1:3574887/1‑71 (MQ=255)
                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                  <  1:3660757/71‑1 (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:1188/1‑53 (MQ=255)
                                         tCTGGTGACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCa                                              >  1:5650/1‑53 (MQ=255)
                                          cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGTCGCGCCAGCACGACGACCGTCTGGGgtg                            <  1:1871833/70‑1 (MQ=255)
                                            ggTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                             <  1:783836/67‑1 (MQ=255)
                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg                                                    <  1:1515063/43‑1 (MQ=255)
                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTa                         >  1:1018001/1‑70 (MQ=255)
                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                        <  1:3367717/71‑1 (MQ=255)
                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                        <  1:1267105/71‑1 (MQ=255)
                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                        <  1:4042575/71‑1 (MQ=255)
                                                aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTc                                   >  1:1406627/1‑57 (MQ=255)
                                                aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCg                     >  1:1236691/1‑71 (MQ=255)
                                                    ccGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGtttt                 <  1:3076819/71‑1 (MQ=255)
                                                     cGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGtttt                 <  1:962110/70‑1 (MQ=255)
                                                     cGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGACAGCACGACGACCGTCTGGGGTGTTACCCGTTTTc                >  1:1836376/1‑71 (MQ=255)
                                                       tGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtg              <  1:4026670/71‑2 (MQ=255)
                                                          aTCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                   <  1:858995/63‑1 (MQ=255)
                                                          aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAc            <  1:2298779/70‑1 (MQ=255)
                                                            cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc           >  1:893114/1‑69 (MQ=255)
                                                            cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa          >  1:3221619/1‑70 (MQ=255)
                                                             gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa        <  1:178805/71‑1 (MQ=255)
                                                                 tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac     >  1:858568/1‑70 (MQ=255)
                                                                 tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac     >  1:285234/1‑70 (MQ=255)
                                                                   gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc  >  1:1806268/1‑71 (MQ=255)
                                                                     |                                                                    
CGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAC  >  minE/608286‑608423

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