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 | 94.5% | 111.8 / ‑0.7 | 36 | 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/1); new base G (17/17); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
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: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC > W3110S.gb/547764‑547896 | tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g < 1:1939924/71‑1 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑g‑g > 1:3028596/1‑70 (MQ=255) cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGa < 1:440141/71‑1 (MQ=255) cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGa < 1:1283608/71‑1 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCa > 1:2948435/1‑71 (MQ=255) ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GCGAAGTGAACTTCATTc < 1:1625386/70‑1 (MQ=255) cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCtt > 1:2511483/1‑71 (MQ=255) aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGt < 1:918651/71‑1 (MQ=255) gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAt > 1:348079/1‑71 (MQ=255) gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAt > 1:2513037/1‑71 (MQ=255) ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCAc > 1:3579193/1‑67 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1247154/71‑1 (MQ=255) cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:2541513/71‑1 (MQ=255) gAAAGATGAAGCGCATGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgt < 1:2842004/69‑1 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGt > 1:292177/1‑53 (MQ=39) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc < 1:3007304/71‑1 (MQ=255) gcgcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:11380/54‑1 (MQ=39) gcgcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:2548129/54‑1 (MQ=39) cgcTGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:420687/53‑1 (MQ=39) tGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:1209504/50‑1 (MQ=255) tGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:3012554/50‑1 (MQ=255) tGGCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc < 1:2171679/56‑1 (MQ=255) ggCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtg > 1:1116702/1‑47 (MQ=255) ggCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtg > 1:2425221/1‑47 (MQ=255) gCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGaa > 1:1882796/1‑50 (MQ=255) gCAATACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg > 1:2191164/1‑71 (MQ=255) aTACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa < 1:49102/61‑1 (MQ=255) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCAc < 1:51615/38‑1 (MQ=38) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:449137/44‑1 (MQ=38) tACTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCAtgtg < 1:814000/44‑1 (MQ=38) tACTGGAGCA‑GGTG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt > 1:2876487/1‑70 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCa > 1:3195245/1‑35 (MQ=37) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc > 1:1782034/1‑71 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc > 1:3682279/1‑71 (MQ=255) cTGGAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc > 1:2490106/1‑71 (MQ=255) ggAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt > 1:459259/1‑70 (MQ=255) gAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc > 1:2946143/1‑71 (MQ=255) gAGCAGGG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc > 1:2548185/1‑71 (MQ=255) | TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑G‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC > W3110S.gb/547764‑547896 |
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: ‑ |