Closed rdocking closed 6 years ago
Somewhat done via merge_annotations.sh
. However, this is solution only works pairwise (i.e. two bedpe files at a time). Moreover, if overlapping intervals exist, there is no function to merge these intervals into a single set of merged coordinates (a la bedtools merge
single observation bed files).
From the existing BEDPE-formatted annotation files, create a new BEDPE-formatted output file that indicates presence/absence for a given fusion in a given annotation source.
E.g., given, say, TCGA and CIViC annotation files, produce a file that looks like either:
Or: