Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,111,290 | 0 | C→G | 100.0% | 0.2 | 2 | R29P (CGC→CCC) | insF‑1 | putative transposase for IS3 |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/0): ref base (0/0): total (2/0) |
AAGCCGCGAGGACAACGCTGTCGCAGTGTTGGCGGAACTGC > REL606/1111258‑1111298 | aaGCCGCGAGGACAACGCTGTCGCAGGGGTGGGgg > 1:5223800/1‑35 (MQ=1) aGCCGCGAGGACAACGCTGTCGCAGTGTTGGGggg > 1:18658744/1‑34 (MQ=0) aGCCGCGAGGACAACGCTGTCGCAGGGTTGGGggg > 1:14567565/1‑34 (MQ=1) aGCCGCGAGGACAACGCTGTCGCAGGGGTGGGggg > 1:333615/1‑34 (MQ=1) gCCGCGAGGACAACGCTGTCGCAGTGTTGGGggga > 1:1968557/1‑33 (MQ=1) gCCGCGAGGACAACGCTGTCGCAGTGTTGGGggga > 1:3496289/1‑33 (MQ=1) ccGCGAGGACAACGCTGTCGCAGTGTTGGGgggaa > 1:4593154/1‑32 (MQ=1) ccGCGAGGACAACGCTGTCGCAGTGTTGGGGGGAc > 1:3933388/1‑35 (MQ=1) ccGCGAGGACAACGCTGTCGCAGGGTTGGGgggaa > 1:5211066/1‑32 (MQ=1) cgAGGACAACGCTGTCGCAGTGTTGGGGGGAatgc > 1:3051971/1‑35 (MQ=1) cgAGGACAACGCTGTCGCAGTGTTGGCGGGAatgc > 1:7794829/1‑35 (MQ=1) cgAGGACAACGCTGTCGCAGGGTTGGCGGGAatgc > 1:2731768/1‑35 (MQ=1) | AAGCCGCGAGGACAACGCTGTCGCAGTGTTGGCGGAACTGC > REL606/1111258‑1111298 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |