| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,849,128 | G→T | 16.2% | Q199K (CAG→AAG) | selD ← | selenophosphate synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,849,128 | 0 | G | T | 16.2% | 107.5 / 6.4 | 37 | Q199K (CAG→AAG) | selD | selenophosphate synthase |
| Reads supporting (aligned to +/- strand): ref base G (19/12); new base T (3/3); total (22/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
CCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTC > W3110S.gb/1849061‑1849192 | ccTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTg < 1:1592135/68‑1 (MQ=255) tCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGt > 1:1105278/1‑69 (MQ=255) tGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTg < 1:2900953/69‑1 (MQ=255) ttGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGt > 1:2618228/1‑69 (MQ=255) aaaGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTc < 1:2660856/67‑1 (MQ=255) aaaGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTc < 1:2517418/67‑1 (MQ=255) aaGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAAc < 1:1313064/69‑1 (MQ=255) aaGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAAc < 1:78705/69‑1 (MQ=255) aaGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAAc < 1:860155/69‑1 (MQ=255) ggACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGttt > 1:3514622/1‑63 (MQ=255) ggACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGttt > 1:3254582/1‑63 (MQ=255) ggACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGttt > 1:2412915/1‑63 (MQ=255) ggACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGttt > 1:1526912/1‑63 (MQ=255) ggACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGttt > 1:1328424/1‑63 (MQ=255) cgcgCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGt > 1:824343/1‑52 (MQ=255) cgcgCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGa > 1:82281/1‑69 (MQ=255) gcgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGAt < 1:64021/69‑1 (MQ=255) tGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGAttttt < 1:1997683/68‑1 (MQ=255) gTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGc < 1:913490/69‑1 (MQ=255) gTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGAt > 1:3036900/1‑58 (MQ=255) gTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGAt > 1:2554187/1‑58 (MQ=255) ccGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGt < 1:2074101/69‑1 (MQ=255) gCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg < 1:2993671/69‑1 (MQ=255) gCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg < 1:2506252/69‑1 (MQ=255) tCACTTCCGTTGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg < 1:2358512/59‑1 (MQ=255) tCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTtct > 1:532349/1‑51 (MQ=255) tCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg > 1:2423463/1‑59 (MQ=255) tCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg > 1:3155656/1‑59 (MQ=255) tCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg > 1:3639543/1‑59 (MQ=255) cACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg > 1:1400779/1‑63 (MQ=255) cACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg > 1:1714896/1‑63 (MQ=255) aCTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGAttttt > 1:975151/1‑46 (MQ=255) aCTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCg > 1:2304426/1‑69 (MQ=255) tCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATc > 1:1763084/1‑69 (MQ=255) ccGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATc > 1:1399708/1‑68 (MQ=255) ccTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTc > 1:1942637/1‑68 (MQ=255) tGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg < 1:3651171/44‑1 (MQ=255) tGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTc < 1:3550028/66‑1 (MQ=255) gATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg < 1:1827541/43‑1 (MQ=255) | CCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTC > W3110S.gb/1849061‑1849192 |
| 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: ‑ |