| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 4578582 | 54 (1.300) | 4 (0.100) | 4/270 | 5.6 | 7.2% | coding (756/1380 nt) | mcrB | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| ? | NC_000913 | = 4578585 | 52 (1.310) | 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/4578438‑4578582‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtcggCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < NC_000913/4578585‑4578448 ATCATGTTCCATTAACATCATCACTTCGCCAAACACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGCCAAAAATTGTAAAATATGCCGTCTTTACGTC < 1:576492/143‑1ATCATGTTCCATAAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCGCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTC < 1:208116/143‑1ATCATGTGCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTC < 1:208223/143‑1 gGTTCGATTAACATCATCACTTCGCCAAATACTTGACTGAGGTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTCACGTCcgaa < 2:132782/142‑5 ccagCATCAGTTCGGCACATACTTCAATGAGATTGGTACGATGGATTTAATCTATAATAAAAATATACGTTTTCACTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACGTAA < 2:42411/139‑1 ATACTGTACTGAGATTGGCACGATTGAGTTCATCGATAGTAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCGGACAAAAACTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTACAA < 1:73652/143‑1 GTAAAGTATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTAC < 1:478389/53‑1 GTAAAGTATGCCGTCTTTACGTCGGCTTCCGACGTAAAGACGGCATATTTTAC > 2:478389/1‑53 CGTCGGCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAA > 1:709784/1‑143 GTCGGCTGCGGACGGAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < 2:222432/143‑1 ATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > NC_000913/4578438‑4578582‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtcggCTTCCGACGTAAAGACGGCATATTTTACAATTTTTGTCAGCAAGCTAAAGAGCAGCCAGAGAAAAAGTATATTTTTATTATAGATGAAATCAATCGTGCCAATCTCAGTAAAGTATTTGGCGAAGTGATGATGTTAAT < NC_000913/4578585‑4578448 |
| 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. |