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
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?
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:
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.
The text was updated successfully, but these errors were encountered: