How to merge several BED files to create a comprehensive overview file
I have a few .bed files from different companies offering exome sequencing kits.
I would like to have a file that summarizes all target regions for all these kits. The .bed file have a basic structure composed for three columns (chr#, Start, End).
I would like to get an output table that shows which genomic regions are covered only by one of these kits, and which regions are covered by more than one (and which ones).
The best way to illustrate this is by an example:
I would like to merge and intersect the files for an output that divides the regions into subregions based on overlap between the input files. It should looks something like this:
|chr#||Start||End||Kit 1||Kit 2||Kit 3|
I would prefer to do this in python, but I could try in R as well. I have created a #pandas dataframe containing ranges from all kits, ordered by 'chr#' and 'Start' coordinates, which looks like this:
Any help would be appreciated.
• 32 views