Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SAM headers should be changed to contain all chromosomes #163

Open
mdelaurentis opened this issue Feb 13, 2013 · 0 comments
Open

SAM headers should be changed to contain all chromosomes #163

mdelaurentis opened this issue Feb 13, 2013 · 0 comments

Comments

@mdelaurentis
Copy link
Contributor

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.

@ghost ghost assigned mdelaurentis Feb 13, 2013
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant