Make truncating the histogram optional in CollectInsertSizeMetrics #1779
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
|
@yfarjoun could you review this small PR? |
yfarjoun
suggested changes
Feb 27, 2022
There was a problem hiding this comment.
Thanks @takutosato for this enhancement.
A few requests to make the changes more inline with Picard "philosophy"....
| import htsjdk.samtools.SAMFileWriter; | ||
| import htsjdk.samtools.SAMFileWriterFactory; | ||
| import htsjdk.samtools.SAMRecordSetBuilder; | ||
| import htsjdk.samtools.*; |
There was a problem hiding this comment.
please bring back the individual imports
src/test/java/picard/analysis/CollectInsertSizeMetricsTest.java
Outdated
Show resolved
Hide resolved
| setBuilder.addPair("outlierRead", 0, start, start + outlierInsertSize - readLength, | ||
| false, false, cigar, cigar, false, true, 60); | ||
|
|
||
| final SAMFileWriter writer = new SAMFileWriterFactory().setCreateIndex(true).makeBAMWriter(setBuilder.getHeader(), false, testSamFile); |
There was a problem hiding this comment.
use try-with-resources here.
| } | ||
|
|
||
| // Add an outlier | ||
| int outlierInsertSize = 2000; |
| final int maxInsertSize = 400; | ||
| int i = 0; | ||
| int start = 1000; | ||
| String cigar = readLength + "M"; |
| final int minInsertSize = 100; | ||
| final int maxInsertSize = 400; | ||
| int i = 0; | ||
| int start = 1000; |
| final MetricsFile<InsertSizeMetrics, Comparable<?>> output = new MetricsFile<>(); | ||
| output.read(new FileReader(outfile)); | ||
|
|
||
| Assert.assertEquals(output.getAllHistograms().get(0).get(outlierInsertSize).getValue(), 1.0); |
There was a problem hiding this comment.
run with other argument to show that this argument makes a difference
| setBuilder.setReadLength(readLength); | ||
| final int minInsertSize = 100; | ||
| final int maxInsertSize = 400; | ||
| int i = 0; |
Co-authored-by: Yossi Farjoun <[email protected]>
Co-authored-by: Yossi Farjoun <[email protected]>
|
@takutosato are you able to complete the last of @yfarjoun 's requested changes for this PR? |
|
From what I can see this looks almost ready to go? @takutosato would you like me to take a look at the last few changes requested or do you have some time to do it? |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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.
Description
It's true that you can adjust some existing parameters to turn off truncating, but sometimes it's easier to just turn it off and the user should have that option.
Checklist (never delete this)
Never delete this, it is our record that procedure was followed. If you find that for whatever reason one of the checklist points doesn't apply to your PR, you can leave it unchecked but please add an explanation below.
Content
Review
For more detailed guidelines, see https://github.com/broadinstitute/picard/wiki/Guidelines-for-pull-requests