| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,847,724 | (C)6→7 | coding (285/1044 nt) | selD ← | selenide, water dikinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,847,718 | 1 | . | C | 79.0% | 45.1 / 12.2 | 19 | coding (291/1044 nt) | selD | selenide, water dikinase |
| Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (6/9); total (7/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01 | |||||||||||
AAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAT > NC_000913/1847669‑1847762 | aaaGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTg < 1:1297972/50‑1 (MQ=255) aaGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c > 1:1794600/1‑50 (MQ=255) tttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑cccc > 1:1293192/1‑50 (MQ=255) tttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑cccc > 1:2432493/1‑50 (MQ=255) ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑ccccccca < 1:1159022/50‑3 (MQ=255) ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑ccccccca < 1:2396132/50‑3 (MQ=255) cGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATc < 1:1222431/50‑1 (MQ=38) gATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTg < 1:1317695/50‑1 (MQ=38) aTCGCAATCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGAt > 1:828892/1‑50 (MQ=255) tCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGAt < 1:1876042/50‑1 (MQ=38) tCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGAt < 1:2119578/50‑1 (MQ=38) tCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGAt < 1:1424456/50‑1 (MQ=38) cGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg > 1:139136/1‑50 (MQ=38) gCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg < 1:383465/49‑1 (MQ=38) aTCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGtt < 1:1045297/50‑1 (MQ=255) ccATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGt > 1:2166464/1‑50 (MQ=255) aTAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcg > 1:319572/1‑50 (MQ=255) tAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgc > 1:2045264/1‑50 (MQ=39) aTCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgccg > 1:26499/1‑50 (MQ=39) tCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCa < 1:1649940/50‑1 (MQ=255) ggTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCaa > 1:2606336/1‑50 (MQ=39) ggTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCa < 1:20446/49‑1 (MQ=38) ggTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCa < 1:2059737/49‑1 (MQ=38) ggTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAt < 1:702627/50‑1 (MQ=255) gTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAt < 1:1992663/50‑1 (MQ=39) | AAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAT > NC_000913/1847669‑1847762 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |