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 | 95.0% | 126.6 / ‑1.0 | 40 | 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/2); new base G (17/21); total (17/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.99e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-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. |
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC > W3110S.gb/547761‑547893 | gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa < 1:2831291/71‑1 (MQ=255) gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg > 1:3566459/1‑71 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCTGGCGCTGGCAATACTGGAGCA‑gg‑g‑g > 1:859205/1‑70 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g‑g > 1:3875944/1‑70 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g‑g > 1:2274222/1‑70 (MQ=255) gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g‑gg > 1:2454616/1‑69 (MQ=255) gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g‑gg > 1:2449164/1‑69 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTg < 1:876204/66‑1 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACtt < 1:3281078/71‑1 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACtt > 1:2928088/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCa < 1:3660080/71‑1 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCa > 1:1316543/1‑71 (MQ=255) tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCAt > 1:465592/1‑70 (MQ=255) ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTc < 1:611123/70‑1 (MQ=255) cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCtt > 1:3080263/1‑71 (MQ=255) gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCAc < 1:1042782/71‑1 (MQ=255) tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g > 1:3035499/1‑42 (MQ=255) ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGTGCA‑gg‑g‑gg < 1:237394/44‑2 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGTG‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1185101/71‑1 (MQ=255) gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:4322588/70‑1 (MQ=255) aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAAt < 1:2106952/71‑1 (MQ=255) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc < 1:2057386/71‑1 (MQ=255) aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc > 1:3012564/1‑71 (MQ=255) gAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt > 1:1930162/1‑71 (MQ=255) gAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt > 1:2899962/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc > 1:2589897/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAAATTCATTCCTTGt > 1:442320/1‑53 (MQ=38) aGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg < 1:1205771/71‑1 (MQ=255) cAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc < 1:2820154/71‑1 (MQ=255) ggCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa < 1:1547981/71‑1 (MQ=255) cgcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg < 1:1955278/71‑1 (MQ=255) cgcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg < 1:3352112/71‑1 (MQ=255) gcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtg > 1:3189984/1‑50 (MQ=39) cTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCAcc < 1:2517227/46‑1 (MQ=39) ggCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg > 1:4204479/1‑71 (MQ=255) gCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg > 1:3402093/1‑71 (MQ=255) gCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg > 1:2875596/1‑71 (MQ=255) aTACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt < 1:4287844/71‑1 (MQ=255) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1209432/44‑1 (MQ=38) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:4419894/44‑1 (MQ=38) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:700613/53‑1 (MQ=255) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:4473427/53‑1 (MQ=255) aCTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc > 1:3917348/1‑64 (MQ=255) aCTGGAGCA‑GG‑GTGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:2974683/52‑1 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGaa > 1:1136175/1‑44 (MQ=38) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt > 1:2929024/1‑50 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc < 1:3096003/71‑1 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc < 1:2936806/71‑1 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc > 1:1198626/1‑71 (MQ=255) | GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC > W3110S.gb/547761‑547893 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |