Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,210,465 | C→T | Y40Y (TAC→TAT) | mcrA → | e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,210,465 | 0 | C | T | 100.0% | 41.9 / NA | 13 | Y40Y (TAC→TAT) | mcrA | e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA |
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (8/5); total (8/5) |
TGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCT > NC_000913/1210388‑1210533 | tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAat < 2:239185/81‑1 (MQ=255) aTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGct < 1:239186/81‑1 (MQ=255) aGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGctct > 1:239209/1‑81 (MQ=255) ttATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACg > 1:239204/1‑81 (MQ=255) gTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACggttgg < 1:239187/81‑1 (MQ=255) aTCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAt > 1:239205/1‑81 (MQ=255) gTGGGGGCCGGGTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATc < 1:239188/81‑1 (MQ=255) gggCCGGGTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa < 1:239189/81‑1 (MQ=255) gTGACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCgtag > 2:239212/1‑81 (MQ=255) gACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCgtagta > 1:239210/1‑81 (MQ=255) gACAGAAACAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCgtagta > 2:239211/1‑81 (MQ=255) aCAGAAACAAAGATTATAATATCGCTCTTGATTATATTTTTGAACGGTTGGTTGATTCCGGTGTATCCCAAGTCgtagta > 1:239207/1‑80 (MQ=255) aaaCAAAGATTATAATATCGCTCTTGATTATATCATTGAACGGTCGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCt > 1:239206/1‑81 (MQ=255) | TGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCT > NC_000913/1210388‑1210533 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |