Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,310,841 | C→T | 40.9% | E110K (GAA→AAA) | pnp ← | polynucleotide phosphorylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,310,841 | 0 | C | T | 40.9% | 3.5 / 14.0 | 20 | E110K (GAA→AAA) | pnp | polynucleotide phosphorylase |
Reads supporting (aligned to +/- strand): ref base C (5/6); new base T (2/6); total (8/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.33e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01 |
CTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATA > NC_000913/3310706‑3310982 | cTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTAGCGGGTTAACAGAAACCACGGTGGCGATAACCTGCACTTCGT‑TGA‑Cg > 2:153422/1‑143 (MQ=255) gcgGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGAGAAGATAGAAGTGTGACACAACTGCCCAACAGACACCACGGTGGCGAGACCCGGCACTTCGT‑GGC‑CGAAGc > 2:67443/1‑143 (MQ=255) ggCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAACACTGCGGAAGCACCAATCATCGCGACGATATCCGGGATAACTAGCGGGATAACAGAAACAACGATGGCGATAACCCGAACCTTGT‑TGA‑CGAAGCCtt > 2:200730/1‑143 (MQ=255) gcagcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTc < 1:236231/143‑1 (MQ=255) aCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGa > 1:207379/1‑143 (MQ=255) ggCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCggg < 2:204606/143‑1 (MQ=255) ggCCAATCAACCGAATACCCGACAGCGACAGCGCTGCCGAAGCACCAATCCTCGCGACGATATCCGGGTTACCTTACCTGGTTATAGAATCCACGGTGGCGATAACCTGACCTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCggg < 2:127991/143‑1 (MQ=255) ttGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAAt > 2:71916/1‑143 (MQ=255) cagacagaGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATc < 1:57445/143‑1 (MQ=255) gaccgggcTGCGGAAGCACCAGATATCGCGACGTTATCCGGGGTAACTGGCGGGTTGACAGAAACGACGGTGGCGATAAGCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGAcgcgc < 2:223637/137‑1 (MQ=255) gAAGCACCAATCAGCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTc < 1:59276/143‑1 (MQ=255) cccAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCtt < 1:179427/142‑1 (MQ=255) ccAATCAGCGCGACGAGATCTGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTc < 2:7485/143‑1 (MQ=255) cAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGCCAAAACCGACGTTATGACGTGCAATTAGTGTGGTTTTTGGCT‑‑CGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGGTTCGCCTTCg > 1:235401/1‑143 (MQ=255) aaTCATCCCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTTGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGc < 2:166683/143‑1 (MQ=255) gacgcTATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGggcggc < 1:156342/143‑1 (MQ=255) gaTATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGAGCAGGGTTTCGCCTTCGCTTGGGCGGCCtt > 2:24478/1‑143 (MQ=255) ttAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGgaa > 1:107715/1‑143 (MQ=255) aCTTTGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTGCGCGTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTg < 1:24478/143‑1 (MQ=255) tCGC‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCCATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATa < 1:134362/143‑1 (MQ=255) | CTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGT‑TGA‑CGAAGCCTTCCGGGAACAGCGGGCGAATCGGGCGGTCAATCAGACGCGCGATCAGGGTTTCGCCTTCGCTTGGGCGGCCTTCACGACGGAAGAAGCTACCCGGGATACGACCAGCAGCGTAGGTACGCTCCTGATA > NC_000913/3310706‑3310982 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |