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 | 90.9% | 27.4 / ‑2.7 | 11 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (3/7); total (4/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.64e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
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. |
GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGC > W3110S.gb/3824473‑3824572 | gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg < 1:897506/71‑1 (MQ=255) cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcgaa > 1:1556389/1‑71 (MQ=255) aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCa < 1:704075/71‑1 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g > 1:1061868/1‑41 (MQ=38) cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg < 1:103367/48‑1 (MQ=255) cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAg > 1:1211786/1‑71 (MQ=255) ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAgg < 1:1530249/71‑1 (MQ=255) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg < 1:1425942/44‑1 (MQ=25) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg < 1:908625/44‑1 (MQ=25) agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc > 1:1337977/1‑71 (MQ=255) tACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCAt < 1:1414918/43‑1 (MQ=25) tgttgGCGCTGGCCCGA‑GGGGAAATCGAATCCATTGTa > 1:19124/1‑38 (MQ=37) | GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGC > W3110S.gb/3824473‑3824572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |