Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,645,013 | C→T | *771* (TAG→TAA) | 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,645,013 | 0 | C | T | 97.6% | 123.3 / ‑7.2 | 41 | *771* (TAG→TAA) | pbpC | penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (22/18); total (23/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC > NC_000913/2644884‑2645148 | tCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGATGACTGTTGGTGAAGGTAACATTTTGTTACAATAAAATGAGATTTAGCAACGATCAAGGACAAAACAGACtctgtctca > 1:130003/1‑131 (MQ=255) tgttgttgtCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGAc < 2:1850292/139‑1 (MQ=255) tGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACt > 2:458434/1‑139 (MQ=255) cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt > 2:1913715/1‑139 (MQ=255) cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt > 1:720288/1‑139 (MQ=255) tACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAttt > 2:498542/1‑139 (MQ=255) cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa < 2:955505/113‑1 (MQ=255) cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa > 1:955505/1‑113 (MQ=255) cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGc > 1:1101432/1‑139 (MQ=255) agaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc < 2:1101432/139‑1 (MQ=255) ggCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg < 1:1913715/139‑1 (MQ=255) taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCAc < 2:720288/139‑1 (MQ=255) gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCaa > 1:236342/1‑139 (MQ=255) ccGTTGACTATTGATGAAGGTAACATTTTTTT‑AAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCa > 1:908565/1‑139 (MQ=255) atgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTc < 1:458434/139‑1 (MQ=255) ggTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg < 1:498542/139‑1 (MQ=255) cATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAggg < 2:236342/139‑1 (MQ=255) ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc < 1:1397136/80‑1 (MQ=255) ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc > 2:1397136/1‑80 (MQ=255) ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACg > 2:1957808/1‑139 (MQ=255) tttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg < 1:1467049/128‑1 (MQ=255) tttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg > 2:1467049/1‑128 (MQ=255) aaataaaatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTgcg < 1:318293/139‑1 (MQ=255) aatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCcgcg < 2:908565/139‑1 (MQ=255) tttAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGtt > 2:118689/1‑139 (MQ=255) tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt < 1:622277/131‑1 (MQ=255) tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt > 2:622277/1‑131 (MQ=255) aGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAAc > 1:1917665/1‑139 (MQ=255) gCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACg < 1:118689/139‑1 (MQ=255) aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCg > 2:560228/1‑139 (MQ=255) aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCg > 2:33341/1‑139 (MQ=255) aTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGc > 1:989623/1‑139 (MQ=255) aGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCAt > 1:1017479/1‑139 (MQ=255) aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt > 2:962180/1‑114 (MQ=255) aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt < 1:962180/114‑1 (MQ=255) aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc < 2:1088012/42‑1 (MQ=255) aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc > 1:1088012/1‑42 (MQ=255) aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATaa > 2:1147640/1‑70 (MQ=255) aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATaa < 1:1147640/70‑1 (MQ=255) aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAac > 2:1754048/1‑139 (MQ=255) aGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAaccac < 1:33341/139‑1 (MQ=255) | TCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC > NC_000913/2644884‑2645148 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |