Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,462,441 | Δ1 bp | 9.4% | intergenic (+16/+18) | PP_4798 → / ← PP_4799 | membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,462,441 | 0 | T | . | 9.4% | 163.6 / 11.3 | 53 | intergenic (+16/+18) | PP_4798/PP_4799 | membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase |
Reads supporting (aligned to +/- strand): ref base T (22/26); new base . (3/2); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01 |
GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACACC > NC_002947/5462294‑5462590 | gatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgc > 1:969039/1‑150 (MQ=255) cGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCgg < 1:187990/150‑1 (MQ=255) gCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCAc < 2:1071728/150‑1 (MQ=255) aTCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAGGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCt < 2:466938/150‑1 (MQ=255) aTCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCt < 1:109657/150‑1 (MQ=255) gCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAAtc < 1:1028649/135‑1 (MQ=255) gCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAAtc > 2:1028649/1‑135 (MQ=255) aTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGAt < 2:249480/148‑1 (MQ=255) aTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGAt > 1:249480/1‑148 (MQ=255) tACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATcc > 2:1129488/1‑149 (MQ=255) tACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATcc < 1:1129488/149‑1 (MQ=255) aaTGCGCAAGGCCCAGCAAGGGCTGGGATGGCCGACGGATGGCTAGCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGGCGCTGC‑GTGCCCCCGGTAATCTCACGCCCTGAACAGACCATGCTCCCGCACCa > 2:86488/1‑150 (MQ=255) cccGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAgt < 2:333158/149‑1 (MQ=255) gCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAATCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGtt > 2:102053/1‑150 (MQ=255) gCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGtt > 2:929921/1‑150 (MQ=255) aTGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGc < 2:276808/150‑1 (MQ=255) tGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGC‑GCAGCGGcccc > 1:638430/1‑102 (MQ=255) tGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGC‑GCAGCGGcccc < 2:638430/102‑1 (MQ=255) ccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGaa > 1:386308/1‑150 (MQ=255) cAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAac > 1:67911/1‑150 (MQ=255) gCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAacac > 2:20447/1‑150 (MQ=255) gATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGc > 2:297849/1‑147 (MQ=255) gATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGc < 1:297849/147‑1 (MQ=255) cTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg < 2:464413/150‑1 (MQ=255) cTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg < 1:20447/150‑1 (MQ=255) tACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAgggg > 1:218745/1‑150 (MQ=255) aCCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGc > 2:646519/1‑150 (MQ=255) cccGACCCACAAGCTGCTTGAGGGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGGTCCAGCAGCAATGGTTGGTGAACACTGTCCAACATCACGCGGgccgccggaggaa > 2:542060/1‑143 (MQ=255) cccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCa > 1:835993/1‑150 (MQ=255) ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCaa < 2:1113548/150‑1 (MQ=255) cGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAg < 2:386308/150‑1 (MQ=255) gACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGc > 1:547064/1‑150 (MQ=255) aCCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa < 1:929921/150‑1 (MQ=255) cacaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATAc < 2:592374/150‑1 (MQ=255) gcaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCAGTAATCTCACGCCCTGCACAGATCCTGCTCCAGAACCAGTGTTTGCTGAACACTGTCCAGCCTCACTCGAGCCCCCAGGGGCAAGCATACa < 1:6999/149‑1 (MQ=255) aaGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGGGCCCCCAGGGGCAAGCATAccta > 2:6999/1‑147 (MQ=255) aGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCATCAGCCGAGCCCCCAGGGGCAAGCATACATTc < 1:1125012/150‑1 (MQ=255) gctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGaa < 2:1194257/73‑1 (MQ=255) gctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGaa > 1:1194257/1‑73 (MQ=255) ctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg < 2:863987/150‑1 (MQ=255) ctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg < 1:171481/150‑1 (MQ=255) agCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTg > 1:641910/1‑150 (MQ=255) cagcagTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTgcg > 1:159859/1‑150 (MQ=255) agcagTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa < 1:766568/122‑1 (MQ=255) agcagTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa > 2:766568/1‑122 (MQ=255) agcagTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTgcgc < 2:625989/150‑1 (MQ=255) gTAATCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAc < 1:102053/150‑1 (MQ=255) aTTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCa < 1:313384/150‑1 (MQ=255) ccGGGGCCGC‑GCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTTCAGCCTCACCCGAGCCCCCAGGGGCaa < 2:605873/106‑1 (MQ=255) ccGGGGCCGC‑GCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTTCAGCCTCACCCGAGCCCCCAGGGGCaa > 1:605873/1‑106 (MQ=255) ccGGGGCCGC‑GCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCa > 2:134271/1‑150 (MQ=255) ggggCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAg > 2:340290/1‑91 (MQ=255) ggggCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAg < 1:340290/91‑1 (MQ=255) tGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACAcc > 2:35016/1‑150 (MQ=255) | GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGC‑GCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACACC > NC_002947/5462294‑5462590 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |