| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 4578583 | 54 (1.250) | 4 (0.100) | 4/270 | 5.8 | 7.2% | coding (756/1380 nt) | mcrB | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| ? | NC_000913 | = 4578586 | 52 (1.250) | coding (753/1380 nt) | mcrB | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB | |||||
| Rejected: Coverage evenness skew score above cutoff. | |||||||||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
ATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > NC_000913/4578439‑4578583‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtcggCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < NC_000913/4578586‑4578449 ATCATGTTCCATTAACATCATCACTTCGCCAAACACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGCCAAAAATTGTAAAATATGCCGTCTTTACGTC < 2:576492/143‑1ATCATGTTCCATAAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCGCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTC < 2:208116/143‑1ATCATGTGCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTC < 2:208223/143‑1 gGTTCGATTAACATCATCACTTCGCCAAATACTTGACTGAGGTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTCACGTCcgaa < 1:132782/142‑5 ccagCATCAGTTCGGCACATACTTCAATGAGATTGGTACGATGGATTTAATCTATAATAAAAATATACGTTTTCACTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACGTAA < 1:42411/139‑1 ATACTGTACTGAGATTGGCACGATTGAGTTCATCGATAGTAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCGGACAAAAACTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTACAA < 2:73652/143‑1 GTAAAGTATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTAC > 1:478389/1‑53 GTAAAGTATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTAC < 2:478389/53‑1 CGTCGGCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAA > 2:709784/1‑143 GTCGGCTGCGGACGGAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < 1:222432/143‑1 ATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > NC_000913/4578439‑4578583‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtcggCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < NC_000913/4578586‑4578449 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ 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. |