Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,528,801 C→A 44.5% Q596K (CAG→AAG)  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,8010CA44.5% 40.4 / 88.8 92Q596K (CAG→AAG) trmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base C (38/13);  new base A (19/22);  total (57/35)
Fisher's exact test for biased strand distribution p-value = 9.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTT  >  minE/1528740‑1528859
                                                             |                                                          
aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:2922332/67‑1 (MQ=255)
aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:1795759/67‑1 (MQ=255)
aaTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGaaa                                                       <  1:2929940/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:1225308/1‑69 (MQ=255)
                cccGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgcgcg                                     >  1:21347/1‑69 (MQ=255)
                ccagATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAgcgcg                                     >  1:2248414/4‑69 (MQ=255)
                     ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:259956/1‑69 (MQ=255)
                     ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:3118411/1‑69 (MQ=255)
                     ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:1178874/1‑69 (MQ=255)
                     ttATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGt                                >  1:2926301/1‑69 (MQ=255)
                      tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:409734/1‑69 (MQ=255)
                      tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:1023028/1‑69 (MQ=255)
                      tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:1823694/1‑69 (MQ=255)
                      tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:2435467/1‑69 (MQ=255)
                      tATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGtt                               >  1:2370080/1‑69 (MQ=255)
                        tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:261086/1‑41 (MQ=255)
                        tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:2926669/1‑41 (MQ=255)
                        tGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGa                                                         >  1:2396239/1‑41 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACg                        >  1:2617552/2‑68 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2404594/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2387709/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:33363/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2669450/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:253670/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:1002584/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:1910005/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:1030685/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:925375/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:142136/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:3074853/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2938930/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGa                       >  1:2881515/2‑69 (MQ=255)
                              tacacTCGTGGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGAGCGCCGGt                                >  1:1791879/2‑60 (MQ=255)
                              tacacTCGTCGAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACg                        >  1:2768968/2‑68 (MQ=255)
                               acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGa                       >  1:1450029/1‑68 (MQ=255)
                               acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGa                       >  1:1925047/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:2888256/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:388413/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:3273257/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:544083/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:101151/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:1470255/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:1563473/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:587700/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:574527/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:1803519/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:1872300/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:452346/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 >  1:2038848/1‑69 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2453372/69‑1 (MQ=255)
                                    cGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2457217/69‑1 (MQ=255)
                                     gTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:2820441/68‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:3201118/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:369607/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:3234009/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:255300/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:743360/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:826959/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:881877/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:362710/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2748182/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1517733/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:154134/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:161139/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:1980522/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:1983571/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2469330/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2670233/69‑1 (MQ=255)
                                        gAATATGCGTCGCTGGCGGAGAAGAAAGATGAAGCGGTAAGCGCGCCGGTTTATGACGATCTCTTTATg             <  1:2822454/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)
                                                        cGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCtt                 <  1:1667721/49‑1 (MQ=255)
                                                             |                                                          
AATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTT  >  minE/1528740‑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: