Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 634,389 T→C 100% T455T (ACT→ACC mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE634,3890TC93.6% 187.3 / 2.6 62T455T (ACT→ACCmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base C (0/58);  total (4/58)
Fisher's exact test for biased strand distribution p-value = 1.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGT  >  minE/634367‑634447
                      |                                                          
aGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATgg                       >  1:99699/1‑60 (MQ=255)
aGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAAc                >  1:2298926/1‑67 (MQ=255)
aGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAAc                >  1:1504419/1‑67 (MQ=255)
aGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAAc                >  1:134665/1‑67 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCAcagtcagt  <  1:354069/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:545273/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2237819/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:52005/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:451403/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:346212/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:324770/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:294732/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:249311/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2301524/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2258572/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2256487/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2251175/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2159707/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:615020/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:641407/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:642228/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:655770/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:772658/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:778345/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:831453/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:8587/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:861725/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:869477/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:958540/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:962290/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:974367/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:980037/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1682301/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1051309/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1057402/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:10696/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1199106/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1233317/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1262836/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1338706/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1371142/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1378537/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1403448/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1534593/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1597586/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1639675/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2191207/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1740729/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:175286/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1756422/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1766586/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1852808/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1856176/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1862950/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1868049/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1960053/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1966590/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1985907/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:2117670/65‑1 (MQ=255)
                aGCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1013105/65‑1 (MQ=255)
                 gCGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:1413227/64‑1 (MQ=255)
                  cGACCGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGt  <  1:193881/63‑1 (MQ=255)
                      |                                                          
AGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGT  >  minE/634367‑634447

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