| 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 | 96.8% | 103.8 / ‑3.5 | 31 | 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 (18/12); total (18/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-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. | |||||||||||
CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT > W3110S.gb/547766‑547890 | cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gggg > 1:717116/1‑70 (MQ=255)cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gggg > 1:3229554/1‑70 (MQ=255)cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gggg > 1:2966940/1‑70 (MQ=255)cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gggg > 1:2270570/1‑70 (MQ=255) gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGaa < 1:3984549/71‑1 (MQ=255) gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGaa < 1:1897560/71‑1 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACtt > 1:629580/1‑70 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACtt > 1:2687139/1‑70 (MQ=255) gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g > 1:2089637/1‑56 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCa > 1:1781601/1‑71 (MQ=255) gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTg > 1:2779886/1‑71 (MQ=255) aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTcc < 1:2646890/66‑1 (MQ=255) ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCAc > 1:1057739/1‑70 (MQ=255) ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCAc > 1:1071362/1‑70 (MQ=255) ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAt < 1:3974069/70‑1 (MQ=255) gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAtgt > 1:2004404/1‑71 (MQ=255) gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAtgtg > 1:4173592/1‑70 (MQ=255) gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAtgtg > 1:2708050/1‑70 (MQ=255) aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCa < 1:2794885/35‑1 (MQ=255) aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGaa > 1:768852/1‑71 (MQ=255) aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1532230/67‑1 (MQ=255) aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATg > 1:4443881/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTcc < 1:1014233/49‑1 (MQ=38) cAGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:861291/71‑1 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc > 1:2244085/1‑70 (MQ=255) ggCGCTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCAt > 1:4028086/1‑39 (MQ=37) gcTGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg > 1:2376156/1‑71 (MQ=255) tGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCAtgt < 1:1756702/49‑1 (MQ=39) tGGCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg < 1:242483/70‑1 (MQ=255) gCAATACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCAcc > 1:3421669/1‑43 (MQ=38) aTACTGGAGCA‑GGGGAAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:2734857/59‑1 (MQ=255) tACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:1617434/58‑1 (MQ=255) tACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt > 1:874137/1‑70 (MQ=255) tACTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt > 1:2280008/1‑70 (MQ=255) cTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGc > 1:1878149/1‑47 (MQ=255) cTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:1301593/51‑1 (MQ=255) cTGGAGCAGGGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:2167006/51‑1 (MQ=255) | CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGGG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT > W3110S.gb/547766‑547890 |
| 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: ‑ |