Merging bam files
Hi, I'm trying to merge these two bam files and I did the command
samtools merge merged.bam file1.bam file2.bam
However, it gives me this error:
[trans_tbl_init] @SQ SN (hg38_gold start=chr14:105786368 strand=+) found in binary header but not text header.
I'm not sure how to get around this. Is there an option I need to turn on to suppress/fix this?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Issue with samtools reheader, bam file
written by nrk_02 •
-
Size of Output BAM file bigger than the SUM Size of Input files after merging with samtools
written by Giggle •
-
Having trouble with samtools cat
written by matthew.dambrosio.1 •
-
Reformat sam header issue
written by pg_canada •
-
GNU parallel on multiple files to re-header bam files
written by pierre.justeau •
-
samtools merge warning
written by lshepard
-
Error sorted paired-end .bam using samtools
written by holbrookjoanna •
-
error in removing some lines in bam file
written by ashkan
-
samtools merge error
written by Hatem Elshazly •
-
Samtools Help needed: index seg fault after reheader
written by akhst7 •
What version are you using?
Version 1.1. Is that too outdated?
Can you update to 1.2 and see if that fixes it?
I updated it and it gives me the same error.
I just checked the code. This should happen whenever you have a chromosome missing between the files that you want to merge. You might be able to get around this with
samtool merge merged.bam file2.bam file1.bam, but if not then it makes sense to reheader at least one of the files to have all of the chromosomes (put the missing ones at the end, don't change the order of any of the existing ones!) and then merge.Or you can try MergeSamFiles commans from picard tools. It is better in handling headers than samtools.
Indeed, picard tends to be a bit friendlier with merging headers.