| 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 | 88.9% | 49.8 / 0.0 | 18 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
| Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (5/11); total (6/12) | |||||||||||
| 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.46e-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. | |||||||||||
CAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC > W3110S.gb/3824468‑3824593 | cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggg‑gg < 1:2112123/71‑2 (MQ=255) gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcgaatc > 1:1300405/1‑71 (MQ=255) cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g > 1:2429069/1‑57 (MQ=255) aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGGTGAATCGAATCCATTg > 1:2041341/1‑71 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGt < 1:347950/71‑1 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGt < 1:2533899/71‑1 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAg < 1:192656/70‑1 (MQ=255) aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAggcg > 1:1658272/1‑71 (MQ=255) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:748659/35‑1 (MQ=37) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:2168904/35‑1 (MQ=37) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:2054062/35‑1 (MQ=37) tCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa < 1:1784528/68‑1 (MQ=255) accTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAAc > 1:911657/3‑69 (MQ=255) tgttgGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCa < 1:1307179/46‑1 (MQ=255) tgttgGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt > 1:2036974/1‑70 (MQ=255) tgttgGCGCTGGCCCGA‑GGGTGAATCGAATCCATTGTAGCGACGCa < 1:695269/46‑1 (MQ=255) ttgGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt < 1:2197652/71‑1 (MQ=255) ggCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa > 1:1973719/1‑45 (MQ=255) gcgctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAg < 1:1668334/68‑1 (MQ=255) gctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc < 1:1893603/70‑1 (MQ=255) gctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc < 1:832276/70‑1 (MQ=255) ctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc < 1:1366897/70‑1 (MQ=255) tGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCa < 1:911128/37‑1 (MQ=38) | CAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC > W3110S.gb/3824468‑3824593 |
| 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: ‑ |