| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,577,958 | C→T | M1M (GTG→ATG) † *460* (TAG→TAA) |
mcrC ← mcrB ← |
5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,577,958 | 0 | C | T | 100.0% | 234.6 / NA | 72 | M1M (GTG→ATG) † *460* (TAG→TAA) | mcrC mcrB | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (40/32); total (40/32) | |||||||||||
ATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT > NC_000913/4577827‑4578091 | atacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTTTACCATTATGaat > 2:358099/1‑137 (MQ=255) acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg > 1:1153030/1‑136 (MQ=255) acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg < 2:1153030/136‑1 (MQ=255) aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc > 2:27893/1‑133 (MQ=255) aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc < 1:27893/133‑1 (MQ=255) aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaa > 1:245143/1‑139 (MQ=255) atCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg < 1:608377/126‑1 (MQ=255) atCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg > 2:608377/1‑126 (MQ=255) cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg < 1:1336458/124‑1 (MQ=255) cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg > 2:1336458/1‑124 (MQ=255) cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg > 2:1058455/1‑139 (MQ=255) cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg > 2:1201362/1‑139 (MQ=255) gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 2:245143/139‑1 (MQ=255) gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTACCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 1:358099/139‑1 (MQ=255) ggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgttg > 2:988305/1‑139 (MQ=255) cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc > 1:46559/1‑99 (MQ=255) cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc < 2:46559/99‑1 (MQ=255) ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaata > 2:128841/1‑101 (MQ=255) ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaata < 1:128841/101‑1 (MQ=255) tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc < 1:407934/92‑1 (MQ=255) tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc > 2:407934/1‑92 (MQ=255) tCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttcttc < 1:455770/139‑1 (MQ=255) aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTaa < 1:1201362/139‑1 (MQ=255) aCCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAgggg < 2:1331263/139‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcccc < 1:495390/65‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcccc > 2:495390/1‑65 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 2:1322799/91‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 1:1322799/1‑91 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt < 2:1252768/88‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt > 1:1252768/1‑88 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGatat > 1:412540/1‑139 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGatat > 1:800570/1‑139 (MQ=255) aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg > 1:932095/1‑53 (MQ=255) aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg < 2:932095/53‑1 (MQ=255) aaGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGtcatca < 2:800570/139‑1 (MQ=255) aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGat > 2:36209/1‑127 (MQ=255) aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGat < 1:36209/127‑1 (MQ=255) aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGtcatc < 2:233179/137‑1 (MQ=255) aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGtcatc > 1:233179/1‑137 (MQ=255) gCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATcaca > 2:464460/1‑139 (MQ=255) gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAAttt > 2:1057008/1‑139 (MQ=255) aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg < 2:94373/72‑1 (MQ=255) aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg > 1:94373/1‑72 (MQ=255) aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCa > 1:1193148/1‑139 (MQ=255) atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCAtt < 2:219292/136‑1 (MQ=255) atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCAtt > 1:219292/1‑136 (MQ=255) ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg > 2:685851/1‑64 (MQ=255) ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg < 1:685851/64‑1 (MQ=255) ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAgggg < 1:551450/108‑1 (MQ=255) ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAgggg > 2:551450/1‑108 (MQ=255) aCGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt > 1:1264967/1‑139 (MQ=255) cGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg < 2:1193148/139‑1 (MQ=255) aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCAtt < 1:1126056/125‑1 (MQ=255) aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCAtt > 2:1126056/1‑125 (MQ=255) aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTAt > 1:141859/1‑139 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATcc > 1:1060861/1‑139 (MQ=255) tCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGga < 2:1060861/139‑1 (MQ=255) cACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgag < 1:988305/139‑1 (MQ=255) cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGagag > 1:1237997/1‑139 (MQ=255) ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 2:1214773/48‑1 (MQ=255) ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 1:1214773/1‑48 (MQ=255) ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGt > 1:492281/1‑139 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcat > 2:1056818/1‑61 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcat < 1:1056818/61‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTg > 1:965946/1‑139 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAgggg < 2:118818/88‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAgggg > 1:118818/1‑88 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATGTCTGTTGTTTATAGGGGTCGTCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTGTCCGGAGAGGTGc > 2:944018/1‑139 (MQ=255) tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc < 2:589672/139‑1 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc > 1:913988/1‑137 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc < 2:913988/137‑1 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAt > 2:1087244/1‑139 (MQ=255) | ATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGA‑GTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT > NC_000913/4577827‑4578091 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |