Question: merge structural variant calls in bedpe format
0
gravatar for Ming Tang
3.3 years ago by
Ming Tang2.5k
Houston/MD Anderson Cancer Center
Ming Tang2.5k wrote:

Hi all,

I have structural calls represented in bedpe format http://bedtools.readthedocs.org/en/latest/content/general-usage.html#bedpe-format

I want to merge the calls with both ends overlapping. similar to this question Merging And De-Duplicating Structural Variant Calls (Bedpe)

I also read https://groups.google.com/forum/#!topic/bedtools-discuss/JXZbJSwVxUo

and it seems to be not an easy job.

Anyone has written tools for this kind of task?

I am aware of svtools https://github.com/hall-lab/svtools, but it is still under development, and not tested yet.

Thanks!

Ming

 

bedpe bedtools • 1.3k views
ADD COMMENTlink modified 3.3 years ago by leekaiinthesky170 • written 3.3 years ago by Ming Tang2.5k
0
gravatar for leekaiinthesky
3.3 years ago by
UCLA
leekaiinthesky170 wrote:

Does bedtools merge work for you?

http://bedtools.readthedocs.org/en/latest/content/tools/merge.html

And if not, what do you need it to do that it doesn't do? (An example might be helpful.)

 

ADD COMMENTlink written 3.3 years ago by leekaiinthesky170

bedpe1:

chr1 100 200 chr1 400 500

chr1 1000 1200 chr1 1500 1600

bedpe2

chr2 40 150   chr1 350 500
chr1 70 150  chr1 450 550

 

The first entry in bedpe1 overlaps with the second entry in bedpe2.

 

---------|---------|-----------------------|----------|-----    

-----|---------|-------------------------|---------------|---

merged to 

-----|-------------|---------------------|---------------|---

This is just an example for intra-chromosomal merging. 

ideally, one can specify, how merging should be performed. minimal overlappings:

--------|------|---------------------------|-----------|-----
ADD REPLYlink written 3.3 years ago by Ming Tang2.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2095 users visited in the last hour