Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,468,126 | Δ1 bp | 8.5% | coding (814/1047 nt) | mltB ← | membrane‑bound lytic murein transglycosylase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,468,126 | 0 | T | . | 8.5% | 168.3 / 10.9 | 59 | coding (814/1047 nt) | mltB | membrane‑bound lytic murein transglycosylase B |
Reads supporting (aligned to +/- strand): ref base T (26/28); new base . (3/2); total (29/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01 |
GCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGG > NC_002947/5467981‑5468271 | gcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCg > 1:1080833/1‑150 (MQ=255) gcgTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCg < 2:1162024/150‑1 (MQ=255) gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAg < 1:47359/150‑1 (MQ=255) gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAg < 2:1055684/150‑1 (MQ=255) tcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGc < 1:696267/150‑1 (MQ=255) gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTCGGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAAGCCTtcgac > 2:303585/1‑148 (MQ=255) gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTtcgtc > 1:888777/1‑150 (MQ=255) ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCggtg > 2:322382/1‑134 (MQ=255) ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCggtg < 1:322382/134‑1 (MQ=255) aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCt > 2:71569/1‑150 (MQ=255) tACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCa > 2:67198/1‑150 (MQ=255) cTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCAcc > 2:1086189/1‑150 (MQ=255) tCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt > 1:1155303/1‑150 (MQ=255) ggggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTc < 1:71569/150‑1 (MQ=255) tctcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGACAGGCCTTAGTCGGCGCGCTCACCTTCCACCCAgg > 1:601485/1‑150 (MQ=255) tcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAg < 2:176063/116‑1 (MQ=255) tcGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAg > 1:176063/1‑116 (MQ=255) cccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGa < 2:1155303/150‑1 (MQ=255) cTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGAcg < 2:875482/150‑2 (MQ=255) cGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGa < 1:316643/150‑1 (MQ=255) gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAc > 2:1053357/1‑150 (MQ=255) gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTAACCTACCACCCAGGCGCAACTGAc < 2:838970/150‑1 (MQ=255) cGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACcaca < 1:67198/150‑1 (MQ=255) cGGCAGATCACCGCGCAGCGAATTATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCAc > 2:878139/1‑150 (MQ=255) ggCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCAcc > 2:783916/1‑150 (MQ=255) tcatcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCg < 2:257130/150‑1 (MQ=255) tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTgg < 2:1037489/150‑1 (MQ=255) tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTgg < 2:66388/150‑1 (MQ=255) gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt > 1:1409301/1‑100 (MQ=255) gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCt < 2:1409301/100‑1 (MQ=255) gcgcAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGgcg > 2:1258258/1‑150 (MQ=255) tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGCTGGTGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg < 1:1086189/150‑1 (MQ=255) tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg < 1:1053357/150‑1 (MQ=255) aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTtcgtc > 1:1072543/1‑71 (MQ=255) aCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTtcgtc < 2:1072543/71‑1 (MQ=255) aaGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt < 2:1240402/150‑1 (MQ=255) aaGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt > 2:50228/1‑150 (MQ=255) aGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCTATGCCTGTGGTCAGGCCTTCGTCGGCTCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg < 2:601485/150‑1 (MQ=255) aCGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATc < 1:1343770/150‑1 (MQ=255) gCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg > 2:402582/1‑125 (MQ=255) gCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg < 1:402582/125‑1 (MQ=255) cccTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg < 2:365682/112‑1 (MQ=255) cccTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTg > 1:365682/1‑112 (MQ=255) ccTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCCGCCTAGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttg > 1:87072/1‑150 (MQ=255) cTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgtttt > 2:1363050/1‑148 (MQ=255) cGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCa > 2:668385/1‑150 (MQ=255) gTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAg > 2:126985/1‑150 (MQ=255) gaccACCGGCATCGA‑GCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa > 1:706007/4‑149 (MQ=255) gaccACCGGCATCGA‑GCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGa < 2:706007/146‑1 (MQ=255) ccACCGGCATCGA‑GCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt < 1:1376499/123‑1 (MQ=255) ccACCGGCATCGA‑GCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt > 2:1376499/2‑124 (MQ=255) cACCGGCATCGA‑GCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt > 1:1209718/1‑150 (MQ=255) ccGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa > 2:10078/1‑150 (MQ=255) cGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCCCCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTTTTCCAGATATTGat < 1:1363050/150‑1 (MQ=255) cGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGat < 1:878139/150‑1 (MQ=255) gC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc < 1:1419096/114‑1 (MQ=255) gC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc > 2:1419096/1‑114 (MQ=255) gATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATgg > 2:29936/1‑148 (MQ=255) gATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATgg < 1:29936/148‑1 (MQ=255) | GCGTAGAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGC‑TCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGG > NC_002947/5467981‑5468271 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |