Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 26,138 G→A 100% E29E (GAG→GAA dapB → dihydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE26,1380GA81.4% 18.9 / ‑0.9 16E29E (GAG→GAAdapBdihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (0/13);  total (2/14)
Fisher's exact test for biased strand distribution p-value = 2.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTT  >  minE/26095‑26181
                                           |                                           
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:1279149/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:1741985/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:1774790/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:2168348/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:2502928/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:2893762/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:3046238/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:3324599/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:3330485/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:357367/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:620110/67‑1 (MQ=255)
gccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:869592/67‑1 (MQ=255)
 ccGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAAGGCGTGCAGCTGGGCGCGGCGCTgg                    <  1:467465/67‑1 (MQ=255)
      tGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCg                >  1:513164/1‑67 (MQ=255)
        ggCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGa             >  1:3416713/1‑68 (MQ=255)
                  tGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATcttctt  <  1:56225/69‑1 (MQ=255)
                                           |                                           
GGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTT  >  minE/26095‑26181

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: