Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 355,755 C→G 49.1% S65T (AGC→ACC)  ybcJ ← predicted RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE355,7550CG49.1% ‑2.3 / 63.9 53S65T (AGC→ACC) ybcJpredicted RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (19/8);  new base G (21/5);  total (40/13)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGC  >  minE/355687‑355819
                                                                    |                                                                
cAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGc                                                                  >  1:2803833/1‑69 (MQ=255)
       cTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGAcc                                                           <  1:3352512/69‑1 (MQ=255)
             agagAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                              <  1:2203642/60‑1 (MQ=255)
             agagAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                              <  1:1226111/60‑1 (MQ=255)
             agagAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                              <  1:2187109/60‑1 (MQ=255)
             agagAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                              <  1:3411872/60‑1 (MQ=255)
                  aCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTc                                                 <  1:3017597/68‑1 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:1741944/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:583547/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:2460893/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:3411830/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:338332/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:2052984/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:1518944/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:3043646/1‑60 (MQ=255)
                     ggCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                      >  1:3020728/1‑60 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:2804041/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:2611145/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:1019049/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:273216/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:1766456/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:790417/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:1384519/1‑59 (MQ=255)
                      gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGGTGTGACCTGCaa                                                      >  1:1032077/1‑59 (MQ=255)
                             gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                     >  1:2812640/1‑69 (MQ=255)
                             gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                     >  1:860499/1‑69 (MQ=255)
                             gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                     >  1:2750624/1‑69 (MQ=255)
                             gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                     >  1:254914/1‑69 (MQ=255)
                             gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                     >  1:1635881/1‑69 (MQ=255)
                                aCAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCAGTGTCTGACCGGc                                   <  1:927083/68‑1 (MQ=255)
                                aCAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGc                                   <  1:3214766/68‑1 (MQ=255)
                                           aaGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCa                       >  1:1827342/1‑69 (MQ=255)
                                               gATCAGTCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGt                   <  1:2272604/69‑1 (MQ=255)
                                               gATCAGGCAACAACCTGTACGGTGTGACCTGCAAAAGTCACTGTCTGACCGGCGACAATTTTGCAGCGt                   <  1:2623015/69‑1 (MQ=255)
                                               gATCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGt                   <  1:2862691/69‑1 (MQ=255)
                                               gATCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGt                   <  1:2154306/69‑1 (MQ=255)
                                               gATCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTAACTGTCTGACCGGCGACAATTTTGCAGCGt                   <  1:1996647/69‑1 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGcgac                                >  1:531737/1‑54 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2342259/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:969996/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2462635/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2939289/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:781122/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:670669/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:584122/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:530836/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:455004/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2458692/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2665598/1‑69 (MQ=255)
                                                 tCAGGCAACAACCTGTACGGTGTGACCTGCAAAACTCACTGTCTGACCGGCGACAATTTTGCAGCGttt                 >  1:2018057/1‑69 (MQ=255)
                                                          aacCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTc        >  1:1761742/1‑69 (MQ=255)
                                                          aacCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTc        >  1:922190/1‑69 (MQ=255)
                                                          aacCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTc        >  1:3318142/1‑69 (MQ=255)
                                                                tACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGc  >  1:1192828/1‑69 (MQ=255)
                                                                    |                                                                
CAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGC  >  minE/355687‑355819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: