Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,031 A→G 79.8% R119R (CGT→CGC yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,0310AG79.8% 62.7 / 9.4 49R119R (CGT→CGCyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/8);  new base G (24/15);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 3.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCG  >  minE/252965‑253098
                                                                  |                                                                   
gAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATaaaatg                                                                   >  1:2472950/1‑67 (MQ=255)
gAAATCCATATCGCGATAGCGCACAAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACg                                                                   >  1:1799779/1‑69 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1002401/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2166070/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2003731/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2589113/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2768636/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2853650/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2855718/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1809219/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2917401/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2975837/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:3194728/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:490115/69‑2 (MQ=255)
                      accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1254750/69‑2 (MQ=255)
                       ccAACCACTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1763570/68‑2 (MQ=255)
                       ccAACCACTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2429672/68‑2 (MQ=255)
                           ccACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                                        <  1:1528940/69‑1 (MQ=255)
                           ccACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                                        <  1:1264403/69‑1 (MQ=255)
                                gCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACtt                                   >  1:1833913/1‑69 (MQ=255)
                                 ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                                               <  1:1432341/40‑1 (MQ=255)
                                 ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                                               <  1:2079556/40‑1 (MQ=255)
                                 ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                                               <  1:2318949/40‑1 (MQ=255)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCg                                                              >  1:2431560/1‑39 (MQ=18)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3180355/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:924416/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:88859/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:666062/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:545916/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:508192/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1258107/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:423/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:395572/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3295866/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:136049/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3189399/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2203408/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3172654/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1615491/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1901592/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1908354/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2727952/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1916792/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2588595/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1931129/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1000458/1‑63 (MQ=19)
                                   ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatc                                       >  1:920606/1‑62 (MQ=19)
                                         acaAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCggg                          <  1:3397077/69‑1 (MQ=255)
                                           aaGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCggg                          <  1:1818901/67‑1 (MQ=255)
                                                           tGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                                        <  1:193841/37‑1 (MQ=255)
                                                                  aCGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCg  >  1:3042297/1‑68 (MQ=255)
                                                                  aCGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCg  >  1:1805621/1‑68 (MQ=255)
                                                                  aCGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCg  >  1:3421974/1‑68 (MQ=255)
                                                                  |                                                                   
GAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCG  >  minE/252965‑253098

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: