You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
A user tried to merge two SAM files from different jobs that were mapped against the same genome. The "samtools merge" command failed, complaining about incompatible headers in the SAM files.
I noticed that each file had a few chromosomes listed in the headers that were not listed in the other file. I advised them to modify the headers so they had the same chromosomes listed, and this fixed the problem.
It looks like we only create header lines for chromosomes that actually have reads that map too them. While this seems sensible, I think in order to have SAM files that are compatible with each other we should probably just create a header line for each annotated chromosome.
The text was updated successfully, but these errors were encountered:
A user tried to merge two SAM files from different jobs that were mapped against the same genome. The "samtools merge" command failed, complaining about incompatible headers in the SAM files.
I noticed that each file had a few chromosomes listed in the headers that were not listed in the other file. I advised them to modify the headers so they had the same chromosomes listed, and this fixed the problem.
It looks like we only create header lines for chromosomes that actually have reads that map too them. While this seems sensible, I think in order to have SAM files that are compatible with each other we should probably just create a header line for each annotated chromosome.
The text was updated successfully, but these errors were encountered: