Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,645,045 | G→A | Q761* (CAA→TAA) | pbpC ← | peptidoglycan glycosyltransferase PbpC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,645,045 | 0 | G | A | 97.2% | 122.5 / ‑4.1 | 36 | Q761* (CAA→TAA) | pbpC | peptidoglycan glycosyltransferase PbpC |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (22/13); total (23/13) | |||||||||||
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.44e-01 |
TTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTT > NC_000913/2644969‑2645119 | tttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTAccc < 1:473023/80‑1 (MQ=255) ttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCa > 2:473058/1‑81 (MQ=255) taaaataaaaTGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCAc > 2:473089/1‑81 (MQ=255) aaaataaaatGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACg > 1:473075/1‑81 (MQ=255) aaaatGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTcatc > 2:473049/1‑81 (MQ=255) aaaatGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTcatc > 2:473059/1‑81 (MQ=255) tGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTcatccatc < 2:473024/81‑1 (MQ=255) gAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCa > 2:473057/1‑81 (MQ=255) gAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCa > 2:473056/1‑81 (MQ=255) gAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCa > 2:473070/1‑81 (MQ=255) ttAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAg > 1:473062/1‑80 (MQ=255) aGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCaa > 2:473083/1‑81 (MQ=255) caacaaTCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAAtt > 1:473066/1‑81 (MQ=255) aacaaTCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTg < 2:473025/81‑1 (MQ=255) aaTCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGAt > 2:473060/1‑80 (MQ=255) aTCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATaa > 1:473068/1‑81 (MQ=255) aGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATc < 2:473026/80‑1 (MQ=255) aGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCa > 2:473064/1‑81 (MQ=255) ggACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCAcc > 1:473051/1‑81 (MQ=255) gACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCt < 2:473027/81‑1 (MQ=255) cAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCttt < 2:473028/81‑1 (MQ=255) cAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCtt > 2:473073/1‑80 (MQ=255) aaaaCAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTa < 1:473029/81‑1 (MQ=255) aaCAGACTATTGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATc < 1:473030/81‑1 (MQ=255) ttGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg > 1:473076/1‑81 (MQ=255) ttGCATGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg < 1:473031/81‑1 (MQ=255) aTGACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAgg < 1:473032/81‑1 (MQ=255) gACAAATTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGt > 2:473084/1‑81 (MQ=255) aaaTTTCACTTTCTCTATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCAc > 1:473080/1‑81 (MQ=255) aTTTCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGt > 2:473078/1‑81 (MQ=255) tCACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTgc > 1:473072/1‑81 (MQ=255) cACTGTCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTgcg > 2:473065/1‑81 (MQ=255) tCGCGATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCcgc < 2:473033/81‑1 (MQ=255) gcgATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCATATTCATGTTCACGTTGCGCCCGCGt > 1:473071/1‑81 (MQ=255) gcgATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt < 1:473034/81‑1 (MQ=255) cgATTTACCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGtt < 1:473035/81‑1 (MQ=255) | TTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTT > NC_000913/2644969‑2645119 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |