Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,528,315 A→G 38.1% N434D (AAT→GAT)  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,3150AG38.1% 67.5 / 84.5 126N434D (AAT→GAT) trmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base A (45/33);  new base G (18/30);  total (63/63)
Fisher's exact test for biased strand distribution p-value = 4.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 126 total aligned reads displayed.

CCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCG  >  minE/1528251‑1528381
                                                                |                                                                  
ccGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAAttt                                                                <  1:828729/69‑1 (MQ=255)
     aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                                            >  1:2250005/1‑68 (MQ=255)
     aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                                            >  1:2159017/1‑68 (MQ=255)
     aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                                            >  1:1992707/1‑68 (MQ=255)
       gtgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAAttt                                                                <  1:1579323/62‑1 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:2269083/1‑69 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:2177990/1‑69 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:681536/1‑69 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:3250972/1‑69 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:107992/1‑69 (MQ=255)
        tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                                        >  1:1078002/1‑69 (MQ=255)
          cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                                      <  1:1118764/69‑1 (MQ=255)
          cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                                      <  1:2481142/69‑1 (MQ=255)
          cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                                      <  1:1445778/69‑1 (MQ=255)
          cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                                      <  1:1032090/69‑1 (MQ=255)
          cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                                      <  1:259559/69‑1 (MQ=255)
                     tcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1749269/68‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGTATGACGgttg                                           <  1:49674/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:3155835/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:2721915/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:850927/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:351291/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:3385901/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:2777374/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:668802/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1016132/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:3293616/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:684119/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:170733/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:729583/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1092342/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:2523679/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1561268/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:829287/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:3183701/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1396556/69‑1 (MQ=255)
                     gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:857359/69‑1 (MQ=255)
                      cAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:2944819/68‑1 (MQ=255)
                      cAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                                           <  1:1614018/68‑1 (MQ=255)
                                gTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAtt                                >  1:2253574/1‑69 (MQ=255)
                                 tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTgttggt                                         >  1:3058894/1‑59 (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:2805720/1‑69 (MQ=255)
                                 tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                               >  1:14288/1‑69 (MQ=255)
                                 tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                               >  1:1041320/1‑69 (MQ=255)
                                 tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt                               >  1:2687867/1‑69 (MQ=255)
                                 tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAtt                                >  1:1125315/1‑68 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:3016364/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:1020790/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:3146655/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:449307/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:500811/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:593478/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:1797580/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:1785630/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:2726260/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:1308462/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:2595780/1‑63 (MQ=255)
                                   tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa                                   >  1:2535958/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:526412/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:45265/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:771285/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:3242739/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:3139708/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:1266579/69‑1 (MQ=255)
                                         tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag                       <  1:2988506/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:2403210/1‑68 (MQ=255)
                                                      tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa           >  1:844237/1‑68 (MQ=255)
                                                      tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa           >  1:2689304/1‑68 (MQ=255)
                                                      tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa           >  1:2105806/1‑68 (MQ=255)
                                                      tCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa           >  1:1787081/1‑68 (MQ=255)
                                                       cAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCagcaagc             >  1:2950585/1‑65 (MQ=255)
                                                       cAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCagcaagc             >  1:3049438/1‑65 (MQ=255)
                                                        aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca         >  1:2016327/1‑68 (MQ=255)
                                                        aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca         >  1:2023199/1‑68 (MQ=255)
                                                        aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca         >  1:1717920/1‑68 (MQ=255)
                                                        aGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca         >  1:144285/1‑68 (MQ=255)
                                                        aGTTACAGAATTTATCCCCTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAaca         >  1:873632/1‑68 (MQ=255)
                                                           tACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAgcaagcaa           >  1:1655704/1‑63 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:1888685/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:578386/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:2380941/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:496199/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:459742/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:2381773/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:2420597/69‑1 (MQ=255)
                                                            aCAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAg    <  1:3099279/69‑1 (MQ=255)
                                                               gAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCg  <  1:264973/68‑1 (MQ=255)
                                                                |                                                                  
CCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAGATCAGCAAGCAACACACAGCG  >  minE/1528251‑1528381

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: