Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,003,834 | G→C | A116P (GCA→CCA) | xdhC → | putative xanthine dehydrogenase iron‑sulfur‑binding subunit XdhC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,003,834 | 0 | G | C | 96.8% | 103.0 / ‑3.4 | 31 | A116P (GCA→CCA) | xdhC | putative xanthine dehydrogenase iron‑sulfur‑binding subunit XdhC |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base C (14/16); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 |
ATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGT > NC_000913/3003757‑3003908 | aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCa < 1:548258/80‑1 (MQ=255) agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGGTTATGGCTACCACGCCAATGCTgg > 2:548283/1‑81 (MQ=255) agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTg < 1:548259/80‑1 (MQ=255) agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTgg > 2:548299/1‑81 (MQ=255) agGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGa < 2:548261/81‑1 (MQ=255) agGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGa < 2:548260/81‑1 (MQ=255) tGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcg > 2:548317/1‑81 (MQ=255) gAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcgc > 2:548310/1‑79 (MQ=255) aaaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcgc > 1:548311/1‑78 (MQ=255) aaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCgaga < 2:548262/81‑1 (MQ=255) aTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGaa < 2:548263/81‑1 (MQ=255) tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAg < 1:548264/81‑1 (MQ=255) cGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGcc < 2:548265/81‑1 (MQ=255) cgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCAtt < 2:548266/81‑1 (MQ=255) cagtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCAt < 2:548267/80‑1 (MQ=255) gtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTAc < 1:548268/81‑1 (MQ=255) gcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACg > 1:548314/1‑80 (MQ=255) cGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACgg > 2:548326/1‑80 (MQ=255) gggTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGa > 1:548304/1‑80 (MQ=255) ggTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGaaa < 1:548269/81‑1 (MQ=255) ttttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAAtt < 1:548270/81‑1 (MQ=255) ttttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAAtt > 2:548322/1‑81 (MQ=255) ttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATtcg > 1:548298/1‑81 (MQ=255) tACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATtcgtc < 1:548271/80‑1 (MQ=255) aCGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgc > 1:548329/1‑81 (MQ=255) aCGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgc > 1:548293/1‑81 (MQ=255) cctgATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGc < 1:548272/81‑1 (MQ=255) ggCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCt > 1:548327/1‑81 (MQ=255) cTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCttt < 2:548273/81‑1 (MQ=255) cTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCttt > 2:548302/1‑81 (MQ=255) ccACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGt < 1:548274/80‑1 (MQ=255) | ATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGT > NC_000913/3003757‑3003908 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |