Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,858,650 G→A 100% L308L (CTG→CTA amiB → N‑acetylmuramoyl‑l‑alanine amidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,858,6500GA86.5% 34.2 / ‑1.3 22L308L (CTG→CTAamiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (0/19);  total (0/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCT  >  minE/2858593‑2858696
                                                         |                                              
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:2821696/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:572866/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:523807/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:696578/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:719865/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:3280250/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:3118653/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:3038242/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:3015509/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:2465831/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:730284/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:201095/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:738665/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:1611630/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:1228625/69‑1 (MQ=255)
tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:118169/69‑1 (MQ=255)
       aaCAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:2068909/62‑1 (MQ=255)
                          ggctggAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:860674/43‑1 (MQ=37)
                          ggctggAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc                                     <  1:1663768/43‑1 (MQ=37)
                             tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCt         <  1:2227231/68‑1 (MQ=255)
                             tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCt         <  1:1192971/68‑1 (MQ=255)
                             tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGGGATGTGCTGGCGAACAGTCAGTCt         <  1:399418/68‑1 (MQ=255)
                                                      gctgctGGGTGGGGCTGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCt  >  1:368973/1‑50 (MQ=255)
                                                         |                                              
TCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCT  >  minE/2858593‑2858696

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: