| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | AP009048 | = 4590440 | 117 (1.440) | 5 (0.070) | 4/88 | NT | 4.3% | coding (1056/3567 nt) | hsdR | endonuclease R |
| ? | AP009048 | = 4590455 | 115 (1.570) | coding (1041/3567 nt) | hsdR | endonuclease R | |||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
CTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > AP009048/4590368‑4590440‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcgCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACCGGATGAAGTG < AP009048/4590455‑4590386 CTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTAT > 1:54467/1‑71CTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTAT > 1:4010097/1‑71CTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTA < 1:6434199/70‑1 CTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATC > 1:277827/1‑70 CACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCG < 1:4410878/55‑1 CATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGCTACAGTAAATGTTTCGATA < 1:763249/70‑1 ATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGCTACAGTAAATGTTTCGATAA < 1:1258644/70‑1 ATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGCTACAGTAAATGTTTCGATAA < 1:2124185/70‑1 CGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGCTACAGTAAATGTTTCGATAAT > 1:1285324/1‑68 TCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGCTACAGTAAATGT < 1:3825518/46‑1 ATCGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAAACTTGCTTGAGCACTACTCACCGGATG > 1:5841177/1‑69 ATCGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAA > 1:4989350/1‑42 TCGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACCGGATGAAG > 1:4944950/1‑71 TCGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACC < 1:3259614/63‑1 TCGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCAC < 1:4977145/55‑1 CGCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACCGGATGAAGT > 1:2068655/1‑71 GCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACCGGATGAAGTG < 1:7829127/71‑1 CTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > AP009048/4590368‑4590440‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcgCTACAGTAAATGTTTCGATAATGGCTTCCTGCGGGAAACCTTGCTTGAGCACTACTCACCGGATGAAGTG < AP009048/4590455‑4590386 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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. |