Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 839,209 | C→T | 8.3% | L164L (CTG→CTA) | prs ← | ribose‑phosphate pyrophosphokinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 839,209 | 0 | C | T | 8.3% | 133.0 / 2.8 | 48 | L164L (CTG→CTA) | prs | ribose‑phosphate pyrophosphokinase |
Reads supporting (aligned to +/- strand): ref base C (22/22); new base T (2/2); total (24/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01 |
GCGCCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATG > NC_002947/839062‑839356 | gcgcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCCCGGCGTGCTGGAGAGGGCGTGCGCGCACGACACCACCGATGTCCGGGGAGTCGATCATCAg < 2:371008/150‑1 (MQ=255) gcgcCCTTCCACGTCGCCTATGATGTGCATAACCTCTGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAg < 2:751927/150‑1 (MQ=255) cgcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGGCGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAgg > 1:412673/1‑150 (MQ=255) gcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGt > 1:424610/1‑150 (MQ=255) cccTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGtt > 1:138570/1‑150 (MQ=255) ttCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc < 2:38134/150‑1 (MQ=255) ttCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc < 2:150650/150‑1 (MQ=255) cGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGt > 2:601554/1‑150 (MQ=255) gtgCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGAt < 2:410416/150‑1 (MQ=255) gCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGc < 2:708574/137‑1 (MQ=255) gCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGc > 1:708574/1‑137 (MQ=255) cATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATgtc > 1:482060/1‑150 (MQ=255) gTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACggg > 1:99369/1‑150 (MQ=255) ttAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGc > 2:565343/1‑150 (MQ=255) aGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGa < 2:138570/150‑1 (MQ=255) cgcgGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGATAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGCTCCTCAGGTTCTCGACACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGa < 2:228492/150‑1 (MQ=255) cGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACgg < 1:628234/150‑1 (MQ=255) gatgatACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGGCTTCGATGTCGTCAACCAGGACTGGCTATCCGTAGTTGTTGTCGACggg > 2:720415/1‑150 (MQ=255) ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACCCCACCAATTTCGGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCga > 2:292974/1‑150 (MQ=255) ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc < 1:555240/82‑1 (MQ=255) ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc > 2:555240/1‑82 (MQ=255) ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCga > 1:33698/1‑150 (MQ=255) ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATACCga > 2:377314/1‑150 (MQ=255) cAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaa < 2:559470/150‑1 (MQ=255) cAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATACCgaa < 1:377314/150‑1 (MQ=255) cccAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGa < 1:367704/150‑1 (MQ=255) ccAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGaa > 1:307159/1‑150 (MQ=255) cAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAAt > 2:259344/1‑150 (MQ=255) ttGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaagaa > 1:156061/1‑134 (MQ=255) ttGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaagaa < 2:156061/134‑1 (MQ=255) ggCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGc > 1:423608/1‑150 (MQ=255) ggCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGc < 2:482060/150‑1 (MQ=255) gACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTg < 1:328659/150‑1 (MQ=255) aGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAg < 1:720415/150‑1 (MQ=255) cgcgTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGt < 1:250339/150‑1 (MQ=255) gcgcACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGt < 1:601554/150‑1 (MQ=255) cgcACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTg < 2:715201/150‑1 (MQ=255) accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg > 2:331572/1‑85 (MQ=255) accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg < 1:331572/85‑1 (MQ=255) accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg < 1:707464/85‑1 (MQ=255) accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg > 2:707464/1‑85 (MQ=255) ccaccGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACgg > 2:343637/1‑150 (MQ=255) accGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTc > 1:534941/1‑150 (MQ=255) accGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTc < 2:307159/150‑1 (MQ=255) cGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGGGGTTGTCGACGGCGTTAACGAAGAATCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTAc > 1:616964/1‑150 (MQ=255) ggggAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCa < 1:175816/118‑1 (MQ=255) ggggAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCa > 2:175816/1‑118 (MQ=255) gagaCGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAAAAGAGTTTAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACa > 2:477128/1‑150 (MQ=255) atcatcAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGc < 2:424610/150‑1 (MQ=255) tcatcaGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCa > 1:207814/1‑150 (MQ=255) atcaGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGGAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATg > 1:617834/1‑150 (MQ=255) | GCGCCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATG > NC_002947/839062‑839356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |