Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 547,694 | A→G | 100% | pseudogene (114/261 nt) pseudogene (114/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,694 | 0 | A | G | 98.1% | 170.9 / ‑4.2 | 52 | pseudogene (114/261 nt) pseudogene (114/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 A (1/0); new base G (19/32); total (20/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-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. |
CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG > W3110S.gb/547629‑547761 | cgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa < 1:2362130/70‑1 (MQ=255) gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa < 1:3239086/70‑1 (MQ=255) cgcgCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa > 1:1362769/1‑70 (MQ=255) gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca < 1:3736428/69‑1 (MQ=255) cTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACt < 1:2458819/70‑1 (MQ=255) cTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg < 1:914918/71‑1 (MQ=255) gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa < 1:3127839/61‑1 (MQ=255) gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa < 1:3907292/61‑1 (MQ=255) tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa < 1:1391026/60‑1 (MQ=255) cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc > 1:3818353/1‑71 (MQ=255) cacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa < 1:981179/58‑1 (MQ=255) acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc < 1:2479374/71‑1 (MQ=255) ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg > 1:1998967/1‑71 (MQ=255) ggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca > 1:2205551/1‑56 (MQ=255) gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg < 1:3743065/69‑1 (MQ=255) aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt > 1:2569607/1‑55 (MQ=255) gtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACTCTGCTTCACGCCGGGCCGCCAAt > 1:2190677/1‑71 (MQ=255) gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg > 1:3853805/1‑70 (MQ=255) cAACCGGCTTCTTAACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg > 1:4148326/1‑71 (MQ=255) cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc > 1:313993/1‑60 (MQ=255) cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa > 1:4036482/1‑62 (MQ=255) ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg > 1:2699560/1‑70 (MQ=255) ggCTTCTTCACTTATCAGCGAACTAAACG‑AGGGAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg > 1:308926/1‑69 (MQ=255) gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa > 1:3608201/1‑47 (MQ=255) cttcttGACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg < 1:1553518/51‑1 (MQ=255) cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg < 1:3917132/65‑1 (MQ=255) cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg < 1:363591/69‑1 (MQ=255) cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga < 1:1183579/70‑1 (MQ=255) ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg > 1:970001/1‑47 (MQ=255) tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg > 1:1934448/1‑62 (MQ=255) tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt > 1:3340120/1‑71 (MQ=255) cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag < 1:1908306/68‑1 (MQ=255) cTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt < 1:320541/36‑1 (MQ=255) cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga < 1:2507883/65‑1 (MQ=255) cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc < 1:2822802/71‑1 (MQ=255) ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt < 1:216826/57‑1 (MQ=255) aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg < 1:1158561/70‑1 (MQ=255) cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa > 1:2730726/1‑71 (MQ=255) cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa > 1:2068920/1‑71 (MQ=255) aaCTAAACCAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa < 1:1920895/60‑1 (MQ=255) aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa < 1:2154985/71‑1 (MQ=255) aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa < 1:2203356/71‑1 (MQ=255) cTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc < 1:1305048/38‑1 (MQ=255) cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa < 1:1322879/71‑1 (MQ=255) cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa < 1:58718/71‑1 (MQ=255) cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa > 1:406718/1‑70 (MQ=255) aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAcc > 1:3381617/1‑58 (MQ=255) aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGc < 1:456805/46‑1 (MQ=255) aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAccc < 1:2226244/63‑1 (MQ=255) aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg < 1:3240741/70‑1 (MQ=255) aCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg < 1:4240133/70‑1 (MQ=255) cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg < 1:1979794/71‑1 (MQ=255) | CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG > W3110S.gb/547629‑547761 |
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: ‑ |