| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,212,540 | C→T | 21.1% | intergenic (+56/+38) | 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,540 | 0 | C | T | 21.1% | 32.2 / 3.9 | 20 | intergenic (+56/+38) | PP_4593/PP_4594 | hypothetical protein/cystathionine gamma‑synthase |
| Reads supporting (aligned to +/- strand): ref base C (9/6); new base T (2/2); total (11/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
GACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCC > NC_002947/5212401‑5212688 | gacgaAGACGTCAACCGCCTGCGCAATGAAAGGTTGCTGATCAAGGACAGGATCGAACGGATCATTCATCCGGCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGT‑GTAGGAGCgg < 2:205907/150‑1 (MQ=255) cGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAg < 1:52010/150‑1 (MQ=255) gATCAAGGACAAGATCGAACGGATCATTCATCAGCCCAAATACTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGCCCCTCAGACAggcc > 2:107321/1‑147 (MQ=38) cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc < 1:196353/142‑1 (MQ=255) cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc > 2:196353/1‑142 (MQ=255) gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCa > 2:178749/1‑150 (MQ=255) aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcc < 1:178749/150‑1 (MQ=255) cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGAGCTTTCAACGCCCCATCCAGGTCCTCGACCa > 1:370300/1‑150 (MQ=255) tCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTc > 2:37905/1‑150 (MQ=255) tCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc < 1:37905/150‑1 (MQ=255) aTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCCATCCCtc > 2:122368/1‑149 (MQ=255) tCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc > 1:78667/1‑150 (MQ=255) tCCACCGACCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAANGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc < 1:107321/150‑1 (MQ=255) cGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACa > 2:471577/1‑150 (MQ=255) cgcgGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcagca < 1:471577/150‑1 (MQ=255) gCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtc > 1:107145/1‑150 (MQ=255) cccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgca < 1:62636/137‑1 (MQ=39) cccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgca > 2:62636/1‑137 (MQ=255) ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGt > 2:110076/1‑103 (MQ=38) ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGt < 1:110076/103‑1 (MQ=38) cgtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtct < 1:80712/133‑1 (MQ=255) cgtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtct > 2:80712/1‑133 (MQ=255) tcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCAcccc < 1:351836/150‑1 (MQ=255) | GACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCC > NC_002947/5212401‑5212688 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |