Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,528,294 T→C 52.8% Y427H (TAT→CAT)  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,2940TC52.8% 8.6 / 75.9 89Y427H (TAT→CAT) trmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (32/10);  new base C (18/29);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 5.40e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 111 total aligned reads displayed.

GTTGGCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAG  >  minE/1528226‑1528360
                                                                    |                                                                  
gtTGGCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGAttt                                                                    >  1:728205/1‑69 (MQ=255)
gtTGGCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGAttt                                                                    >  1:3276679/1‑69 (MQ=255)
gtTGGCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGAttt                                                                    >  1:825276/1‑69 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCATGGTTTGCAGATttatt                                                                 <  1:1340622/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:1332541/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:610732/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:215414/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:831007/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:1097251/69‑1 (MQ=255)
   ggCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:2646534/69‑1 (MQ=255)
            cagcagAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATc                                                        <  1:710814/69‑1 (MQ=255)
                agAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGt                                                     >  1:799877/1‑68 (MQ=255)
                agAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGt                                                     >  1:328500/1‑68 (MQ=255)
                agAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGt                                                     >  1:714566/1‑68 (MQ=255)
                         ccGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAAttt                                           <  1:845198/69‑1 (MQ=255)
                         ccGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAAttt                                           <  1:828729/69‑1 (MQ=255)
                           gTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATttatt                                                                 <  1:2323515/45‑1 (MQ=255)
                              aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                       >  1:2159017/1‑68 (MQ=255)
                              aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                       >  1:2250005/1‑68 (MQ=255)
                              aTGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATcc                                       >  1:1992707/1‑68 (MQ=255)
                                gtgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAAttt                                           <  1:1579323/62‑1 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:681536/1‑69 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:2269083/1‑69 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:526431/1‑69 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:1078002/1‑69 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:107992/1‑69 (MQ=255)
                                 tgCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTa                                   >  1:1542790/1‑69 (MQ=255)
                                   cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                 <  1:2743004/69‑1 (MQ=255)
                                   cTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAg                                 <  1:1032090/69‑1 (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:259559/69‑1 (MQ=255)
                                              tcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1749269/68‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGTATGACGgttg                      <  1:49674/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:857359/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:2777374/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:3155835/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1016132/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:3183701/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:2721915/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:850927/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:829287/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:729583/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:3293616/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:3385901/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:351291/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:684119/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:668802/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1396556/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1364319/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1092342/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:2523679/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:170733/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1120028/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1561268/69‑1 (MQ=255)
                                              gcAACAGCAGGGTTTGCAGATTCATTATCAACATCACTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1840711/69‑1 (MQ=255)
                                               cAACAGCAGGGTTTGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGgttg                      <  1:1614018/68‑1 (MQ=255)
                                                         gTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAtt           >  1:2253574/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTgttggt                    >  1:3058894/1‑59 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGaa             >  1:234649/1‑66 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:2134056/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:1752284/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:2687867/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:2153109/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:55845/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:14288/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:3055663/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:2805720/1‑69 (MQ=255)
                                                          tttGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAAttt          >  1:1041320/1‑69 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:1334801/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:1308462/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:829882/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:674221/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:1020790/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:593478/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:500811/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:2050544/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:2595780/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:2726260/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:2535958/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:3016364/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:3146655/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:3269554/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:1785630/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:449307/1‑63 (MQ=255)
                                                            tGCAGATTCATTATCAACATCAGTTACAGGATTTATCCCGTAAGGATGACGGTTGGTTGTTGa              >  1:1797580/1‑63 (MQ=255)
                                                              cAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGc       >  1:17011/1‑68 (MQ=255)
                                                               aGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGc       >  1:3426675/1‑67 (MQ=255)
                                                               aGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGc       >  1:68141/1‑67 (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:1266579/69‑1 (MQ=255)
                                                                  tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag  <  1:45265/69‑1 (MQ=255)
                                                                  tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag  <  1:2988506/69‑1 (MQ=255)
                                                                  tttATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGgag  <  1:603113/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)
                                                                    |                                                                  
GTTGGCTGTGCCCAGCAGAACTGACCCGTAATGTGCTGGAACTGGCGCAACAGCAGGGTTTGCAGATTTATTATCAATATCAGTTACAGAATTTATCCCGTAAGGATGACTGTTGGTTGTTGAATTTTGCAGGAG  >  minE/1528226‑1528360

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: