Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,130 C→T 25.4% S196S (AGC→AGT yjbN → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,1300CT25.4% 100.4 / 27.7 82S196S (AGC→AGTyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (39/22);  new base T (10/11);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTT  >  minE/2801069‑2801195
                                                             |                                                                 
gCAAAGGTGAGTGTTAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:1629207/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGtcc                                                                 >  1:2139107/1‑64 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:904947/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:2336126/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:2029360/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:3146562/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:1680356/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:1598756/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:1483991/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                                             >  1:8806/1‑68 (MQ=255)
gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTCCTGGGTTAAGTCCgaaa                                                             >  1:1541880/1‑68 (MQ=255)
       tgAGTGTGAGATTTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:1974096/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGTGTTAAGTCCGAAAGAAAACCg                                                     <  1:812852/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:29659/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:789172/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:3106924/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:2159748/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:1221878/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:837290/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:3246985/68‑1 (MQ=255)
       tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:3320278/68‑1 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:1818830/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:199643/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:1814857/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:1646999/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:1327805/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:3399036/1‑66 (MQ=255)
             tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:1574665/1‑66 (MQ=255)
                      catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTACGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       >  1:503203/1‑68 (MQ=255)
                      catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                  >  1:3128552/1‑57 (MQ=255)
                      catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       >  1:2072228/1‑68 (MQ=255)
                      catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       >  1:1189686/1‑68 (MQ=255)
                      catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       >  1:869426/1‑68 (MQ=255)
                      catcatCAACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       >  1:3330604/1‑68 (MQ=255)
                         catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                       <  1:3076096/65‑1 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  >  1:1298103/1‑69 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  <  1:2716586/69‑1 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  <  1:235353/69‑1 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  >  1:2336674/1‑69 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  >  1:4444/1‑69 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  <  1:345671/69‑1 (MQ=255)
                          atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  <  1:894102/69‑1 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaaga                                                           >  1:2270166/1‑41 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:2345056/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:2506274/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:2473255/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:124327/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:1452942/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:1529714/1‑64 (MQ=255)
                             cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                    >  1:1689932/1‑64 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:1146362/69‑1 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:2763306/69‑1 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:667225/69‑1 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:1090979/69‑1 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:1594360/69‑1 (MQ=255)
                                gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                            <  1:3178197/69‑1 (MQ=255)
                                 cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc                           >  1:1745841/1‑69 (MQ=255)
                                 cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc                           >  1:2078953/1‑69 (MQ=255)
                                 cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc                           >  1:301857/1‑69 (MQ=255)
                                     tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                  <  1:2903919/58‑1 (MQ=255)
                                        aGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                                     <  1:1454968/36‑1 (MQ=38)
                                         gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa                    <  1:2505829/68‑1 (MQ=255)
                                         gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa                    <  1:1503254/68‑1 (MQ=255)
                                           cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                 <  1:3168713/69‑1 (MQ=255)
                                                      ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt       <  1:1473910/68‑1 (MQ=255)
                                                      ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt       <  1:2485308/68‑1 (MQ=255)
                                                      ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg      <  1:1616277/69‑1 (MQ=255)
                                                       gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg      <  1:1432987/68‑1 (MQ=255)
                                                       gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg      <  1:465513/68‑1 (MQ=255)
                                                         tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg      <  1:515709/66‑1 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGATATCCCGCCGCTCGAt                                  <  1:3296632/37‑1 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:1551725/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:2375659/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:2624876/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:714008/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:2875394/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:358763/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:3386255/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:228269/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:1767680/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:2258923/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:3150335/1‑69 (MQ=255)
                                                          aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt  >  1:1938248/1‑69 (MQ=255)
                                                             |                                                                 
GCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTT  >  minE/2801069‑2801195

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: