Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,174 A→G 38.7% A52A (GCA→GCG secA → preprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE93,1740AG38.7% 10.2 / 17.3 31A52A (GCA→GCGsecApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base A (11/8);  new base G (5/7);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

GAACCGGAGATGGAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCAA  >  minE/93112‑93241
                                                              |                                                                   
gAACCGGAGATGGAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCt                                                                >  1:3297490/1‑68 (MQ=255)
gAACCGGAGATGGAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCt                                                                >  1:517903/1‑68 (MQ=255)
gAACCGGAGATGGAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCt                                                                >  1:2084275/1‑68 (MQ=255)
           ggAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:3300890/69‑1 (MQ=255)
           ggAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:1714502/69‑1 (MQ=255)
           ggAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:1963669/69‑1 (MQ=255)
           ggAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:2595457/69‑1 (MQ=255)
           ggAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:3427086/69‑1 (MQ=255)
            gAAAAACTCTCCGACGAAGAACTGAAATGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:3040216/68‑1 (MQ=255)
            gAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGaa                                                   >  1:822229/1‑69 (MQ=255)
            gAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGaa                                                   >  1:2152640/1‑69 (MQ=255)
            gAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:2315153/68‑1 (MQ=255)
            gAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGa                                                    <  1:3305696/68‑1 (MQ=255)
                  ctctCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCt                                              >  1:1295887/1‑68 (MQ=255)
                              gaaCTGAAAGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATc                                 >  1:517654/1‑69 (MQ=255)
                              gaaCTGAAAGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATc                                 >  1:75167/1‑69 (MQ=255)
                              gaaCTGAAAGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATc                                 >  1:2652373/1‑69 (MQ=255)
                              gaaCTGAAAGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATc                                 >  1:153115/1‑69 (MQ=255)
                              gaaCTGAAAGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCCGATc                                 >  1:813588/1‑69 (MQ=255)
                                 cTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATccc                               >  1:2068564/1‑68 (MQ=255)
                                 cTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATccc                               >  1:2203816/1‑68 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:515231/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:3319245/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:2830105/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:3229320/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:3019485/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:2023724/62‑1 (MQ=255)
                                    aaaGGGAAAACCGCAGAGTTTCGTGCGCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGAt                                  <  1:1694466/62‑1 (MQ=255)
                                             aCCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGcc                  >  1:3311630/1‑69 (MQ=255)
                                                             cACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCa   >  1:3368281/1‑68 (MQ=255)
                                                             cACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCa   >  1:2047897/1‑68 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAg                          >  1:604150/1‑44 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAgg     >  1:2288744/1‑65 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCaa  >  1:2375181/1‑68 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCaa  >  1:388374/1‑68 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCaa  >  1:1990508/1‑68 (MQ=255)
                                                              aCGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCaa  >  1:1427289/1‑68 (MQ=255)
                                                              |                                                                   
GAACCGGAGATGGAAAAACTCTCCGACGAAGAACTGAAAGGGAAAACCGCAGAGTTTCGTGCACGTCTGGAAAAAGGCGAAGTGCTGGAAAATCTGATCCCGGAAGCTTTCGCCGTGGTACGTGAGGCAA  >  minE/93112‑93241

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: