Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,858,656 T→C 100% G310G (GGT→GGC 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,6560TC82.9% 33.0 / ‑0.1 23G310G (GGT→GGCamiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (0/19);  total (1/22)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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: