Skip to content

Fix the iteration over the CIGAR segments#1202

Merged
daviesrob merged 1 commit intosamtools:developfrom
valeriuo:overlap_fix
Jan 15, 2021
Merged

Fix the iteration over the CIGAR segments#1202
daviesrob merged 1 commit intosamtools:developfrom
valeriuo:overlap_fix

Conversation

@valeriuo
Copy link
Copy Markdown
Contributor

When encountering a new CIGAR segment, the index icig, which indicates the position of a base inside the current segment, has to be reset. This fix changes the reset value to -1 (instead of 0), so that the next first matching base starts from icig=0 in the calling method (tweak_overlap_quality).

Replaces #1196
Fixes #1195

… indels, clips and skips, so that the next MATCH starts from 0 index.

Co-authored-by: Valeriu Ohan <vo2@sanger.ac.uk>
@daviesrob
Copy link
Copy Markdown
Member

Looks good, although the operation of tweak_overlap_quality() is certainly somewhat obscure.

@daviesrob daviesrob merged commit 481bd22 into samtools:develop Jan 15, 2021
@valeriuo valeriuo deleted the overlap_fix branch January 18, 2021 10:34
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

Successfully merging this pull request may close these issues.

during tweak quality of overlapped paired of reads, the quality of some bases in the end of overlap region are not tweaked

2 participants