Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,053,407 G→A 52.1% P268P (CCG→CCA cca → fused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,053,4070GA52.1% 0.1 / 62.1 71P268P (CCG→CCAccafused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base A (17/20);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGC  >  minE/2053340‑2053462
                                                                   |                                                       
tCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgc                                                        <  1:151279/69‑1 (MQ=255)
    tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                                     <  1:232610/68‑1 (MQ=255)
    tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                                     <  1:2287795/68‑1 (MQ=255)
    tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                                     <  1:1348057/68‑1 (MQ=255)
    tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                                     <  1:464950/68‑1 (MQ=255)
    tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                                     <  1:725648/68‑1 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  <  1:3257565/69‑1 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  <  1:3178340/69‑1 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  >  1:1960090/1‑69 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  <  1:2203032/69‑1 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  >  1:1539891/1‑69 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  >  1:778673/1‑69 (MQ=255)
      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  <  1:1178357/69‑1 (MQ=255)
       ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                  <  1:1478942/68‑1 (MQ=255)
        cGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatca                                                 <  1:3047859/68‑1 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:2737251/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:2616315/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:2801904/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:2266793/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:1032969/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:1964394/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                            >  1:794124/1‑69 (MQ=255)
            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                             <  1:2050657/68‑1 (MQ=255)
                    ttATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcg                                    >  1:3294204/1‑69 (MQ=255)
                             gATCTCTGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:441068/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGtcat                                                   >  1:1349049/1‑45 (MQ=37)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:3056728/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:601998/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:753797/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:925286/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2440356/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2436517/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:942928/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2329782/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:210378/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:211063/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2247627/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:1745103/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2302922/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2325875/1‑67 (MQ=255)
                             gATCTCGGTAAAGGACTGACGCCGCCAGAAATTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2635322/1‑67 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGTCGGGTGTGAAGTTAg                    <  1:811403/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:702239/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1504074/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:180411/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1334773/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:182081/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1980920/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:904269/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1615137/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1290842/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2953406/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2839180/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2118679/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2530560/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:237321/68‑1 (MQ=255)
                                     tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGGGTGAAGTTAg                    <  1:2977978/68‑1 (MQ=255)
                                      aaaGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2529298/67‑1 (MQ=255)
                                      aaaGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2399757/67‑1 (MQ=255)
                                             tGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                             <  1:2445201/35‑1 (MQ=255)
                                             tGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac           <  1:1761707/69‑1 (MQ=255)
                                                 gccgccAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:2037406/56‑1 (MQ=38)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:922687/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:898440/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:1016192/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:715245/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:3255664/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:3204799/1‑47 (MQ=255)
                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                             >  1:2056496/1‑47 (MQ=255)
                                                      cAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGc  <  1:345814/69‑1 (MQ=255)
                                                      cAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGc  <  1:2603323/69‑1 (MQ=255)
                                                               ggCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg                    <  1:1678414/42‑1 (MQ=37)
                                                                   |                                                       
TCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGC  >  minE/2053340‑2053462

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: