| 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 | 98.3% | 190.7 / ‑3.9 | 59 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (24/34); total (24/35) | |||||||||||
| 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.63e-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 indel expands or contracts a homopolymer stretch. | |||||||||||
GCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑G‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824467‑3824597 | gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑ggg > 1:2090835/1‑71 (MQ=255) aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑g < 1:320288/67‑1 (MQ=255) cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcgaa < 1:3976126/71‑1 (MQ=255) gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcgaatc < 1:1468296/71‑1 (MQ=255) aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCa < 1:1019521/71‑1 (MQ=255) aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:1823944/54‑1 (MQ=255) aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTg < 1:2992418/71‑1 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGt < 1:459663/71‑1 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGt < 1:952348/71‑1 (MQ=255) gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGa > 1:3968067/1‑71 (MQ=255) gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGa > 1:3596064/1‑71 (MQ=255) gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATcc > 1:4191603/1‑58 (MQ=255) gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACg > 1:741107/1‑70 (MQ=255) ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCa < 1:1419014/71‑1 (MQ=255) ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCa < 1:4434611/71‑1 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGc > 1:393889/1‑63 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAg < 1:2680082/70‑1 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑gg > 1:15250/1‑43 (MQ=38) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑gg > 1:2206608/1‑43 (MQ=38) cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAg < 1:480963/71‑1 (MQ=255) ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGaa > 1:3925294/1‑69 (MQ=255) ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAgg < 1:4197056/71‑1 (MQ=255) cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggc < 1:2837832/71‑1 (MQ=255) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcgg < 1:15012/71‑1 (MQ=255) gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcggc > 1:471214/1‑71 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg > 1:2356664/1‑71 (MQ=255) gagCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt < 1:4356356/71‑1 (MQ=255) agCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcgaat > 1:1979395/1‑42 (MQ=25) agCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcgaat > 1:1870504/1‑42 (MQ=25) agCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTg < 1:3520350/48‑1 (MQ=255) gCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg > 1:871215/1‑71 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCa < 1:1349279/56‑1 (MQ=255) cTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCa < 1:3813840/56‑1 (MQ=255) cTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa > 1:3997264/1‑71 (MQ=255) cTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa < 1:3481175/71‑1 (MQ=255) tCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCAt > 1:864592/1‑40 (MQ=25) tCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGa > 1:4263843/1‑48 (MQ=255) caccTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt > 1:458988/4‑71 (MQ=255) accTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt > 1:4326285/3‑70 (MQ=255) accTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt > 1:125311/3‑70 (MQ=255) accTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCAATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt > 1:3976496/3‑70 (MQ=255) ccTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTg < 1:2365548/38‑1 (MQ=25) ccTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAg < 1:1558109/49‑1 (MQ=255) cTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt < 1:2176766/71‑1 (MQ=255) ttgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGaa > 1:1159887/1‑50 (MQ=255) ttgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGaa > 1:1164930/1‑50 (MQ=255) ttgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggc < 1:2430966/53‑1 (MQ=255) ttgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg < 1:1179984/68‑1 (MQ=255) ttgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt > 1:1862383/1‑71 (MQ=255) tgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcgg > 1:3640854/1‑54 (MQ=255) tgttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcgg > 1:1352416/1‑54 (MQ=255) gttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGACAAACTGAttt > 1:3947978/1‑70 (MQ=255) ttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAg < 1:2833094/45‑1 (MQ=255) ttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAg < 1:2468947/45‑1 (MQ=255) ttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt < 1:4259466/71‑1 (MQ=255) ttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt < 1:3103152/71‑1 (MQ=255) ttgGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCCGCGCTGGCAAACTGAttttt < 1:490716/71‑1 (MQ=255) ttgGCGCTGGCCCGA‑GTGGGAATCGAATCCATTGTAGCGACGCAg < 1:1692655/45‑1 (MQ=255) gcgctGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAg < 1:3066463/68‑1 (MQ=255) ctGGCCCGAGG‑GGGAATCGAATCGATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCg < 1:2511444/71‑1 (MQ=255) ctGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCg < 1:2912757/71‑1 (MQ=255) ctGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCg < 1:2285509/71‑1 (MQ=255) gCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:516990/71‑1 (MQ=255) gCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt < 1:3289068/71‑1 (MQ=255) | GCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑G‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT > W3110S.gb/3824467‑3824597 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |