| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 978,430 | A→G | 100% | E561E (GAA→GAG) | mukB → | fused chromosome partitioning proteins |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 978,430 | 0 | A | G | 90.6% | 95.6 / 1.5 | 32 | E561E (GAA→GAG) | mukB | fused chromosome partitioning proteins |
| Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (12/17); total (14/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.68e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCAGATTTCTGCAAACGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATGG > W3110S.gb/978365‑978491 | gCAGATTTCTGCAAACGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTgg < 1:1616409/70‑1 (MQ=255)gCAGATTTCTGCAAACGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTgg < 1:216922/70‑1 (MQ=255) aCGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCct < 1:550185/71‑1 (MQ=255) aCGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCct > 1:2370061/1‑71 (MQ=255) cGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCctc > 1:3364450/1‑71 (MQ=255) cGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCctc > 1:3126431/1‑71 (MQ=255) cGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCctc > 1:3507479/1‑71 (MQ=255) cGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCctc > 1:1849910/1‑71 (MQ=255) gCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTcc < 1:2765371/71‑1 (MQ=255) gCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTcc < 1:2224547/71‑1 (MQ=255) aTTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGc > 1:678118/1‑50 (MQ=255) aTTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGc > 1:3371635/1‑71 (MQ=255) aTTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGc > 1:1396006/1‑71 (MQ=255) ttttGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCAtt < 1:104225/52‑1 (MQ=255) ttttGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCAtt < 1:145266/52‑1 (MQ=255) ttttGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCAtt < 1:1359673/52‑1 (MQ=255) ttttGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCAtt < 1:3321679/52‑1 (MQ=255) ttttGATATCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCAtt < 1:1189543/52‑1 (MQ=255) atatCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt < 1:1391065/71‑1 (MQ=255) atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt < 1:3711262/71‑1 (MQ=255) tatCGACGAACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt < 1:2173426/70‑1 (MQ=255) gacgaACTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGc > 1:2375269/1‑71 (MQ=255) aCTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTg < 1:2021641/71‑1 (MQ=255) aCTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTg < 1:2020625/71‑1 (MQ=255) aCTGGAAGCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTg < 1:57136/71‑1 (MQ=255) gCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg < 1:2319524/71‑1 (MQ=255) gCCCTGCATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg < 1:3102407/71‑1 (MQ=255) cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg > 1:2078237/1‑70 (MQ=255) cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg > 1:2010736/1‑70 (MQ=255) cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg > 1:1247997/1‑70 (MQ=255) cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg > 1:821482/1‑70 (MQ=255) cATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg > 1:262471/1‑70 (MQ=255) | GCAGATTTCTGCAAACGTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATGG > W3110S.gb/978365‑978491 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |