Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,076,841 | G→C | 7.2% | R72P (CGC→CCC) | mreC → | cell shape‑determining protein MreC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,076,841 | 0 | G | C | 7.2% | 173.1 / 3.1 | 56 | R72P (CGC→CCC) | mreC | cell shape‑determining protein MreC |
Reads supporting (aligned to +/- strand): ref base G (26/25); new base C (2/2); total (28/28) | |||||||||||
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.89e-01 |
GGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAA > NC_002947/1076698‑1076988 | gggCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCa < 2:35564/150‑1 (MQ=255) gCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTTCTGCAGGCGCCCCTGCAgc > 1:458014/1‑149 (MQ=255) gTGCTGATGGAGTCGTACTGGATCATCCGCTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTTCTGCAGGCGCCCCTGCagcagcag < 1:1071579/150‑3 (MQ=255) tGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcgg > 2:941590/1‑150 (MQ=255) aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc > 1:1060206/1‑150 (MQ=255) aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc > 2:465924/1‑150 (MQ=255) tGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGccc < 2:48393/150‑1 (MQ=255) ggAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGACGCTAAAGACCGAAGCCCTGCTGCTGCAGGGGCTCCTGCAGAAGCTGGCGGCCCt < 2:1231606/150‑1 (MQ=255) gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGaa > 1:1288848/1‑150 (MQ=255) aTCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAGACGATCAGCACAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAAc > 1:798845/1‑150 (MQ=255) aCCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCAGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTg > 1:104294/1‑150 (MQ=255) cGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCg > 2:1126262/1‑150 (MQ=255) ccACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTgcgc > 1:1399019/1‑150 (MQ=255) acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGa > 1:425406/1‑150 (MQ=255) gCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGtt < 1:289670/150‑1 (MQ=255) aGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTc < 2:1399019/150‑1 (MQ=255) tgtgGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGGAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGgc > 1:329016/1‑150 (MQ=255) tgtgGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGgc > 1:1005692/1‑150 (MQ=255) tgtgGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCCAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGgc > 1:889864/1‑150 (MQ=255) gccggccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTg < 1:485204/150‑1 (MQ=255) cggccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGTGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGt > 2:664259/1‑150 (MQ=255) ccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGa < 1:411847/102‑1 (MQ=255) ccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGa > 2:411847/1‑102 (MQ=255) cAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggc > 1:1225479/1‑89 (MQ=255) cAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggc < 2:1225479/89‑1 (MQ=255) cAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAAc < 2:1288848/150‑1 (MQ=255) aGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGtt < 1:1126262/150‑1 (MQ=255) cGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAg > 1:781761/1‑111 (MQ=255) cGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAg < 2:781761/111‑1 (MQ=255) cACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGt > 1:207812/1‑150 (MQ=255) aCCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTg > 2:1131024/1‑150 (MQ=255) ccGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTgg > 2:322927/1‑150 (MQ=255) cGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGTGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGc < 1:664259/150‑1 (MQ=255) aGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCg > 2:1253833/1‑150 (MQ=255) gCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGa < 1:1326583/150‑1 (MQ=255) tGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGc < 1:941590/150‑1 (MQ=255) tGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCACAACGTGCGGCTGCGCGAGTTGCTCACCTCTTCGGCGCTGGTCAACGACAAGGTTGAAGTGGCCCAGc < 2:433777/150‑1 (MQ=255) aTTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAg > 1:451006/1‑139 (MQ=255) aTTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAg < 2:451006/139‑1 (MQ=255) ttGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGTGCAGAACGTGCGGCTGCGCTAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGa < 1:465924/150‑1 (MQ=255) gCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATc < 2:207812/150‑1 (MQ=255) gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGgt < 1:1253833/150‑1 (MQ=255) aaaaCGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGt > 2:1292279/1‑93 (MQ=255) aaaaCGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGt < 1:1292279/93‑1 (MQ=255) aaaCGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGgtgt > 1:1030743/1‑150 (MQ=255) gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTgcgc < 1:801192/84‑1 (MQ=255) gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTgcgc > 2:801192/1‑84 (MQ=255) gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCg > 2:449331/1‑142 (MQ=255) gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCg < 1:449331/142‑1 (MQ=255) cTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTTTCGACCCCTAc < 2:798845/150‑1 (MQ=255) ccGAAGCCCTGCTGCCGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTa < 1:161653/150‑1 (MQ=255) gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa > 2:84927/3‑150 (MQ=255) gCCCTGCTGCTGCAGGGGCCCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCAc > 1:1294479/1‑150 (MQ=255) tgctgcAGGGGCGCCTGCAGACCCTGACGGCCGTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTc < 2:889864/150‑1 (MQ=255) gcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGAt < 2:329016/150‑1 (MQ=255) cAGGCGCCCCTGCAGCAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAg < 1:630777/105‑1 (MQ=255) gcgcCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATcaacaa < 1:1131024/150‑1 (MQ=255) | GGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCTGATCAACAA > NC_002947/1076698‑1076988 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |