Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,399,985 | (C)7→8 | coding (252/516 nt) | ogt ← | O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,399,978 | 1 | . | C | 94.1% | 49.8 / ‑2.9 | 17 | coding (259/516 nt) | ogt | O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (7/9); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01 |
ATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑C‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT > NC_000913/1399849‑1400100 | aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑cccccccg < 2:959471/139‑3 (MQ=255) aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑cccccccg < 1:760986/139‑3 (MQ=255) aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑cccccccg < 2:321279/139‑3 (MQ=255) aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CTCCCCCCg < 1:438161/139‑1 (MQ=255) tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCa < 1:1035658/139‑1 (MQ=255) gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGa > 1:545147/1‑139 (MQ=255) aaCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGt < 2:545147/139‑1 (MQ=255) aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGGAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGGGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCaat < 2:65222/139‑1 (MQ=255) ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt < 2:990246/139‑1 (MQ=255) ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt < 1:738315/139‑1 (MQ=255) tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCaa > 1:1234766/1‑101 (MQ=255) tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCaa < 2:1234766/101‑1 (MQ=255) tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTg > 1:1166086/1‑139 (MQ=255) aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa < 2:1166086/139‑1 (MQ=255) gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑cccccccg < 2:290238/61‑3 (MQ=255) gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑cccccccg > 1:290238/1‑59 (MQ=255) gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑ccccc > 1:772789/1‑58 (MQ=255) gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c‑ccccc < 2:772789/58‑1 (MQ=255) ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTggg < 1:671326/69‑1 (MQ=255) ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTggg > 2:671326/1‑69 (MQ=255) ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg > 1:66162/1‑139 (MQ=255) tttCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCaat > 1:857147/1‑58 (MQ=255) tttCCAGACTTCGCGCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCaat < 2:857147/58‑1 (MQ=255) cgcgCTGAAATGGCGTCC‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt > 2:953315/1‑139 (MQ=255) | ATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑C‑CCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT > NC_000913/1399849‑1400100 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |