| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,824,533:1 | +G | 100% | pseudogene (666/687 nt) | rph → | ECK3633:JW3618:b3643; ribonuclease PH |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,824,533 | 1 | . | G | 96.6% | 102.3 / ‑3.1 | 29 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (9/19); total (9/20) | |||||||||||
| 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: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
GGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824483‑3824597 | ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg > 1:508068/1‑70 (MQ=255)ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg > 1:1122289/1‑70 (MQ=255)ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg > 1:1468248/1‑70 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAg > 1:898460/1‑70 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAg > 1:1223769/1‑70 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAg > 1:1378064/1‑70 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAggcggc < 1:1635870/70‑1 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAggcggc < 1:36103/70‑1 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAggcggc < 1:1859102/70‑1 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAggcggc > 1:882222/1‑70 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg > 1:1843286/1‑70 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc < 1:1554768/71‑1 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc < 1:1973190/71‑1 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc < 1:888908/71‑1 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc < 1:709106/71‑1 (MQ=255) tACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc < 1:1816501/70‑1 (MQ=255) ttgGCGCTGGCCCGA‑GGGGTAATCGAATCCATTGTAGCGACGCAGAAggcggc < 1:1857188/53‑1 (MQ=255) tgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc > 1:80371/1‑54 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGTCAAACTGATTTTTAAGGCGACt < 1:1556225/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:988898/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:781490/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1050715/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:441626/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:312442/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1794545/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1348414/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1199079/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1195136/71‑1 (MQ=255) gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1084244/71‑1 (MQ=255) | GGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824483‑3824597 |
| 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: ‑ |