| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 93,173 | A→C | R587S (AGA→AGC) D3A (GAT→GCT) |
ftsI → murE → |
transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3 UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 93,173 | 0 | A | C | 96.3% | 77.0 / NA | 27 | R587S (AGA→AGC) D3A (GAT→GCT) | ftsI murE | transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3 UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (9/17); total (10/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGATCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGATA > NC_000913/93030‑93317 | tactacGGCGGCGCCGTTTCCGGGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTa < 1:341199/149‑1 (MQ=255) ccGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACct < 2:60757/149‑1 (MQ=255) gTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACcttc < 1:355362/149‑1 (MQ=255) tttCCGCGCCGTTCTTGGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACcttct < 2:191238/149‑1 (MQ=255) gcTTGGGTGCCACCATGGGCGGCGTATTGCGTACCACAAACATCGAGCCGGATGCGCTGACAACGGGCAATAAACATGCCCTTATGCCTACTCAAGGCTAGGGTACAGGTGCCACCTCGTAATTTGCGCGACATTCTTGCTCCGTGGGt < 2:189188/148‑1 (MQ=255) tgttgTGCCATCATGTGCGGCGTATTGCGTACCATGAACATCGGGCCGGATGCGCTGACAACGCGGGAGAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTCTCGCGACCTTCTTGCTCCGTGGGTGc < 2:139876/145‑1 (MQ=255) gTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGATCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGa > 2:229109/1‑149 (MQ=255) ccATCAGGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGc < 2:360886/149‑1 (MQ=255) atGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGCGGGTGCCAGACGCACCtt > 2:365417/1‑149 (MQ=255) gggCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGATAGGTTGCAGCTAGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTcg < 1:441082/149‑1 (MQ=255) tgtGCGTGCCGTGAACAGCGGGGCGGATGCGCTGACAGCGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTg < 2:109042/147‑1 (MQ=255) gTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCgaga < 2:127693/149‑1 (MQ=255) aTGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGa < 1:38206/149‑1 (MQ=255) aTGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGa > 2:435637/1‑149 (MQ=255) gATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGAGGACCCTCGACAGCCgtg > 2:236135/1‑149 (MQ=255) ggCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCggcgggcg > 2:208010/1‑149 (MQ=255) gCGATAAAAATGAATTTGGGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGa < 2:125996/149‑1 (MQ=255) aacaTAAAAATGAAGTTGTGATTAGTCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGa < 2:188902/146‑1 (MQ=255) gaaaaTGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGCACTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGGGCGATGACACTCGACAGCCGTCTGGCTGCGGCCGGCGAtctcg > 2:434193/2‑148 (MQ=255) aaGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGtcatca > 2:62238/1‑149 (MQ=255) gCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGGGCACTGCGAGAGATGACACTCGACAGCCGGGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGtcacagggc > 2:266877/1‑143 (MQ=255) gCGAGGGGACAGGTGGCAGCTAGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACGCGACAGCCGTGTGGCTGCGGCGGGCGATCTGTTTGGAGCTGTAGTAGGTCATCAGGc > 2:303824/1‑149 (MQ=255) aGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGt < 2:78042/149‑1 (MQ=255) gTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGTCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCg < 1:320311/149‑1 (MQ=255) gTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCg > 2:230235/1‑149 (MQ=255) ggCAGCTCGTCATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGat < 2:432374/149‑1 (MQ=255) gCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGata < 1:208010/149‑1 (MQ=255) | TACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGATCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGATA > NC_000913/93030‑93317 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |