| 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. | |||||||||||
GCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCA > NZ_CP009273/3757352‑3757505 | gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc < 1:689453/90‑1 (MQ=35)gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc < 1:617848/90‑1 (MQ=35)gccctgcgtggtacggacggcctttccccacgcatcgtactgtgtcccgtttctgctgccgtccgTACCACGCAGGGCGGGCTAACgcgc < 1:550954‑M2/25‑1 (MQ=255) cgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAg < 2:720586/90‑1 (MQ=35) cgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAg < 2:619597/90‑1 (MQ=35) cctgcgtggtacggacggcctttccccacgcatcgtactgtgtcccgtttctgctgccgtccgTACCACGCAGGGCGGGCTAACGCGCAg < 2:287745‑M2/27‑1 (MQ=255) gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc < 2:574326/87‑1 (MQ=33) gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc > 1:574326/1‑87 (MQ=33) cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt < 1:721605/90‑1 (MQ=35) cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt < 2:836925/90‑1 (MQ=33) ttGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATgg < 1:412660/90‑1 (MQ=35) ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAAt < 1:205017/90‑1 (MQ=35) gAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGAc < 2:17920/90‑1 (MQ=35) cacgcatcgtactgtgtcccgtttctgctgccgtccgTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACggg > 1:328138‑M2/38‑90 (MQ=255) cacaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc < 1:361531/90‑1 (MQ=33) caGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAt < 1:389738/69‑1 (MQ=37) caGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAt > 2:389738/1‑69 (MQ=35) gTACGATGCGTGGGGAAAGGCCGTCCGTACCGCGCAGGGCGGGCTaa < 2:461716/47‑1 (MQ=255) gTACGATGCGTGGGGAAAGGCCGTCCGTACCGCGCAGGGCGGGCTaa > 1:461716/1‑47 (MQ=255) gTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt > 1:681773/1‑90 (MQ=33) tACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg > 1:656434/1‑90 (MQ=33) tgtgtcccgtttctgctgccgtccgTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa < 2:328138‑M2/65‑1 (MQ=255) aTGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCa > 2:158486/1‑90 (MQ=33) cGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTg > 1:833389/1‑90 (MQ=33) gTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa > 2:590768/1‑90 (MQ=33) gTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa > 2:641795/1‑90 (MQ=33) tGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc > 1:523885/1‑42 (MQ=38) tGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc < 2:523885/42‑1 (MQ=38) tGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt > 1:818539/1‑47 (MQ=38) tGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt < 2:818539/47‑1 (MQ=40) tGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaa > 1:541997/1‑90 (MQ=33) ggggAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaaa > 1:59925/1‑90 (MQ=33) aGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCGGGACGGGTCATCCGCCGGACCGGTGAAAAGGGCa > 1:327009/1‑90 (MQ=11) cgtccgTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCca < 1:467550‑M2/84‑1 (MQ=255) gtccgTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc > 2:776397‑M2/6‑67 (MQ=255) gtccgTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc < 1:776397‑M2/62‑1 (MQ=255) gTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcacca > 2:664600‑M2/2‑90 (MQ=255) | GCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCA > NZ_CP009273/3757352‑3757505 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ 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. |