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

Do not trim reads when both ends are clipped in consensus calling #1026

Open
wants to merge 3 commits into
base: main
Choose a base branch
from

Conversation

nh13
Copy link
Member

@nh13 nh13 commented Feb 6, 2025

No description provided.

Copy link

codecov bot commented Feb 6, 2025

Codecov Report

Attention: Patch coverage is 80.00000% with 3 lines in your changes missing coverage. Please review.

Project coverage is 95.63%. Comparing base (f93fdfb) to head (4423d72).

Files with missing lines Patch % Lines
...ala/com/fulcrumgenomics/bam/SamRecordClipper.scala 80.00% 3 Missing ⚠️
Additional details and impacted files
@@            Coverage Diff             @@
##             main    #1026      +/-   ##
==========================================
- Coverage   95.65%   95.63%   -0.02%     
==========================================
  Files         126      126              
  Lines        7403     7405       +2     
  Branches      521      499      -22     
==========================================
+ Hits         7081     7082       +1     
- Misses        322      323       +1     
Flag Coverage Δ
unittests 95.63% <80.00%> (-0.02%) ⬇️

Flags with carried forward coverage won't be shown. Click here to find out more.

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@nh13 nh13 linked an issue Feb 6, 2025 that may be closed by this pull request
@nh13 nh13 changed the title test: show behavior of issue #1106 Do not trim reads when both ends are clipped in consensus calling Feb 6, 2025
@nh13 nh13 marked this pull request as ready for review February 6, 2025 21:29
@nh13 nh13 requested a review from tfenne as a code owner February 6, 2025 21:29
case true if rec.end >= mateEnd =>
case true if rec.end >= mateEnd =>
// account for any soft-clipping on the same side of the template
val trailingClippedBases = if (rec.end == mateEnd) Math.min(rec.cigar.trailingClippedBases, rec.mateCigar.get.trailingClippedBases) else 0
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this is wrong, mate unclipped end should be used...

@nh13 nh13 marked this pull request as draft February 6, 2025 21:53
@nh13 nh13 temporarily deployed to github-actions February 6, 2025 23:59 — with GitHub Actions Inactive
@nh13 nh13 temporarily deployed to github-actions February 6, 2025 23:59 — with GitHub Actions Inactive
@nh13 nh13 temporarily deployed to github-actions February 7, 2025 00:01 — with GitHub Actions Inactive
@nh13 nh13 temporarily deployed to github-actions February 7, 2025 00:01 — with GitHub Actions Inactive
@nh13 nh13 marked this pull request as ready for review February 7, 2025 00:01
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.

Soft clip disappeared after CallMolecularConsensusReads
1 participant