Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,528,792 T→C 43.6% L593L (TTG→CTG)  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,7920TC43.6% 20.2 / 74.3 94L593L (TTG→CTG) trmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (41/12);  new base C (19/22);  total (60/34)
Fisher's exact test for biased strand distribution p-value = 2.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

ATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTT  >  minE/1528730‑1528859
                                                              |                                                                   
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:2813222/1‑69 (MQ=255)
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:2778441/1‑69 (MQ=255)
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:1137544/1‑69 (MQ=255)
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:2562230/1‑69 (MQ=255)
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:3285504/1‑69 (MQ=255)
atcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCgg                                                               >  1:1742285/1‑69 (MQ=255)
          aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:2922332/67‑1 (MQ=255)
          aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:2929940/67‑1 (MQ=255)
          aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:1795759/67‑1 (MQ=255)
          aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:917263/67‑1 (MQ=255)
            tGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATg                                                   >  1:2632481/1‑69 (MQ=255)
                       gTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgc                                        >  1:1998179/1‑69 (MQ=255)
                        ttCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgc                                        >  1:3033419/1‑68 (MQ=255)
                          cccGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgcgcg                                     >  1:21347/1‑69 (MQ=255)
                          cccGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgcgcg                                     >  1:1225308/1‑69 (MQ=255)
                          ccagATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgcgcg                                     >  1:2248414/4‑69 (MQ=255)
                               ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:1178874/1‑69 (MQ=255)
                               ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:3118411/1‑69 (MQ=255)
                               ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:259956/1‑69 (MQ=255)
                               ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:2926301/1‑69 (MQ=255)
                                tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:1823694/1‑69 (MQ=255)
                                tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:2370080/1‑69 (MQ=255)
                                tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:2435467/1‑69 (MQ=255)
                                tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:409734/1‑69 (MQ=255)
                                tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:1023028/1‑69 (MQ=255)
                                  tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:2926669/1‑41 (MQ=255)
                                  tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:261086/1‑41 (MQ=255)
                                  tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:2396239/1‑41 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACg                        >  1:2617552/2‑68 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:253670/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:1002584/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2938930/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2404594/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2387709/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2669450/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:1030685/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2881515/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:3074853/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2109277/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:60663/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:3185257/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:142136/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:925375/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:1910005/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:33363/2‑69 (MQ=255)
                                        tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGAGCGCCGGt                                >  1:1791879/2‑60 (MQ=255)
                                        tacacTCGTCGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACg                        >  1:2768968/2‑68 (MQ=255)
                                         acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGa                       >  1:1925047/1‑68 (MQ=255)
                                         acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGa                       >  1:1450029/1‑68 (MQ=255)
                                         acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                      >  1:133425/1‑69 (MQ=255)
                                         acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                      >  1:1190555/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATg                                                   >  1:2345800/1‑35 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATctct                  >  1:28483/1‑68 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:388413/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:3273257/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:452346/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:544083/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:574527/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:2888256/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:587700/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:101151/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2453372/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:1803519/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:1872300/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:2038848/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:1563473/1‑69 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2457217/69‑1 (MQ=255)
                                              cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:1470255/69‑1 (MQ=255)
                                               gTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2820441/68‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1517733/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:154134/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:743360/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2822454/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:826959/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:161139/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:881877/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1654586/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1660701/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:177720/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:369607/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1980522/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2748182/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2670233/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:255300/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2469330/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1983571/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:3201118/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:3234009/69‑1 (MQ=255)
                                                  gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:362710/69‑1 (MQ=255)
                                                   aaTATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2125768/68‑1 (MQ=255)
                                                   aaTATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2280702/68‑1 (MQ=255)
                                                       tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTg         >  1:849960/1‑68 (MQ=255)
                                                       tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTg         >  1:2729682/1‑68 (MQ=255)
                                                       tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTg         >  1:1358345/1‑68 (MQ=255)
                                                         cgtcgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCgg      >  1:849133/1‑69 (MQ=255)
                                                             gtTGGCGGAGCAGAAAGATGAGGCGGTATGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCttt  >  1:2722848/1‑69 (MQ=255)
                                                             gtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCttt  >  1:999579/1‑69 (MQ=255)
                                                              |                                                                   
ATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTT  >  minE/1528730‑1528859

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: