Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,212,539 | T→G | 14.6% | intergenic (+55/+39) | PP_4593 → / ← PP_4594 | hypothetical protein/cystathionine gamma‑synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,212,539 | 0 | T | G | 14.6% | 86.6 / 7.4 | 41 | intergenic (+55/+39) | PP_4593/PP_4594 | hypothetical protein/cystathionine gamma‑synthase |
Reads supporting (aligned to +/- strand): ref base T (16/19); new base G (2/4); total (18/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.79e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01 |
CGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACC > NC_002947/5212391‑5212686 | cGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACcgtc < 1:176050/150‑1 (MQ=255) aCCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg < 1:312326/150‑1 (MQ=255) cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTTACTCCACCTGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCg < 1:139055/150‑1 (MQ=255) acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGT‑GTAGGAGCggg < 2:686999/150‑1 (MQ=255) ccGTCAACCGCCTGCGCAAAGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCgggcttgcc < 1:496974/150‑7 (MQ=37) gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCcg > 2:881399/1‑150 (MQ=255) gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCcg > 2:881075/1‑150 (MQ=255) ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGaaa > 1:863008/1‑150 (MQ=255) gcgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGCTTGCCCGCGaagaggccg < 1:855940/150‑8 (MQ=25) cagagcagAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGATAAAAC‑GCCGTCTGAGGCGCGGGTTTACCCGCGAAAGGGCCt > 2:820412/8‑150 (MQ=255) ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc < 1:829264/150‑1 (MQ=255) tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑CCCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca < 2:993477/150‑1 (MQ=255) tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca < 1:881399/150‑1 (MQ=255) tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca < 1:881075/150‑1 (MQ=255) gACAAGATCGAACGGATCATTCATCCGCCCAAATCGTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcgc > 1:366364/1‑150 (MQ=255) aCAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt > 2:962402/1‑150 (MQ=255) cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCtt < 2:464725/150‑1 (MQ=255) aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCATCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGctt > 1:779961/1‑106 (MQ=25) aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCATCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGctt < 2:779961/106‑1 (MQ=25) gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg < 2:535561/142‑1 (MQ=255) gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg > 1:535561/1‑142 (MQ=255) aCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTCTTGACACCTACGGCGTCTTCGCTGGCAAGCCCGAAACAAC‑TCCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGccc < 1:520085/150‑1 (MQ=37) tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGAGCTTTCAGCGCCACTTcc > 2:793951/1‑150 (MQ=255) ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAgg > 2:156973/1‑150 (MQ=255) ccGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc < 1:156973/150‑1 (MQ=255) aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg < 2:863008/150‑1 (MQ=255) cTGACAGGCATCTATTGAATCCACCGGCCTCTTCCAGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCAc < 2:366364/150‑1 (MQ=255) gACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACTGGCGAAAGGGCCGGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg > 1:411485/1‑150 (MQ=255) gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg > 1:936034/1‑79 (MQ=255) gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg < 2:936034/79‑1 (MQ=255) tGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATcc < 1:795017/150‑1 (MQ=255) tGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATcc < 1:793951/150‑1 (MQ=255) cACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTGTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCg > 2:739139/1‑150 (MQ=255) gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATGCTCTGTCCCTACCGACAACCGCAGCAGGTTCTCGc > 2:88328/1‑150 (MQ=255) ccGCTCCTACGACGGT‑GTAGGAGCGGGCTTACCCGTGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc > 2:1075478/1‑69 (MQ=17) ccGCTCCTACGACGGT‑GTAGGAGCGGGCTTACCCGTGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc < 1:1075478/69‑1 (MQ=17) ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg < 2:437359/107‑1 (MQ=38) ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg > 1:437359/1‑107 (MQ=38) ccGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCGCTAACCCAACCGACAACCGCAGCAGCTTCCCGCTg > 2:456427/1‑150 (MQ=39) gCTCCTAC‑GACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcagca < 1:375864/137‑1 (MQ=38) tAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCa < 1:962402/150‑1 (MQ=255) acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCa > 1:227587/1‑150 (MQ=255) cgtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCAcc > 1:750278/1‑150 (MQ=255) gtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa > 1:760736/1‑51 (MQ=255) gtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa < 2:760736/51‑1 (MQ=255) tcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGa < 2:127748/137‑1 (MQ=255) tcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGa > 1:127748/1‑137 (MQ=255) | CGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACC > NC_002947/5212391‑5212686 |
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: ‑ |