Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,134,852 G→T 48.7% L264L (CTG→CTT)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,134,8520GT48.7% 9.5 / 82.9 76L264L (CTG→CTT) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base G (14/25);  new base T (22/15);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGC  >  minE/1134786‑1134917
                                                                  |                                                                 
gCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGaa                                                                 >  1:826538/1‑69 (MQ=255)
               aaGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                  >  1:2809592/1‑69 (MQ=255)
               aaGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                  >  1:941121/1‑69 (MQ=255)
                  ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGCCACTGAAACTGGCAGAAAATGCctc                                               >  1:294127/1‑69 (MQ=255)
                  ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCctc                                               >  1:2057717/1‑69 (MQ=255)
                  ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCctc                                               >  1:574776/1‑69 (MQ=255)
                  ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCctc                                               >  1:377479/1‑69 (MQ=255)
                       gTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          <  1:3395680/69‑1 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATTGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:3332530/1‑66 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTg                                                             >  1:186239/1‑47 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:1019723/1‑66 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:2377063/1‑66 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:2933121/1‑66 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:2298999/1‑66 (MQ=255)
                          ccGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                          >  1:446582/1‑66 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:2724301/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:3270131/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:2454254/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1973284/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:2826113/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1891041/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:378441/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1528459/69‑1 (MQ=255)
                              tGGCGCTGGTGACAGCATGGTCGGCGCCATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:401268/69‑1 (MQ=255)
                               ggCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1899550/68‑1 (MQ=255)
                               ggCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1549053/68‑1 (MQ=255)
                                            gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCg                     <  1:1732995/69‑1 (MQ=255)
                                            gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCg                     <  1:1676882/69‑1 (MQ=255)
                                            gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCg                     <  1:3374952/69‑1 (MQ=255)
                                            gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCg                     <  1:1190402/69‑1 (MQ=255)
                                            gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCg                     <  1:2059465/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:864940/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:965343/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:3205669/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:800241/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:1932710/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:1503592/69‑1 (MQ=255)
                                             cATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGt                    <  1:1664327/69‑1 (MQ=255)
                                                   cGGCGCTATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:673467/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATTCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:946416/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:2980194/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1119784/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:3157269/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:57079/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:2537665/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1228901/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:2116402/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:2141953/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1517651/69‑1 (MQ=255)
                                                   cGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1424284/69‑1 (MQ=255)
                                                   cGGCGCGATGACACCTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1474269/69‑1 (MQ=255)
                                                   cGGCGCGATGAAACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:1965446/69‑1 (MQ=255)
                                                    ggCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGc              <  1:415268/68‑1 (MQ=255)
                                                     gcgcGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGAt                                   <  1:1695661/46‑1 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGagcagc  >  1:3051312/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1441551/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1469823/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:494193/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:457460/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1285642/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1781491/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:3427802/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1791778/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1276333/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:1207678/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2128571/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2506630/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2746254/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:277105/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2435953/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2280585/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2278955/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:3178560/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:3356277/1‑69 (MQ=255)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTTGAAAAGATGGTCCGt                            >  1:3168690/1‑43 (MQ=37)
                                                               acTTAAACTGGCAGAAAATGCCTCTCTT‑AAGAGATGGTCCGTTTTGGCGTAGCTGCGGGGAGTGCAGc  >  1:2452090/1‑68 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGGTGGTTCGTTTTGGCGTAGCTGCgggg          >  1:718935/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:2983361/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:3088137/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:233560/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:491337/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:2226233/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:531145/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:13223/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:2607699/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:2606490/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:120885/1‑58 (MQ=255)
                                                                  gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgggg          >  1:2533518/1‑58 (MQ=255)
                                                                  |                                                                 
GCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGC  >  minE/1134786‑1134917

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