| 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. | |||||||||||
AAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTT > NZ_CP009273/4567830‑4568007 | aaGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCcac > 1:4654/1‑91 (MQ=255) gAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCcacaa > 1:1057679/1‑90 (MQ=255) aaaTATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCACGCCAcagc < 1:127637/91‑1 (MQ=255) aaaTATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCAcagc < 1:73890/91‑1 (MQ=255) aaTATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCAcagca < 2:854604/91‑1 (MQ=255) tcttcGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGt < 1:145902/91‑1 (MQ=255) tcttcGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGt < 1:572930/91‑1 (MQ=255) gAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTaatat > 2:919385/1‑88 (MQ=255) aaGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGccc < 1:11358/91‑1 (MQ=255) gCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCaa > 1:105594/1‑88 (MQ=255) gCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCaa < 2:105594/88‑1 (MQ=255) gATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCgg > 2:97574/1‑91 (MQ=255) cGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCcac < 2:877160/59‑1 (MQ=255) cGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCcac > 1:877160/1‑59 (MQ=255) gTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 2:993148/91‑1 (MQ=255) gTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 2:4654/91‑1 (MQ=255) gTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 1:934995/91‑1 (MQ=255) gTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 1:914507/91‑1 (MQ=255) tcatcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCACCCCACAGCAGAAGAAACTATGCCCAATGCGGAATCtttc > 1:525669/1‑90 (MQ=255) tcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCAcagc < 1:601662/57‑1 (MQ=255) tcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCAcagc > 2:601662/1‑57 (MQ=255) tAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGc < 1:385233/91‑1 (MQ=255) aTTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCttt < 1:997590/91‑1 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCaa < 1:945688/56‑1 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCaa > 2:945688/1‑56 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 1:947945/68‑1 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt > 2:947945/1‑68 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCa > 1:859435/1‑73 (MQ=255) ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCa < 2:859435/73‑1 (MQ=255) tGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGc < 1:687349/91‑1 (MQ=255) tATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGAt > 2:941395/1‑91 (MQ=255) tATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGAt < 2:866806/91‑1 (MQ=255) aTCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTGTCCCAAGGATAGTGGCCTCTTTGCTGAtt > 2:670123/1‑91 (MQ=255) gTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt < 2:598021/52‑1 (MQ=255) gTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCtt > 1:598021/1‑52 (MQ=255) aTCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCa > 1:203260/1‑86 (MQ=255) aTCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCa < 2:203260/86‑1 (MQ=255) ccAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTcc > 2:202222/1‑74 (MQ=255) ccAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTcc < 1:202222/74‑1 (MQ=255) ccAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCAt > 2:89792/1‑91 (MQ=255) ccACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCAtttttt > 1:196990/1‑91 (MQ=255) | AAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTT > NZ_CP009273/4567830‑4568007 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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. |