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 | 95.7% | 66.6 / ‑4.8 | 23 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
Reads supporting (aligned to +/- strand): ref base . (0/1); major base G (8/13); minor base T (0/1); total (8/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
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 base substitution creates a homopolymer stretch. |
GCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824467‑3824597 | gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggg > 1:1048568/1‑71 (MQ=255) cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg < 1:2834136/71‑2 (MQ=255) aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggga < 1:465506/71‑3 (MQ=255) ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa < 1:1731962/71‑4 (MQ=255) gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc < 1:2690513/71‑1 (MQ=255) aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc > 1:812206/1‑66 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg > 1:3065803/1‑70 (MQ=255) gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggg < 1:1418018/51‑1 (MQ=255) aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg < 1:2648002/71‑1 (MQ=255) gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGATGGGGAATCGAATCCATTGTAGCGACGc < 1:3021076/71‑1 (MQ=255) ttCACCCATGAAGAGCTACTGACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg > 1:3299930/1‑69 (MQ=255) tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg < 1:2030430/46‑2 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg > 1:43609/1‑70 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg < 1:3361837/70‑1 (MQ=255) cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc < 1:1495939/71‑1 (MQ=255) cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g > 1:580136/1‑37 (MQ=38) cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGGgaa < 1:1805752/42‑1 (MQ=38) cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt > 1:3276853/1‑44 (MQ=25) tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc < 1:981878/46‑1 (MQ=255) cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg < 1:1924751/43‑1 (MQ=37) cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt < 1:3259505/71‑1 (MQ=255) gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg > 1:2641128/1‑35 (MQ=25) gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg > 1:2595243/1‑35 (MQ=25) ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg < 1:3279515/45‑1 (MQ=255) ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt < 1:177848/71‑1 (MQ=255) ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttt < 1:248638/70‑1 (MQ=255) ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa < 1:3617448/70‑1 (MQ=255) ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa < 1:915051/70‑1 (MQ=255) gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt > 1:80549/1‑71 (MQ=255) gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:1218701/71‑1 (MQ=255) | GCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824467‑3824597 |
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: ‑ |