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.8% | 148.6 / ‑3.8 | 45 | 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 . (0/1); new base G (25/19); total (25/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-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. |
TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC > W3110S.gb/547764‑547896 | tgcgtgGGCGCATGTCTGTTCGAAGGGTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg < 1:2058533/70‑1 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑gg > 1:329156/1‑70 (MQ=255) gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTg < 1:1303692/71‑1 (MQ=255) gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:2944908/1‑71 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt > 1:3356087/1‑70 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt > 1:2770058/1‑70 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt < 1:2741238/69‑1 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:1433054/1‑71 (MQ=255) cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCAt > 1:1532131/1‑71 (MQ=255) gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTc > 1:890993/1‑71 (MQ=255) cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCtt > 1:99020/1‑71 (MQ=255) gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCtt < 1:3245510/70‑1 (MQ=255) aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa > 1:578966/1‑58 (MQ=39) ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCa < 1:2961631/71‑1 (MQ=255) gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGa < 1:2751005/54‑1 (MQ=39) gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa < 1:782757/57‑1 (MQ=39) gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGa < 1:3859079/71‑1 (MQ=255) gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGa < 1:1205412/71‑1 (MQ=255) aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g < 1:1564423/38‑1 (MQ=255) aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAAt > 1:435818/1‑71 (MQ=255) aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGTGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt < 1:507008/71‑1 (MQ=255) aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCAtt > 1:2809300/1‑52 (MQ=39) aGATGAAGCGCAGGCGATGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTc > 1:1435108/1‑53 (MQ=38) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTcc > 1:1335456/1‑52 (MQ=39) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTcc > 1:406112/1‑52 (MQ=39) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc < 1:658914/71‑1 (MQ=255) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc > 1:2289460/1‑71 (MQ=255) tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg > 1:3100808/1‑71 (MQ=255) gAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt > 1:3447320/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt < 1:4281890/70‑1 (MQ=255) cgcAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg > 1:2298528/1‑71 (MQ=255) gcAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc > 1:4061273/1‑71 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc > 1:2988435/1‑70 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:804547/70‑1 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:975290/70‑1 (MQ=255) ggCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCt > 1:640872/1‑43 (MQ=38) cgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg < 1:4092114/71‑1 (MQ=255) cgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg < 1:3885135/71‑1 (MQ=255) gcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg > 1:1968557/1‑71 (MQ=255) gcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg > 1:1402404/1‑71 (MQ=255) tGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCtt < 1:3240763/39‑1 (MQ=37) tGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg > 1:3410575/1‑70 (MQ=255) gCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg < 1:664843/71‑1 (MQ=255) cAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt < 1:2167466/71‑1 (MQ=255) aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATACCGTCGGGCCAATGGGCGGTAt < 1:3970213/71‑1 (MQ=255) tACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt < 1:536404/71‑1 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCa > 1:4202770/1‑38 (MQ=38) gAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc > 1:4221320/1‑71 (MQ=255) | TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC > W3110S.gb/547764‑547896 |
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: ‑ |