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

Incorrect bin field set in bam file #144

Open
jrderuiter opened this issue Jun 19, 2017 · 0 comments
Open

Incorrect bin field set in bam file #144

jrderuiter opened this issue Jun 19, 2017 · 0 comments

Comments

@jrderuiter
Copy link

jrderuiter commented Jun 19, 2017

I'm running into issues when using picard tools to process bam files produced by ngs-disambiguate (which uses bamtools). The error suggests that the bin field is not being set correctly in the bam file that is being produced:

Exception in thread "main" htsjdk.samtools.SAMFormatException: SAM validation error: ERROR: Record 898366, Read name HWI-ST867:205:C2MWWACXX:5:1315:5951:11990, bin field of BAM record does not equal value computed based on alignment start and end, and length of sequence to which read is aligned

This error does not occur with the input bam files.

On biostars (https://www.biostars.org/p/110004/), it was suggested that this error might be due to bamtools not treating edge cases correctly. Is this indeed the case?

See AstraZeneca-NGS/disambiguate#6 for the original issue.

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