Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,304,678 | C→A | 23.5% | G141C (GGC→TGC) | recD ← | multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,304,678 | 0 | C | A | 23.5% | 36.5 / 3.6 | 17 | G141C (GGC→TGC) | recD | multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease |
Reads supporting (aligned to +/- strand): ref base C (8/5); new base A (2/2); total (10/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01 |
TTCGGCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACA > NC_002947/5304565‑5304812 | ttCGGCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCAAGGCGTATCACCTGCGGccacc > 1:340410/1‑136 (MQ=255) ttCGGCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCAAGGCGTATCACCTGCGGccacc > 1:495595/1‑136 (MQ=255) gCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCAAGGCGTATCACCTGCGGCCACCAGCg < 2:495595/136‑1 (MQ=255) gCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCAAGGCGTATCACCTGCGGCCACCAGCg < 2:340410/136‑1 (MQ=255) ggCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGc < 2:330073/136‑1 (MQ=255) gATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCAttt > 2:284408/1‑136 (MQ=255) tccgcTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTAGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCAGGCGCTGATGCGCTGAAGCCACGCATTTAAGTc > 1:19622/2‑136 (MQ=255) gccgcTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTc > 1:19626/1‑136 (MQ=255) ccgcTGGAAACTCC‑TTAACGGCGGATGGAAATGCGGCGCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCa < 1:71030/136‑1 (MQ=255) cAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGc > 2:478218/1‑136 (MQ=255) c‑gcTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAg < 2:308897/136‑1 (MQ=255) ccACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGa > 1:99401/1‑136 (MQ=255) ccACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGa > 1:368738/1‑136 (MQ=255) gCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCAcc < 2:372952/136‑1 (MQ=255) gcgcTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCcgg > 1:323322/1‑136 (MQ=255) cTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAAc < 2:328467/136‑1 (MQ=255) ggCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTNCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCCNGCGGCAACGa > 2:87608/1‑136 (MQ=255) ccAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACa < 2:458809/136‑1 (MQ=255) | TTCGGCCTGACTCAGGCGCTGGCGCAGGGCCTGGTCGATGCGCCGCTCGTAACTCCAGTAACGGCGCAGGTAAAGGC‑GCTCGCCACTGCGTACAAGCGGCCGCGCTTGCTGGCCAGGCGTATCACCTGCGGCCACCAGCGGGCTGGCGCTGATGCGCTGAAGCCAGGCATTTAAGTCAAGGTTGGCCAGCAGTTGCGATGGCAGTAGCAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACA > NC_002947/5304565‑5304812 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |