Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,195,945 | 0 | A | C | 71.0% | 55.5 / 24.4 | 31 | I275L (ATC→CTC) | icd | isocitrate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (6/3); new base C (12/10); total (18/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGTGAAGAGTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTAT > NC_000913/1195868‑1196018 | cgcgTGAAGAGTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCg > 1:219559/1‑81 (MQ=255) aaGAGTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCAtt < 1:219528/80‑1 (MQ=255) gagTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTaaa < 2:219530/81‑1 (MQ=255) gagTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTa > 1:219557/1‑79 (MQ=255) gagTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTaaa < 2:219529/81‑1 (MQ=255) gTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGa < 2:219531/81‑1 (MQ=255) ggCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGt < 2:219532/80‑1 (MQ=255) ggCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGt > 2:219579/1‑80 (MQ=255) gAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCt > 2:219574/1‑81 (MQ=255) aTCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATg < 2:219533/79‑1 (MQ=255) gACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCAtt < 1:219534/80‑1 (MQ=255) gACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTc > 1:219572/1‑81 (MQ=255) gACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTc > 1:219569/1‑81 (MQ=255) gACGGTGGCCCCTTGCTGAAAGTTAAAAACCCGAACACTCGCCAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTc > 2:219565/1‑81 (MQ=255) cccGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCaa > 2:411482/1‑79 (MQ=255) gTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCaa > 2:219568/1‑76 (MQ=255) gAAAGTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCctgct > 2:219580/1‑81 (MQ=255) aaaGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCctgc < 2:219535/79‑1 (MQ=255) gTTAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGt < 2:219536/81‑1 (MQ=255) ttAAAAACCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTc > 1:219566/1‑81 (MQ=255) aaaaaCCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTcc < 1:219537/80‑1 (MQ=255) aaaaaCCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTcc < 2:219538/80‑1 (MQ=255) aaaaCCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCgg < 2:219539/81‑1 (MQ=255) aCCCGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCt > 2:219564/1‑80 (MQ=255) cccGAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTg > 1:219575/1‑80 (MQ=255) gAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAata > 2:219563/1‑81 (MQ=255) gAACACTGGCAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAata > 2:219570/1‑81 (MQ=255) ggCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGt > 1:219567/1‑80 (MQ=255) cAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTAt > 1:219578/1‑81 (MQ=255) cAAAGAGCTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTAt < 1:219540/81‑1 (MQ=255) cAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTAt > 1:219583/1‑81 (MQ=255) | CGCGTGAAGAGTTTGGCGGTGAACTGATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTAT > NC_000913/1195868‑1196018 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |