| 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% | 98.6 / ‑3.6 | 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 (17/13); total (17/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.52e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-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. | |||||||||||
CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC > W3110S.gb/547762‑547893 | cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g > 1:2338849/1‑71 (MQ=255) cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑ggg < 1:1924865/71‑2 (MQ=255) tgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑ggga < 1:246173/70‑3 (MQ=255) aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACt > 1:880913/1‑71 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTg < 1:266829/66‑1 (MQ=255) tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt < 1:2196796/71‑1 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa > 1:316441/1‑71 (MQ=255) cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt < 1:747939/68‑1 (MQ=255) ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTcc > 1:2182119/1‑71 (MQ=255) gAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCAtt > 1:620191/1‑48 (MQ=38) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCa > 1:96179/1‑55 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGaa > 1:393188/1‑64 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc > 1:1798187/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc > 1:429647/1‑71 (MQ=255) aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc > 1:1084545/1‑71 (MQ=255) gcgcAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg < 1:713784/71‑1 (MQ=255) aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCa < 1:600239/52‑1 (MQ=39) gcgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGACGTCGGGCCAAt < 1:2257042/71‑1 (MQ=255) cgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg > 1:708241/1‑71 (MQ=255) gcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCAtt < 1:1180278/37‑1 (MQ=37) gcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCAtt < 1:813278/37‑1 (MQ=37) cTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg < 1:2018193/71‑1 (MQ=255) cTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg < 1:1384785/71‑1 (MQ=255) aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGc < 1:793775/50‑1 (MQ=255) aTACTGGAGCA‑GGCGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt < 1:2041228/71‑1 (MQ=255) tACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt > 1:1237958/1‑70 (MQ=255) aCTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg > 1:2535928/1‑65 (MQ=255) aCTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt < 1:2202539/70‑1 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt > 1:1295271/1‑50 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa > 1:759834/1‑58 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc > 1:2204247/1‑63 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc < 1:520887/71‑1 (MQ=255) cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc > 1:2055774/1‑71 (MQ=255) gAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa > 1:31521/1‑55 (MQ=255) | CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC > W3110S.gb/547762‑547893 |
| 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: ‑ |