Is there any option in bedtools or bedops that after merging show the original regions of the merged interval as well?
For example for the following input
chr1 500 1000 chr1 750 1000 chr2 200 500
can I have the output somewhat like this
chr1 500 1000 chr1 500 1000 chr1 750 1000 chr2 200 500
For every merged region, I want to see the original regions along with it.