Problem With Sorting Bam File: Samtools
1
0
Entering edit mode
9.3 years ago
lyfsa ▴ 30

Hi!

I have three aligned pair end reads in BAM format. Two of them are of same sample run in different lanes. I first merged the files from same sample with read groups and added read group to third one separately. After that I merge all three bam files into one, also changing header to required one. I then tried to sort the merged bam file. But this is what I get as a sorted file. File_RG.sorted.0000.bam , File_RG.sorted.0001.bam, File_RG.sorted.0002.bam and so on....

Why am I not able to sort the BAM file? I am using samtools for merging and sorting.

bam samtools merge • 5.6k views
ADD COMMENT
4
Entering edit mode
9.3 years ago
matted 7.5k

Are you waiting for the command to finish, or is it crashing early for some reason? samtools sort will make those intermediate files (0000, 0001, etc.) since it cannot hold all reads in memory during the sort. It will do a final mergesort that produces the single bam you want, and should then delete those intermediate files.

ADD COMMENT
1
Entering edit mode

You don't even need to delete them, they are merged back to the output file prefix you gave, while launching the command.

ADD REPLY
0
Entering edit mode

I did not wait for command to finish. I was not aware of this piece of information for sorting. Thank you for answering :)

ADD REPLY
0
Entering edit mode

I am having the same error .I am generating 25 samples ,some of the sorted bam files get automerged after sorting,but some of them dont get merged.

Should I use the Samtools "merge" command to merge those sorted bam files ? or I have to give memory options?

ADD REPLY

Login before adding your answer.

Traffic: 1786 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6