Biostar Beta. Not for public use.
Bed File Bedpe Format
1
Entering edit mode
16 months ago
Assa Yeroslaviz ♦ 1.2k
Munich

Hi,

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

workflow:
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?

Thanks

A.

ADD COMMENTlink
2
Entering edit mode
4.7 years ago
iw9oel_ad 6.0k

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
0
Entering edit mode

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

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1