This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Observing merged genomic intervals along with their original unmerged intervals

Hello everyone,

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.

Thank you

bedtools merging
• 1,100 views • link
updated by Kevin Blighe Kevin Blighe • written by arsala521 • arsala521

1 answer

There is definitely an awk solution here; however, you may try out BEDTools:

cat test.bed 
chr1    500 1000
chr1    750 1000
chr2    200 500

bedtools merge -i test.bed -c 1,2,3 -o distinct -delim " "
chr1    500 1000    chr1    500 750 1000
chr2    200 500     chr2    200 500

Give that a try.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this