This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bam file merge

Hello. I have numerous bam files which I need to merge into a single bam file. How do I know whether I need to first index and sort the bam files?

bam merge
• 3,431 views • link
updated by h.mon h.mon • written by genomics Newbie • genomics Newbie

In http://www.htslib.org/doc/samtools-1.2.html, it says,

-h is specified the @SQ headers of input files will be merged into the specified header, otherwise they will be merged into a composite header created from the input headers. If in the process of merging @SQ lines for coordinate sorted input files, a conflict arises as to the order (for example input1.bam has @SQ for a,b,c and input2.bam has b,a,c) then the resulting output file will need to be re-sorted back into coordinate order.

• 1 views • link
written by tiago211287 tiago211287

Dear genomics Newbie, Hi

Please have a look at To sort or not to Sort?.

I guess you can first use "samtools merge" and then "samtools sort -o".

• 1 views • link
written by Farbod Farbod

1 answer

You should sort your files, as samtools merge expects sorted bam files, either by position (default) or by name (samtools merge -n). You can check if you bam are already sorted by looking at the SO flag:

samtools view -H file.bam | grep "SO:"

@HD VN:1.0 SO:coordinate

• 1 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this