Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,121 G→T 33.8% S193S (TCG→TCT yjbN → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,1210GT33.8% 68.7 / 42.9 65S193S (TCG→TCTyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (28/15);  new base T (11/11);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCA  >  minE/2801056‑2801180
                                                                 |                                                           
aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                          <  1:2686664/69‑1 (MQ=255)
aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                          <  1:1251096/69‑1 (MQ=255)
aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                          <  1:2456766/69‑1 (MQ=255)
    ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTaa                                                      >  1:3191300/1‑69 (MQ=255)
             gCAAAGGTGAGTGTTAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:1629207/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGtcc                                                  >  1:2139107/1‑64 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:904947/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:3146562/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:2336126/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:8806/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:1483991/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:2029360/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:1680356/1‑68 (MQ=255)
             gCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCgaaa                                              >  1:1598756/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:3106924/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:29659/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:1221878/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:3246985/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:3320278/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:789172/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:2159748/68‑1 (MQ=255)
                    tgAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:837290/68‑1 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:1327805/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:3399036/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:1574665/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:199643/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:1646999/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:1818830/1‑66 (MQ=255)
                          tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:1814857/1‑66 (MQ=255)
                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTACGCCCGAAAGAAAACCGTGAAATCCcgccgc                        >  1:503203/1‑68 (MQ=255)
                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                   >  1:3128552/1‑57 (MQ=255)
                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                        >  1:869426/1‑68 (MQ=255)
                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                        >  1:1189686/1‑68 (MQ=255)
                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                        >  1:2072228/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:4444/1‑69 (MQ=255)
                                       atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                   <  1:894102/69‑1 (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:345671/69‑1 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaaga                                            >  1:2270166/1‑41 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                     >  1:2345056/1‑64 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                     >  1:124327/1‑64 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                     >  1:2506274/1‑64 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                     >  1:1452942/1‑64 (MQ=255)
                                          cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                     >  1:2473255/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:3178197/69‑1 (MQ=255)
                                             gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg             <  1:2763306/69‑1 (MQ=255)
                                             gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg             <  1:1090979/69‑1 (MQ=255)
                                             gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg             <  1:1594360/69‑1 (MQ=255)
                                             gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg             <  1:667225/69‑1 (MQ=255)
                                              cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc            >  1:1745841/1‑69 (MQ=255)
                                              cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc            >  1:301857/1‑69 (MQ=255)
                                              cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc            >  1:2078953/1‑69 (MQ=255)
                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                   <  1:2903919/58‑1 (MQ=255)
                                                     aGCCTGGCTTTCTGGGTTAAGTCCGAAAGAAAACCg                                      <  1:1454968/36‑1 (MQ=38)
                                                      gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa     <  1:1503254/68‑1 (MQ=255)
                                                      gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa     <  1:2505829/68‑1 (MQ=255)
                                                        cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa  <  1:3168713/69‑1 (MQ=255)
                                                                 |                                                           
AATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCA  >  minE/2801056‑2801180

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: