New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP009273 | 3614374 = | 6 (0.210) | 4 (0.150) | 4/154 | 1.8 | 18.9% | coding (1823/4236 nt) | rhsB | rhs element protein RhsB |
? | NZ_CP009273 | 3757417 = | 29 (1.100) | coding (1875/4134 nt) | rhsA | RHS element protein RhsA | |||||
Rejected: Frequency below/above cutoff threshold. |
TGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAT > NZ_CP009273/3614285‑3614460 | tggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtac > 2:664600‑M1/90‑90 (MQ=255) tggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtacggacg < 1:467550‑M1/6‑1 (MQ=255) tcaggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGAcaca < 2:328138‑M1/25‑1 (MQ=255) ggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtacggac < 1:776397‑M1/5‑1 (MQ=255) ggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtacggac > 2:776397‑M1/63‑67 (MQ=255) cccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTg > 1:328138‑M1/54‑90 (MQ=255) ctgcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAgg < 2:287745‑M1/63‑1 (MQ=255) gcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAgggc < 1:550954‑M1/65‑1 (MQ=255) gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc < 1:220394/90‑1 (MQ=35) gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc < 1:419143/90‑1 (MQ=35) gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc < 2:656888/90‑1 (MQ=35) tCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa < 1:430126/90‑1 (MQ=35) cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTAt < 1:743304/90‑1 (MQ=35) ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAt < 2:579511/90‑1 (MQ=35) | TGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAT > NZ_CP009273/3614285‑3614460 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ 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. |