Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 636,766 A→G 100% I1248V (ATT→GTT)  mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE636,7660AG90.9% 28.1 / ‑3.1 11I1248V (ATT→GTT) mukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/0);  total (11/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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.

CCCGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCATTCAGCGCGAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGG  >  minE/636707‑636827
                                                           |                                                             
cccGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCa                                                               >  1:137348/1‑60 (MQ=255)
cccGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCa                                                               >  1:567955/1‑60 (MQ=255)
cccGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCa                                                               >  1:366250/1‑60 (MQ=255)
cccGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCa                                                               >  1:304686/1‑60 (MQ=255)
cccGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCa                                                               >  1:248731/1‑60 (MQ=255)
                      cAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCATTCAGCGCGAGCAGAAccg                                            >  1:11891/1‑57 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:169819/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:309341/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:312511/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:316189/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:137688/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:398821/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:447948/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:506222/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:547837/1‑62 (MQ=255)
                                  cGTGGCGAACATCATTCGCAAAACCGTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAAcc                           >  1:567424/1‑62 (MQ=255)
                                                           aTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTgg  <  1:261681/62‑1 (MQ=255)
                                                           aTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTgg  <  1:383697/62‑1 (MQ=255)
                                                           aTTCAGCGCGAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTgg  <  1:453827/62‑1 (MQ=255)
                                                           |                                                             
CCCGTGAACAGAAACTGGCGATCAGTTCCCGCAGCGTGGCGAACATCATTCGCAAAACCATTCAGCGCGAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGG  >  minE/636707‑636827

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