| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,646,583 | (G)6→7 | coding (743/2313 nt) | pbpC ← | peptidoglycan glycosyltransferase PbpC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,646,577 | 1 | . | G | 92.5% | 133.6 / 3.8 | 40 | coding (749/2313 nt) | pbpC | peptidoglycan glycosyltransferase PbpC |
| Reads supporting (aligned to +/- strand): ref base . (2/1); new base G (10/27); total (12/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.09e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01 | |||||||||||
CCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAG > NC_000913/2646448‑2646710 | ccTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACAGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑gggggggcca < 2:493404/140‑5 (MQ=255) tttCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGgccagc > 2:512620/1‑140 (MQ=255) cAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGa < 2:235547/140‑1 (MQ=255) gcgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttc < 2:30593/140‑1 (MQ=255) cgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGttctt > 1:267456/1‑140 (MQ=255) cAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTcc < 1:854775/140‑1 (MQ=255) aCGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtttt < 2:800909/140‑1 (MQ=255) aCGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGGCAGGCAGATGGGGTCTTCCCTTTACTCctttt > 2:201734/1‑139 (MQ=255) gtcgtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc < 1:482044/140‑1 (MQ=255) tcgtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc > 2:436073/1‑140 (MQ=255) gtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTg < 1:26830/140‑1 (MQ=255) gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg < 2:773687/140‑1 (MQ=255) gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg < 1:791832/140‑1 (MQ=255) gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg < 1:506353/140‑1 (MQ=255) cATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACca < 1:543609/140‑1 (MQ=255) tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcaca > 2:350583/1‑140 (MQ=255) cAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacaca < 1:316247/140‑1 (MQ=255) agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGc < 1:331348/140‑1 (MQ=255) agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGc > 2:839222/1‑140 (MQ=255) tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGTTCACGGGACCACACACCTTGCAc > 1:459171/1‑140 (MQ=255) gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGcc < 1:350583/140‑1 (MQ=255) aTTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCa < 2:283136/140‑1 (MQ=255) ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTg < 2:255417/131‑1 (MQ=255) ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTg > 1:255417/1‑131 (MQ=255) tttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc < 1:436073/140‑1 (MQ=255) tttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATc < 1:170474/140‑1 (MQ=255) gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt < 2:593418/126‑1 (MQ=255) gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt > 1:593418/1‑126 (MQ=255) gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGt > 2:55367/1‑140 (MQ=255) cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc < 2:305513/140‑1 (MQ=255) cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGa < 1:456733/140‑1 (MQ=255) cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGa < 1:2770/140‑1 (MQ=255) cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGa < 2:313657/140‑1 (MQ=255) cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAg < 2:267456/140‑1 (MQ=255) tttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGc < 1:915234/140‑1 (MQ=255) catGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCt < 2:160550/140‑1 (MQ=255) atGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCtt > 1:63377/1‑140 (MQ=255) aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGTACTTTATTACGCGCGGCTTCGGCACGCt < 1:185829/140‑1 (MQ=255) aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCt < 2:267004/140‑1 (MQ=255) aaCAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCt < 1:512620/140‑1 (MQ=255) cGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCggg > 1:855482/1‑140 (MQ=255) tttGTCGGGGGGGCCAGCCAGATGGGCTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCGCTTTATTACGCGCGGCTTCGGCACGCTTCGTCCAACGATCCGGGTGAAg < 2:100785/140‑1 (MQ=255) | CCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAG > NC_000913/2646448‑2646710 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |