Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,480,838 C→A 18.8% L292L (CTG→CTT yidZ ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,480,8380CA18.8% 116.3 / 14.4 58L292L (CTG→CTTyidZpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/14);  new base A (5/6);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TAAAGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAAGA  >  minE/2480770‑2480895
                                                                    |                                                         
tAAAGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagc                                                           <  1:730213/69‑1 (MQ=255)
tAAAGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagc                                                           <  1:2282562/69‑1 (MQ=255)
 aaaGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagc                                                           <  1:1892778/68‑1 (MQ=255)
  aaGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagcag                                                         >  1:2516515/1‑69 (MQ=255)
  aaGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagcag                                                         >  1:2546268/1‑69 (MQ=255)
  aaGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCAcagcag                                                         >  1:1791906/1‑69 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGGCACAGCAGGGTAAAAGGa                                               <  1:2413912/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:1179073/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:2531365/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:2781006/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:286105/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:2215445/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:2194652/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:1995002/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:1697771/68‑1 (MQ=255)
             tGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:1222913/68‑1 (MQ=255)
              ggTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:684052/67‑1 (MQ=255)
              ggTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGa                                               <  1:2715292/67‑1 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:3380724/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:3304162/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:3407761/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         <  1:288195/62‑1 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:361470/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:524699/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:551465/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         <  1:826629/62‑1 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:922312/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:270516/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:2372326/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:1669578/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:1673496/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:2335817/1‑62 (MQ=255)
                         aCCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTcc                                         >  1:1914955/1‑62 (MQ=255)
                              aCGATCTTCGGATTACGGCTGTTCCGTTTATGCCACAGAAGGGTAAAAgg                                                >  1:2512016/1‑50 (MQ=255)
                              aCGATCTTCGGATTACGGCTGTTCCGTTTATGCCACAGAAGGGTAAAAgg                                                >  1:715453/1‑50 (MQ=255)
                              aCGATCTTCGGATTACGGCTGTTCCGTTTATGCCACAGAAGGGTAAAAgg                                                >  1:1637372/1‑50 (MQ=255)
                              aCGATCTTCGGATTACGGCTGTTCCGTTTATGCCACAGAAGGGTAAAAgg                                                >  1:895086/1‑50 (MQ=255)
                              aCGATCTTCGGATTACGGCTGTTCCGTTTATGCCACAGAAGGGTAAAAgg                                                >  1:1915919/1‑50 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:2156361/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:872210/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:3309059/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:2785576/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:2781040/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:1799434/1‑69 (MQ=255)
                                          ttATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:2343468/1‑69 (MQ=255)
                                           tacggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:2944905/66‑1 (MQ=255)
                                           tacggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:1228365/66‑1 (MQ=255)
                                           tacggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:948519/66‑1 (MQ=255)
                                           tacggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:2640301/66‑1 (MQ=255)
                                           tacggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:697212/66‑1 (MQ=255)
                                            acggCTGTTCCGTTTATGCCACAGAAGGGTAAAAGGGACTTCCAGCTTTTTTTGCTGGCTTTCATCaa                <  1:704574/66‑1 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:1475232/1‑58 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:522805/1‑58 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:3142864/1‑58 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCa                 >  1:911002/1‑58 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGgaga        >  1:984685/1‑67 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGgaga        >  1:1891406/1‑67 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGgaga        >  1:3039990/1‑67 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGgaga        >  1:1122028/1‑67 (MQ=255)
                                                     ccGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGgaga        >  1:2495395/1‑67 (MQ=255)
                                                          tATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAaga  >  1:810857/1‑68 (MQ=255)
                                                          tATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAaga  >  1:1540132/1‑68 (MQ=255)
                                                          tATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAaga  >  1:917385/1‑68 (MQ=255)
                                                          tATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAaga  >  1:2240524/1‑68 (MQ=255)
                                                                    |                                                         
TAAAGGTTTTTAATGGTTTCCCGTAACCAGACGATCTTCGGATTATGGCTGTTCCGTTTATGCCACAGCAGGGTAAAAGGAACTTCCAGCTTTTTTTGCTGGCTTTCGTCAAACGGGAGAGGAAGA  >  minE/2480770‑2480895

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