Question: Bed File Bedpe Format
gravatar for Assa Yeroslaviz
9.2 years ago by
Assa Yeroslaviz1.4k
Assa Yeroslaviz1.4k wrote:


I'm having trouble with converting the bam file into bed -bedpe using the bedtools.

samtools sort -n mut.bam mut.Namesorted
bamTobed -i mut.Namesorted.bam -bedpe > dilpMerged_bedpe.bed

After sorting the file by read name (option -n) I run the bamTobed command. but it gives me an error message after running a few lines:

*ERROR: -bedpe requires BAM to be sorted/grouped by query name.

What am I doing wrong here?



bedtools bed bam conversion • 7.2k views
ADD COMMENTlink modified 2.6 years ago by Biostar ♦♦ 20 • written 9.2 years ago by Assa Yeroslaviz1.4k
gravatar for biobot 0.0.77.a.1099
9.2 years ago by
biobot 0.0.77.a.10996.1k wrote:

I suspect that one of two things is happening: either a) your SAM header (which will be unchanged by samtools) says that your BAM file is coordinate sorted (or unsorted) after your sort procedure and BEDTools is noticing this or 2) you have orphan (i.e. unpaired reads) and BEDTools is using a naive pairwise iteration that does not account for intervening orphans.

I can discount 1) by looking at the source (confirmed by experiment).

I can confirm 2) by experiment. I removed a single member of a pair of reads in a name-sorted BAM file and obtained the same error message from bamToBed as you did.

So I suspect that you have orphans that you need to remove.

ADD COMMENTlink modified 9.2 years ago • written 9.2 years ago by biobot 0.0.77.a.10996.1k

That's right, the current version expects alignments to come in pairs. I've not yet implemented logic for interleaved orphans.

ADD REPLYlink written 9.0 years ago by Aaronquinlan11k
Please log in to add an answer.


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