Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,051 A→G 81.0% L113L (TTG→CTG) ‡ yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,0510AG81.0% 62.4 / 8.1 47L113L (TTG→CTG) ‡yajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (23/15);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

ACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCG  >  minE/252987‑253098
                                                                |                                               
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1002401/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:490115/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:3194728/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2975837/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2917401/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2855718/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2853650/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2768636/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2589113/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2166070/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:2003731/69‑2 (MQ=255)
accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                                             <  1:1809219/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:1264403/69‑1 (MQ=255)
     ccACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                                        <  1:1528940/69‑1 (MQ=255)
          gCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACtt                                   >  1:1833913/1‑69 (MQ=255)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3172654/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:924416/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3180355/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3189399/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:3295866/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:395572/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:423/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1000458/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:508192/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:545916/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:666062/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:88859/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1258107/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1901592/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1908354/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1916792/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1931129/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2203408/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2588595/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:2727952/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:1615491/1‑63 (MQ=19)
             ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca                                      >  1:136049/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:3421974/1‑68 (MQ=255)
                                            aCGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCg  >  1:1805621/1‑68 (MQ=255)
                                            aCGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCg  >  1:3042297/1‑68 (MQ=255)
                                                                |                                               
ACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCG  >  minE/252987‑253098

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