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

unmapped.bam lacks BAM flag 0x8 for read pairs where both mates are unmapped #3

Open
cbrueffer opened this issue Apr 9, 2015 · 0 comments

Comments

@cbrueffer
Copy link
Contributor

TopHat up to version 2.0.14 doesn't set the 0x8 (mate unmapped) flag for read pairs where both mates are unmapped. This causes downstream software, e.g., samtools merge to error on the data.

The bug can be reproduced by running TopHat on the test_data set from the TopHat tutorial (http://ccb.jhu.edu/software/tophat/downloads/test_data.tar.gz)

samtools view -f 0x8 unmapped.bam returns no output, where it should yield 16 reads with this bit set.

Running the tophat-recondition script from https://github.com/cbrueffer/tophat-recondition restores the missing flags.

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