| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,304,643 | T→C | 11.8% | E886G (GAG→GGG) | recC ← | multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,304,643 | 0 | T | C | 11.8% | 85.4 / 3.6 | 34 | E886G (GAG→GGG) | recC | multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit |
| Reads supporting (aligned to +/- strand): ref base T (12/18); new base C (2/2); total (14/20) | |||||||||||
| 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.56e-01 | |||||||||||
TAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGGC > NC_002947/5304504‑5304786 | tAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGAGCCCTTCGTcctc < 2:206222/150‑3 (MQ=255) gCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgggg > 1:41772/1‑150 (MQ=255) agACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCt > 2:166729/1‑150 (MQ=255) cAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgc < 1:166729/150‑1 (MQ=255) ggCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc < 2:356483/150‑1 (MQ=255) ggCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc < 1:516399/150‑1 (MQ=255) tCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTacac < 2:175075/150‑1 (MQ=255) cAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc > 2:334879/1‑99 (MQ=255) cAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc < 1:334879/99‑1 (MQ=255) aaGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGCATCGCGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg > 2:244488/1‑150 (MQ=255) cGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCa > 2:94355/1‑105 (MQ=255) cGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCa < 1:94355/105‑1 (MQ=255) tCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgcgc < 1:761008/150‑1 (MQ=255) ccAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaaga > 2:729708/1‑150 (MQ=255) aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaag < 2:159447/150‑1 (MQ=255) cAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCt < 1:510447/150‑1 (MQ=255) gCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgc > 2:350007/1‑129 (MQ=255) gCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgc < 1:350007/129‑1 (MQ=255) ggCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCa < 2:41772/150‑1 (MQ=255) tttCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGAc > 1:396655/1‑150 (MQ=255) aCTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAg < 2:724287/150‑1 (MQ=255) gTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGc < 2:785497/150‑1 (MQ=255) tAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa > 2:357858/1‑135 (MQ=255) tAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa < 1:357858/135‑1 (MQ=255) gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGTAGACTCAGGGCTTCAGCGTCCTGATAAg > 1:563728/1‑150 (MQ=255) ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg < 2:443527/76‑1 (MQ=255) ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg > 1:443527/1‑76 (MQ=255) tATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGTAGACTCAGGGCTTCAGCGTCCTGATAAggcg < 2:563728/150‑1 (MQ=255) tCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAggcggc < 2:396655/150‑1 (MQ=255) cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggc > 1:143110/4‑69 (MQ=25) cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggc < 2:143110/66‑1 (MQ=25) aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAg > 1:407203/1‑150 (MQ=255) ggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTTCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGc > 2:328679/1‑150 (MQ=255) aaGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGgc < 1:244488/150‑1 (MQ=255) | TAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGGC > NC_002947/5304504‑5304786 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |