Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,214,266 G→A 58.1% A332V (GCT→GTT)  mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,214,2660GA58.1% 5.4 / 19.9 31A332V (GCT→GTT) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base G (2/11);  new base A (6/12);  total (8/23)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCA  >  minE/2214212‑2214325
                                                      |                                                           
gCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                                                         <  1:901206/59‑1 (MQ=255)
gCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                                                         <  1:2048911/59‑1 (MQ=255)
gCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                                                         <  1:2175407/59‑1 (MQ=255)
gCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                                                         <  1:3416409/59‑1 (MQ=255)
gCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                                                         <  1:324462/59‑1 (MQ=255)
      tgtgGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt                                         >  1:120451/1‑69 (MQ=255)
        tgGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt                                         <  1:2185799/67‑1 (MQ=255)
        tgGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt                                         <  1:2999717/67‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2552050/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:562978/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:400657/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:390920/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:380958/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:338777/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2542801/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2347449/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2055366/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTGGGAGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2380940/53‑1 (MQ=255)
            gagCGATTAGCAGCGGGCTGTGTAGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2292368/53‑1 (MQ=255)
             agCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgctg                                    >  1:2139044/1‑67 (MQ=255)
                    agcagcGGGCTGTGTGGGCGCAGCAGGCGCAGTAACTGCATTTcc                                                   <  1:2689434/45‑1 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:258192/4‑69 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:899512/4‑69 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:2803706/4‑69 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:2244299/4‑69 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:2183825/4‑69 (MQ=255)
                                 gtgggCGCAGCAGGCGCAGTAACTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGTGATcccc              >  1:862252/4‑69 (MQ=255)
                                  tCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTg                        <  1:182082/58‑1 (MQ=255)
                                             ggCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCa  <  1:1348586/69‑1 (MQ=255)
                                             ggCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCa  <  1:2733450/69‑1 (MQ=255)
                                              gagcAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCa  <  1:2455134/66‑1 (MQ=255)
                                                      |                                                           
GCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCA  >  minE/2214212‑2214325

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