Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,756,156 A→G 35.5% P221P (CCT→CCC rsmB ← 16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,756,1560AG35.5% 70.1 / 64.6 107P221P (CCT→CCCrsmB16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent
Reads supporting (aligned to +/- strand):  ref base A (29/40);  new base G (5/33);  total (34/73)
Fisher's exact test for biased strand distribution p-value = 2.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 108 total aligned reads displayed.

GCAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTAATCCG  >  minE/2756091‑2756211
                                                                 |                                                       
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:1000503/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:2865107/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:2717053/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:3326545/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:357468/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:515749/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:616172/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:1458389/69‑1 (MQ=255)
gcAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGt                                                      <  1:1439190/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2625146/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:3431311/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2097613/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:870370/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:328940/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:3240311/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2228095/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2474129/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:3408238/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2698562/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:881964/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2735265/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2979640/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2952755/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:293502/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2919001/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1352984/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:758492/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:149872/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:671936/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1547547/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1572879/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1596225/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1600073/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:579142/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1353340/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1669913/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1711707/69‑1 (MQ=255)
     ccAGGTCATGCAACCTTGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:3436124/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:54247/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:3247305/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:3199868/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:799998/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:3259717/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:3076031/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2990371/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:635697/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:387506/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:342605/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:1020285/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2058731/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2061969/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2763333/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2650459/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:1586368/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2376420/69‑1 (MQ=255)
        ggTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:22202/69‑1 (MQ=255)
         gTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:1500834/68‑1 (MQ=255)
           cATGCACCCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAc                                           >  1:924900/1‑69 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCa                                                         >  1:2980336/1‑55 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:1699539/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:1618064/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:631225/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:2165489/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:651231/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:3101814/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:753999/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:773759/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:78954/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:1385948/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:844548/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:922157/1‑56 (MQ=255)
           cATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAg                                                        >  1:2930941/1‑56 (MQ=255)
             tGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAc                                           <  1:2872049/67‑1 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:964276/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:930280/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:2824027/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:1860649/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:707373/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:2207503/1‑68 (MQ=255)
               cAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcagg                                        >  1:3378746/1‑68 (MQ=255)
                 aCCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtg                                      <  1:1461935/68‑1 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:2100844/1‑69 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:1600697/1‑69 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:1964002/1‑69 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:1623866/1‑69 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:3029681/1‑69 (MQ=255)
                    ttGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAAcaggtgcagg                                  >  1:1285766/1‑69 (MQ=255)
                    ttGTGCTGAAGCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:2094889/54‑1 (MQ=255)
                             agCATCCTGAACGGTAACCCATCCTTCTTCAAAACCGGGTAgcgc                                                 <  1:1346512/44‑1 (MQ=38)
                              gCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGa                        <  1:2891188/69‑1 (MQ=255)
                                      aaCGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATg                                              <  1:2206233/39‑1 (MQ=255)
                                               ccATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTa       <  1:2951591/69‑1 (MQ=255)
                                               ccATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTa       <  1:2814446/69‑1 (MQ=255)
                                               ccATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTa       <  1:1642125/69‑1 (MQ=255)
                                               ccATCCGTCTTCAAAACCAGGTAGAGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTa       <  1:944062/69‑1 (MQ=255)
                                               ccATCCGTCTTCAAAACCAGGTAGAGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTa       <  1:962060/69‑1 (MQ=255)
                                                    cttCTTCAAAACCGGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGCTGTACAGCATCAGGGTAATCCg  >  1:869287/3‑69 (MQ=255)
                                                    cttCTTCAAAACCGGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGCTGTACAGCATCAGGGTAATCCg  >  1:515811/3‑69 (MQ=255)
                                                    cttCTTCAAAACCGGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGCTGTACAGCATCAGGGTAATCCg  >  1:3047243/3‑69 (MQ=255)
                                                    cttCTTCAAAACCGGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGCTGTACAGCATCAGGGTAATCCg  >  1:2644426/3‑69 (MQ=255)
                                                                 |                                                       
GCAAGCCAGGTCATGCAACCTTGTGCTGATGCATCCTGAACGGTAACCCATCCGTCTTCAAAACCAGGTAGCGCATGAACAGGTGCAGGTGTTTCCAGACGTACAGCATCAGGGTAATCCG  >  minE/2756091‑2756211

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: