Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,053,383 G→A 43.7% G260G (GGG→GGA 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,3830GA43.7% 21.6 / 60.9 80G260G (GGG→GGAccafused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (17/18);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAG  >  minE/2053317‑2053444
                                                                  |                                                             
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:848570/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:3178288/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:3152080/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2960182/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2870916/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2802085/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2794888/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2429855/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2260144/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:2061596/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:1840916/69‑1 (MQ=255)
ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:1745191/69‑1 (MQ=255)
 gcggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCt                                                             <  1:1096241/68‑1 (MQ=255)
                       tCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgc                                      <  1:151279/69‑1 (MQ=255)
                           tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                   <  1:1348057/68‑1 (MQ=255)
                           tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                   <  1:232610/68‑1 (MQ=255)
                           tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                   <  1:2287795/68‑1 (MQ=255)
                           tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                   <  1:464950/68‑1 (MQ=255)
                           tGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtc                                   <  1:725648/68‑1 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:1617521/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:1296212/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:3101815/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:3105890/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:1458505/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:453486/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:800622/1‑52 (MQ=255)
                            gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAg                                                  >  1:1884507/1‑52 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                >  1:778673/1‑69 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                >  1:1539891/1‑69 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                <  1:1178357/69‑1 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                <  1:2203032/69‑1 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                <  1:3178340/69‑1 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                <  1:3257565/69‑1 (MQ=255)
                             tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                >  1:1960090/1‑69 (MQ=255)
                              ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                <  1:1478942/68‑1 (MQ=255)
                               cGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatca                               <  1:3047859/68‑1 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:2616315/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:1964394/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:2266793/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:794124/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:2737251/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:1032969/1‑69 (MQ=255)
                                   tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                          >  1:2801904/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:942928/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:925286/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:3056728/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:601998/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:753797/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2436517/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:211063/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:210378/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2247627/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2302922/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2325875/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2329782/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2440356/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:1745103/1‑67 (MQ=255)
                                                    gATCTCGGTAAAGGACTGACGCCGCCAGAAATTTGGCCACGTCATCATGGTCATGGCCCGGCGGgtg           >  1:2635322/1‑67 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGTCGGGTGTGAAGTTAg  <  1:811403/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:182081/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:1615137/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:1980920/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:1504074/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:180411/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:904269/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:1334773/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:702239/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2118679/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:237321/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:1290842/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2530560/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2839180/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2953406/68‑1 (MQ=255)
                                                            tAAAGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGGGTGAAGTTAg  <  1:2977978/68‑1 (MQ=255)
                                                             aaaGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2399757/67‑1 (MQ=255)
                                                             aaaGGACTGACGCCGCCAGAACTTTGGCCACGTCATCATGGTCATGGCCCGGCGGGTGTGAAGTTAg  <  1:2529298/67‑1 (MQ=255)
                                                                  |                                                             
GGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAG  >  minE/2053317‑2053444

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: