| 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. | |||||||||||
GCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > NZ_CP009273/4567723‑4567814‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACGGATATCGCCCCTTTACTCGAAGAATATTTCTT < NZ_CP009273/4567918‑4567830 GCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTG < 2:483017/91‑1 GATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTGGGTTGGAAGATGGCACCT > 2:1057679/1‑91 CAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTG < 1:287619/60‑1 CAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTG > 2:287619/1‑60 GCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTGGGTTGGAAGATGGCACCTCT < 1:351588/70‑1 GCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTGGGTTGGAAGATGGCACCTCT > 2:351588/1‑70 GTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTG < 1:519444/45‑1 GTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTG > 2:519444/1‑45 TTGTGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACGGATATCGCCCCTTTACTCGAAGAATATTTC < 1:1057679/91‑1 GTGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACGGATATCGCCCCTTTACTCGAAGAATATTTCTT < 1:4654/91‑1 GTGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACG < 1:877160/59‑1 GTGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACG > 2:877160/1‑59 GCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > NZ_CP009273/4567723‑4567814‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtGGGTTGGAAGATGGCACCTCTCCGGATACGCAATGGCTTAATGAAATTGTGATGACGGATATCGCCCCTTTACTCGAAGAATATTTCTT < NZ_CP009273/4567918‑4567830 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ 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. |