New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | = 2680447 | 34 (1.540) | 3 (0.150) | 3/134 | 2.6 | 10.8% | coding (1519/1557 nt) | ECB_02510 | Recombinase |
? | REL606 | = 2680655 | 19 (0.960) | coding (181/897 nt) | ECB_02511 | ParB family protein | |||||
Rejected: Frequency below/above cutoff threshold. |
AAAGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAGAGA > REL606/2680581‑2680728 | aaaGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTt < 2:427329‑M2/76‑2 (MQ=255) aaaaGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTc > 2:328131/1‑76 (MQ=255) aTTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatca > 2:237780‑M2/1‑56 (MQ=255) aTTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGgaga > 1:480004/1‑76 (MQ=255) ttAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGaa > 1:539329/1‑76 (MQ=255) tAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAAt < 2:419303/76‑1 (MQ=255) tAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAAt < 1:608510/76‑1 (MQ=255) agcaagcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatcaagat > 2:95478‑M2/1‑52 (MQ=255) caagcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatcaagatta > 2:64102‑M2/1‑50 (MQ=255) gcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGct > 1:335960/1‑76 (MQ=255) tataAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCtctc < 1:402698/76‑1 (MQ=255) taAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGt > 2:339750/1‑76 (MQ=255) gTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCAt > 2:408009/1‑76 (MQ=255) gTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCAt > 2:408012/1‑76 (MQ=255) agagCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAg < 2:400013/76‑1 (MQ=255) gatgaAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACta > 2:657896/1‑76 (MQ=255) gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg < 1:443985/76‑1 (MQ=255) gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg < 2:603335/76‑1 (MQ=255) gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg < 1:147845/76‑1 (MQ=255) aaaTATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATta < 2:681966/76‑1 (MQ=255) aaTATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATtat > 1:437444/1‑76 (MQ=255) gCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAagag < 1:275616/76‑1 (MQ=255) cTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAgaga > 1:559770/1‑76 (MQ=255) | AAAGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAGAGA > REL606/2680581‑2680728 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ 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. |