Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 547,831:1 | +G | 100% | pseudogene (251/261 nt) pseudogene (251/1259 nt) |
ylbE → ylbE → |
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment ECK0512:JW0508+JW0507:b4572; hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 547,831 | 1 | . | G | 97.1% | 119.6 / ‑3.5 | 35 | pseudogene (251/261 nt) pseudogene (251/1259 nt) | ylbE ylbE | ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment ECK0512:JW0508+JW0507:b4572; hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (13/21); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.00e-01 | |||||||||||
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. |
CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT > W3110S.gb/547766‑547890 | cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑gg > 1:2256068/1‑70 (MQ=255) cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTg < 1:3662671/68‑1 (MQ=255) cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTg < 1:2224007/68‑1 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:2361196/1‑68 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGTGGAAGTGaa > 1:1856796/1‑67 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:2975863/1‑66 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:3400840/1‑66 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:1940543/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:546647/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:2017979/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:2142896/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:1054324/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:1230189/1‑71 (MQ=255) ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:2003518/1‑62 (MQ=255) ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:3177352/1‑62 (MQ=255) gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGt < 1:3440991/68‑1 (MQ=255) gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGt < 1:3018572/68‑1 (MQ=255) ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCAcc < 1:1094619/71‑1 (MQ=255) gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCa < 1:2862061/70‑1 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:2253861/71‑1 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:3288975/71‑1 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1349390/71‑1 (MQ=255) aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATg < 1:1848688/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:609363/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:956895/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:2603472/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:2430982/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:2369274/71‑1 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa < 1:3664046/71‑1 (MQ=255) gcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg > 1:3498302/1‑71 (MQ=255) gCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg < 1:1982505/71‑1 (MQ=255) gCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg < 1:1804667/71‑1 (MQ=255) aaTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTa > 1:2345383/1‑71 (MQ=255) aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt < 1:3146746/71‑1 (MQ=255) aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt < 1:2679645/71‑1 (MQ=255) tACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:1465360/58‑1 (MQ=255) | CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT > W3110S.gb/547766‑547890 |
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: ‑ |