Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,517,528 | C→T | 7.8% | R240C (CGT→TGT) | murE → | UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,517,528 | 0 | C | T | 7.8% | 141.6 / 2.2 | 51 | R240C (CGT→TGT) | murE | UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase |
Reads supporting (aligned to +/- strand): ref base C (25/22); new base T (2/2); total (27/24) | |||||||||||
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.24e-01 |
GGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGT > NC_002947/1517380‑1517673 | ggAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGgtc > 4:100531/1‑149 (MQ=255) aTTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGgtcgtc < 7:315198/149‑1 (MQ=255) tCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTAAGAGGCCGCCAAGGCCAAGCTGTCACGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTggc < 1:272268/149‑1 (MQ=255) gcgGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCCTCTggcggc > 4:93415/1‑149 (MQ=255) cgGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGcc > 4:5657/1‑149 (MQ=255) aCCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgc < 5:272404/149‑1 (MQ=255) aCCAACCTGTCCCGCGACCACCTCGACTATCACCGCAGCATGGAAGCCTATGATGCCGCCAAGGCCAAGCTGT‑TAGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgc > 6:192429/1‑149 (MQ=255) cAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgcgc > 1:34496/1‑149 (MQ=255) cAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGACTGGCGGCCGACTTcgcgc > 3:45819/1‑149 (MQ=255) ccGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCa < 3:100531/148‑1 (MQ=255) cGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTGGAccttt > 8:254443/1‑146 (MQ=255) cGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCata > 8:259139/1‑149 (MQ=255) aCTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc > 1:206864/1‑149 (MQ=255) tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCg > 4:166978/1‑148 (MQ=255) tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGCCTGGCGGCCGACTTCGCGCGTCGCCCaa > 4:116925/1‑133 (MQ=255) tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGCCTGGCGGCCGACTTCGCGCGTCGCCCaa < 3:116925/133‑1 (MQ=255) tCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGaaa > 3:208142/1‑149 (MQ=255) cACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAAc > 8:101005/1‑149 (MQ=255) cACGGCAGCATGGAAGCCTATGAGGCCGCCAACGCCAAGCAGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCCTATCGAAAc > 7:211654/1‑149 (MQ=255) aTGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGCCTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGgctgct > 7:289465/1‑149 (MQ=255) aTGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGgctgct > 4:63251/1‑149 (MQ=255) gAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGa < 7:101005/148‑1 (MQ=255) aaGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGc > 7:117756/1‑149 (MQ=255) aaGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGc < 4:208142/149‑1 (MQ=255) gCCTCTGCGGCCGCCAAGGCCATGCTGT‑TCGCCTGGCCGTGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTa < 8:211654/149‑1 (MQ=255) cTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc > 8:206996/1‑127 (MQ=255) cTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc < 7:206996/127‑1 (MQ=255) cTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACa < 3:63251/149‑1 (MQ=255) cTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACa < 3:166978/149‑1 (MQ=255) tATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAg > 5:277037/1‑149 (MQ=255) aGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCt < 2:206864/148‑1 (MQ=255) aGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTg < 3:5657/149‑1 (MQ=255) caagCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTGCAGCCTGGAAAGCCCAGAAg > 4:83339/1‑149 (MQ=255) cTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTc < 7:259139/149‑1 (MQ=255) cTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTc < 7:254443/149‑1 (MQ=255) t‑tCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgtt < 2:34496/149‑1 (MQ=255) t‑tCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgt < 3:93415/148‑1 (MQ=255) cagcttGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgttg < 7:201665/144‑1 (MQ=255) ccTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTc < 8:117756/149‑1 (MQ=255) gCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGcc < 5:192429/149‑1 (MQ=255) ccTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAg > 7:231511/1‑149 (MQ=255) ccTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAg > 8:141063/1‑149 (MQ=255) tGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCTTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCAAGAAGCCTAGTTGTTCTGCCGAGAAGcc > 5:319907/1‑149 (MQ=255) cGTTGCCAGTTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCg < 8:231511/148‑1 (MQ=255) ccAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTAtt > 2:341352/1‑147 (MQ=255) gtCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGa > 6:37248/1‑149 (MQ=255) aaCCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATg < 1:341352/148‑1 (MQ=255) ccTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCg > 2:165737/1‑149 (MQ=255) ggATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTgc > 7:257526/1‑149 (MQ=255) atgatgACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTgcgc < 6:277037/149‑1 (MQ=255) tgatgaCTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCg < 6:111898/106‑1 (MQ=255) cgatgaCTTGGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCg < 6:111899/105‑1 (MQ=255) tttGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCt > 8:257170/1‑149 (MQ=255) gtcgtcgtcTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGt < 7:141063/148‑1 (MQ=255) | GGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGT‑TCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGT > NC_002947/1517380‑1517673 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |