New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP009273 | = 4567814 | 30 (0.970) | 4 (0.130) | 4/154 | 3.0 | 10.8% | coding (1341/1380 nt) | mcrB | 5‑methylcytosine‑specific restriction endonuclease subunit McrB |
? | NZ_CP009273 | = 4567918 | 37 (1.240) | coding (1237/1380 nt) | mcrB | 5‑methylcytosine‑specific restriction endonuclease subunit McrB | |||||
Rejected: Coverage evenness skew score above cutoff. | |||||||||||
Rejected: Frequency below/above cutoff threshold. |
AGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGA > NZ_CP009273/4567727‑4567900 | aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta > 1:486063/1‑91 (MQ=255) ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat < 1:1038317/91‑1 (MQ=255) aaGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAgg > 1:57828/1‑91 (MQ=255) tatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGaaa > 1:854604/1‑91 (MQ=255) ttACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAata > 1:993148/1‑91 (MQ=255) ggTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat < 2:1035003/63‑1 (MQ=255) ggTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat > 1:1035003/1‑63 (MQ=255) gggCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg > 1:388027/1‑91 (MQ=255) gCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg > 2:685931/1‑91 (MQ=255) gCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg > 2:385233/1‑91 (MQ=255) gCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg > 2:106005/1‑91 (MQ=255) gCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg > 1:623766/1‑91 (MQ=255) cTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGa > 1:239666/1‑91 (MQ=255) aTGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat > 2:687349/1‑91 (MQ=255) gAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaa > 2:775784/1‑52 (MQ=255) gAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaa < 1:775784/52‑1 (MQ=255) tCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACaa > 2:542989/1‑91 (MQ=255) cccTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc < 1:570653/85‑1 (MQ=255) cccTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc > 2:570653/1‑85 (MQ=255) ccTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAttt > 1:608976/1‑91 (MQ=255) taatTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTa < 1:542989/91‑1 (MQ=255) ttgttgGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAtt > 1:228719/1‑81 (MQ=255) ttgttgGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAtt < 2:228719/81‑1 (MQ=255) ttgttgGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc < 1:685931/91‑1 (MQ=255) tgGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg > 2:934995/1‑91 (MQ=255) ggTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCAAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc > 2:436151/1‑91 (MQ=255) ccATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTa < 2:388576/91‑1 (MQ=255) tgttgtTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgag < 2:608976/91‑1 (MQ=255) tgttgtTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgag < 1:106005/91‑1 (MQ=255) gttgttTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgaga < 2:629826/91‑1 (MQ=255) | AGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGA > NZ_CP009273/4567727‑4567900 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ 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. |