Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,528,336 T→G 43.7% C441G (TGT→GGT)  trmC → fused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,528,3360TG43.7% 19.6 / 122.7 110C441G (TGT→GGT) trmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (38/24);  new base G (18/30);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 2.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAGTGGTACTGGCGA  >  minE/1528272‑1528396
                                                                |                                                            
tcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1749269/68‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGTATGACGgttg                                                          <  1:49674/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1396556/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:460367/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:3155835/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:668802/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:351291/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:2523679/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:3293616/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1561268/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:684119/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:3183701/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1364319/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:729583/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:829287/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:2549034/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:2777374/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1120028/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1092342/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:850927/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:857359/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:2721915/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1016132/69‑1 (MQ=255)
gcAACAGCAGGGTTTGCAGATTCATTATCAACATCACTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1840711/69‑1 (MQ=255)
 cAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:1614018/68‑1 (MQ=255)
 cAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:2944819/68‑1 (MQ=255)
           gTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAtt                                               >  1:2253574/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGaa                                                 >  1:234649/1‑66 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:3055663/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:2687867/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:2805720/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:2153109/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:2134056/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:55845/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:1752284/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:14288/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                                              >  1:1041320/1‑69 (MQ=255)
            tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAtt                                               >  1:1125315/1‑68 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:3269554/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:449307/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:3146655/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:593478/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:3016364/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:671871/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:674221/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:829882/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:1020790/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:1785630/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:2595780/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:2535958/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:1797580/1‑63 (MQ=255)
              tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                                  >  1:1334801/1‑63 (MQ=255)
                cAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGc                                           >  1:17011/1‑68 (MQ=255)
                 aGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGc                                           >  1:3426675/1‑67 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:603113/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:1696796/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:3242739/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:3139708/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:771285/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:45265/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:1266579/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:2988506/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:526412/69‑1 (MQ=255)
                    tttATTATCAATATCAGTTAAAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:676237/69‑1 (MQ=255)
                    tttATTATCAATATCAGATACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                                      <  1:2503105/69‑1 (MQ=255)
                                 tCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                                          <  1:3034704/36‑1 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:2689304/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:642622/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:1787081/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:2105806/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:3152465/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:2403210/1‑68 (MQ=255)
                                 tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:844237/1‑68 (MQ=255)
                                  cAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCagcaagc                            >  1:3049438/1‑65 (MQ=255)
                                  cAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCagcaagc                            >  1:2950585/1‑65 (MQ=255)
                                  cAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCagcaagc                            >  1:1062262/1‑65 (MQ=255)
                                   aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:2023199/1‑68 (MQ=255)
                                   aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:2163134/1‑68 (MQ=255)
                                   aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:2016327/1‑68 (MQ=255)
                                   aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:144285/1‑68 (MQ=255)
                                   aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:1717920/1‑68 (MQ=255)
                                   aGTTACAGAATTTATCCCCTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca                        >  1:873632/1‑68 (MQ=255)
                                      tACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa                          >  1:1655704/1‑63 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAAGACACAg                   <  1:436845/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:2381773/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:2668562/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:2965787/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:2420597/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:792767/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:3099279/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:578386/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:2380941/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:1888685/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:459742/69‑1 (MQ=255)
                                       aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg                   <  1:496199/69‑1 (MQ=255)
                                          gAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCg                 <  1:264973/68‑1 (MQ=255)
                                             tttATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAGt             >  1:328164/1‑69 (MQ=255)
                                                        aGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAg              >  1:1078504/1‑57 (MQ=255)
                                                        aGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAGTGGTACTGGCGa  >  1:3385944/1‑69 (MQ=255)
                                                        aGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAGTGGTACTGGCGa  >  1:1157562/1‑69 (MQ=255)
                                                                |                                                            
GCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCGTAGTGGTACTGGCGA  >  minE/1528272‑1528396

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: