Skip to content

Conversation

nh13
Copy link
Member

@nh13 nh13 commented Feb 20, 2025

…n base qualities are missing.

Fixes: #1031

This has the following improvements:

  1. A better error message is displayed, with the name of the read. Previously, a generic ArrayIndexOutOfBoundsException with on mention of which read.
  2. If the --missing-base-quality is not given, then the exception occurs (no change in current behavior)
  3. If the --missing-base-quality is given, then it is used (gracefully)

My only question is if we should have a non-None value by default. See #1031 for a discussion on why that might be a bad idea.

@@ -113,6 +113,7 @@ class CallDuplexConsensusReads
val maxReadsPerStrand: Option[Int] = None,
@arg(doc="The number of threads to use while consensus calling.") val threads: Int = 1,
@arg(doc="Consensus call overlapping bases in mapped paired end reads") val consensusCallOverlappingBases: Boolean = true,
@arg(doc = "The phred-scaled base quality to use when base qualities are missing in the input") val missingBaseQuality: Option[PhredScore] = None,
Copy link
Member Author

Choose a reason for hiding this comment

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

Suggested change
@arg(doc = "The phred-scaled base quality to use when base qualities are missing in the input") val missingBaseQuality: Option[PhredScore] = None,
@arg(doc="The phred-scaled base quality to use when base qualities are missing in the input") val missingBaseQuality: Option[PhredScore] = None,

Copy link

codecov bot commented Feb 20, 2025

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 95.68%. Comparing base (f0fd24c) to head (d7efcd2).

Additional details and impacted files
@@           Coverage Diff           @@
##             main    #1032   +/-   ##
=======================================
  Coverage   95.68%   95.68%           
=======================================
  Files         126      126           
  Lines        7410     7415    +5     
  Branches      500      524   +24     
=======================================
+ Hits         7090     7095    +5     
  Misses        320      320           
Flag Coverage Δ
unittests 95.68% <100.00%> (+<0.01%) ⬆️

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 assigned nh13 and unassigned tfenne Mar 5, 2025
@nh13 nh13 marked this pull request as draft March 5, 2025 18:23
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

How should we handle missing base qualities in the consensus caller?
2 participants