Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,217,545 G→A 100% G396G (GGC→GGT)  prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,217,5450GA88.6% 44.3 / ‑1.6 26G396G (GGC→GGT) prccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (23/0);  total (25/1)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAG  >  minE/1217479‑1217604
                                                                  |                                                           
cTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGcc                                                           <  1:760826/69‑1 (MQ=255)
                         aaCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGAc                                   >  1:1529776/1‑68 (MQ=255)
                                                       aaCGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCa            >  1:1411875/1‑61 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGc    >  1:3173286/1‑67 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:43195/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:3295081/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:3210064/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2949555/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2922054/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2891379/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2840885/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:50174/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2542358/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:60318/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2344469/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:220389/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2078323/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:1889207/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:1807302/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:668774/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATCACACTGCTAACATTCTGTTTTTCCAGTTTCTGCaa  >  1:1536760/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATAACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2769537/1‑68 (MQ=255)
                                                         cGCCCCACCACCATTGCTACGCAGGTCGATAATAACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:2765406/1‑68 (MQ=255)
                                                         cGCCCCACAACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:249807/1‑68 (MQ=255)
                                                         cGCCCCACAACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCaa  >  1:1520927/1‑68 (MQ=255)
                                                         cGCCCAACCACCATTGCTACGCAGGTCGATAATCACACTGCTGACATTCTGTTTTTCCAGTTTCTGCa   >  1:1663122/1‑68 (MQ=255)
                                                                  |                                                           
CTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAG  >  minE/1217479‑1217604

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: ‑