Fix decoding of long CIGAR ops in BAM #173
Merged
+46
−5
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Problem
cljam.io.bam
fails to decode alignments with CIGAR operations longer than 134217728Cause
A CIGAR is stored as an array of
op_len<<4|op
inuint32_t[n_cigar_op]
.Thus
op_len
can be as long as 28bit0x0FFFFFFF
but it is not handled properly due to type promotion from int to long.Changes
Use
Integer/toUnsignedLong
to get a positive long value.Tests
lein check
🆗lein test :all
🆗lein cloverage
🆗lein eastwood
🆗