| 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% | 223.2 / NA | 68 | 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 (36/32); total (36/32) | |||||||||||
ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT > NC_000913/4577829‑4578091 | acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg > 1:617420/1‑136 (MQ=255)acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg < 2:617420/136‑1 (MQ=255) ggATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg > 2:17605/1‑128 (MQ=255) ggATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg < 1:17605/128‑1 (MQ=255) cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg > 1:465388/1‑139 (MQ=255) gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 1:275386/139‑1 (MQ=255) cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt < 2:1597498/121‑1 (MQ=255) cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt > 1:1597498/1‑121 (MQ=255) tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTAt > 2:81984/1‑87 (MQ=255) tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTAt < 1:81984/87‑1 (MQ=255) tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaa > 1:201760/1‑139 (MQ=255) aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa > 1:1428853/1‑139 (MQ=255) cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg > 1:180220/1‑91 (MQ=255) cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg < 2:180220/91‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg > 1:1835972/1‑60 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg < 2:1835972/60‑1 (MQ=255) tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGCTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGatat > 2:226417/1‑139 (MQ=255) gCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 1:1800056/90‑1 (MQ=255) gCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 2:1800056/1‑90 (MQ=255) gTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat < 1:1867684/139‑1 (MQ=255) gCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 1:380639/80‑1 (MQ=255) gCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 2:380639/1‑80 (MQ=255) cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACaa > 1:233222/1‑139 (MQ=255) cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAATCTTCGAGTAAAGGGGCGATATCCGTCATCACaa < 2:1727729/139‑1 (MQ=255) aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt < 2:201760/139‑1 (MQ=255) gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttct > 1:177888/1‑105 (MQ=255) gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttct < 2:177888/105‑1 (MQ=255) gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAttt > 1:93985/1‑139 (MQ=255) gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAttt > 1:1915390/1‑139 (MQ=255) tAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg < 2:1656248/42‑1 (MQ=255) tAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg > 1:1656248/1‑42 (MQ=255) aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCa < 2:465388/139‑1 (MQ=255) aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCa > 2:76022/1‑139 (MQ=255) atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAggg > 1:945257/1‑110 (MQ=255) atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAggg < 2:945257/110‑1 (MQ=255) atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAg > 2:1351027/1‑139 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGt > 2:1043303/1‑35 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGt < 1:1043303/35‑1 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCta < 2:1982044/41‑1 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCta > 1:1982044/1‑41 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc < 2:1522063/121‑1 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc > 1:1522063/1‑121 (MQ=255) tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt > 1:503867/1‑126 (MQ=255) tACGGACAGGTATCACGGGCTGTGCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt < 2:503867/126‑1 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt < 1:1335464/55‑1 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt > 2:1335464/1‑55 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt < 2:567656/123‑1 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt > 1:567656/1‑123 (MQ=255) gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATcc > 2:674244/1‑139 (MQ=255) tATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCg < 2:1428853/139‑1 (MQ=255) cACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgag < 2:233222/139‑1 (MQ=255) ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgg > 2:849406/1‑88 (MQ=255) ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgg < 1:849406/88‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt > 1:541855/1‑41 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt < 2:541855/41‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 1:1075026/1‑46 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 2:1075026/46‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaaa < 1:1805888/65‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaaa > 2:1805888/1‑65 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat < 1:1359257/101‑1 (MQ=255) cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat > 2:1359257/1‑101 (MQ=255) tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt < 2:1893791/38‑1 (MQ=255) tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt > 1:1893791/1‑38 (MQ=255) tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt < 2:956773/45‑1 (MQ=255) tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt > 1:956773/1‑45 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa < 2:1061479/79‑1 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa > 1:1061479/1‑79 (MQ=255) ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAt > 1:558576/1‑139 (MQ=255) | ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT > NC_000913/4577829‑4578091 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |