Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,646,583 | (G)6→7 | coding (743/2313 nt) | pbpC ← | penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,646,577 | 1 | . | G | 90.0% | 64.8 / 1.3 | 20 | coding (749/2313 nt) | pbpC | penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (7/11); total (8/12) | |||||||||||
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.05e-01 |
GCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGC > NC_000913/2646440‑2646683 | gCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑g < 1:188812/139‑1 (MQ=255) ttGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATggg < 2:677205/139‑1 (MQ=255) gcgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt < 1:32153/139‑1 (MQ=255) gcgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt < 1:363443/139‑1 (MQ=255) tcttcCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCtt < 2:375419/139‑1 (MQ=255) tcgtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc < 2:86792/139‑1 (MQ=255) gtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttc > 1:432415/1‑119 (MQ=255) gtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttc < 2:432415/119‑1 (MQ=255) cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc > 1:517006/1‑128 (MQ=255) cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc < 2:517006/128‑1 (MQ=255) gCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAc < 1:710136/139‑1 (MQ=255) tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcac < 2:237791/139‑1 (MQ=255) aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt > 1:627324/1‑100 (MQ=255) aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt < 2:627324/100‑1 (MQ=255) tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg > 1:687550/1‑121 (MQ=255) tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg < 2:687550/121‑1 (MQ=255) gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCg > 2:195727/1‑139 (MQ=255) tttgcttttACCGAGCATCATGCGCGAACACACCAGTTCCCGTCGCGGCCTTTGTCGGGGGGGCCAGCCAGATGGGCTGCTCCCTTGACTCTTTTACCTGCTCACCGCACCACACACCTCGCACTGCCACCCGTCCGAg > 2:385551/1‑139 (MQ=255) catcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttct > 1:77380/1‑65 (MQ=255) catcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttct < 2:77380/65‑1 (MQ=255) tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCtt > 1:106915/1‑139 (MQ=255) cgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGc > 1:404946/1‑139 (MQ=255) | GCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGC > NC_000913/2646440‑2646683 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |