-
Notifications
You must be signed in to change notification settings - Fork 308
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
Add quality score binner #1485
Add quality score binner #1485
Conversation
Test FAILed. Build result: FAILURE[...truncated 16 lines...] > /home/jenkins/git2/bin/git rev-parse origin/pr/1485/merge^{commit} # timeout=10 > /home/jenkins/git2/bin/git branch -a --contains d0fb3edd7021762f878ef7a5ab66ab37dd25c38e # timeout=10 > /home/jenkins/git2/bin/git rev-parse remotes/origin/pr/1485/merge^{commit} # timeout=10Checking out Revision d0fb3edd7021762f878ef7a5ab66ab37dd25c38e (origin/pr/1485/merge) > /home/jenkins/git2/bin/git config core.sparsecheckout # timeout=10 > /home/jenkins/git2/bin/git checkout -f d0fb3edd7021762f878ef7a5ab66ab37dd25c38eFirst time build. Skipping changelog.Triggering ADAM-prb ? 2.3.0,2.11,1.6.1,centosTriggering ADAM-prb ? 2.6.0,2.10,1.6.1,centosTriggering ADAM-prb ? 2.3.0,2.11,2.0.0,centosTriggering ADAM-prb ? 2.3.0,2.10,2.0.0,centosTriggering ADAM-prb ? 2.6.0,2.11,2.0.0,centosTriggering ADAM-prb ? 2.6.0,2.10,2.0.0,centosTriggering ADAM-prb ? 2.3.0,2.10,1.6.1,centosTriggering ADAM-prb ? 2.6.0,2.11,1.6.1,centosADAM-prb ? 2.3.0,2.11,1.6.1,centos completed with result SUCCESSADAM-prb ? 2.6.0,2.10,1.6.1,centos completed with result FAILUREADAM-prb ? 2.3.0,2.11,2.0.0,centos completed with result SUCCESSADAM-prb ? 2.3.0,2.10,2.0.0,centos completed with result FAILUREADAM-prb ? 2.6.0,2.11,2.0.0,centos completed with result SUCCESSADAM-prb ? 2.6.0,2.10,2.0.0,centos completed with result FAILUREADAM-prb ? 2.3.0,2.10,1.6.1,centos completed with result FAILUREADAM-prb ? 2.6.0,2.11,1.6.1,centos completed with result SUCCESSNotifying endpoint 'HTTP:https://webhooks.gitter.im/e/ac8bb6e9f53357bc8aa8'Test FAILed. |
Jenkins, retest this please. |
Test PASSed. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good! Nice test coverage. Does it help performance much?
@@ -115,6 +115,8 @@ class TransformArgs extends Args4jBase with ADAMSaveAnyArgs with ParquetArgs { | |||
var mdTagsFragmentSize: Long = 1000000L | |||
@Args4jOption(required = false, name = "-md_tag_overwrite", usage = "When adding MD tags to reads, overwrite existing incorrect tags.") | |||
var mdTagsOverwrite: Boolean = false | |||
@Args4jOption(required = false, name = "-bin_quality_scores", usage = "Rewrites quality scores of reads into bins. User must provide bin description string.") |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
How about
Rewrites quality scores of reads into bins from a string of bin descriptions, e.g. 0,20,10;20,40,30.
// if we have multiple bins, validate them | ||
// - check that we don't have gaps between bins | ||
// - check that we don't have overlapping bins | ||
if (bins.size > 1) { |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
would this check be better in QualityScoreBin.apply
?
I'd be ok merging this with the doc string fix suggested above |
Resolves bigdatagenomics#1462. Adds an API for binning quality scores to `AlignmentRecordRDD` and `FragmentRDD`. Exposes this on the `Transform` CLI.
2d9cb2c
to
a089670
Compare
@heuermh addressed the doc issue. This is good to go from my side. |
Test PASSed. |
Thank you, @fnothaft |
…formFragments` Resolves bigdatagenomics#1359. Also, propegates the quality score binner (bigdatagenomics#1485) up to the `transformFragments` CLI.
…formFragments` Resolves bigdatagenomics#1359. Also, propegates the quality score binner (bigdatagenomics#1485) up to the `transformFragments` CLI.
Resolves #1462. WIP. Still need to: