-
Notifications
You must be signed in to change notification settings - Fork 589
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
handle unmarked duplicates with mate MQ = 0 in Mutect2 #5734
Conversation
Codecov Report
@@ Coverage Diff @@
## master #5734 +/- ##
==============================================
+ Coverage 87.061% 87.071% +0.01%
- Complexity 31880 31902 +22
==============================================
Files 1940 1940
Lines 146804 146832 +28
Branches 16234 16238 +4
==============================================
+ Hits 127809 127848 +39
+ Misses 13073 13064 -9
+ Partials 5922 5920 -2
|
@davidbenjamin Results on FC looked really good. Reviewed manually. |
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.
Minor change and one question. I will approve on the assumption that you can address without my need to re-review.
assemblyRegion.removeAll(readStubs); | ||
} | ||
|
||
private void removeUnmarkedDuplicates(final AssemblyRegion assemblyRegion) { |
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.
This looks like it considers:
- SVs: This code checks whether the mates are clustered, even if on a different contig from the supporting reads of the variant.
- Duplicate start and end positions of supporting reads.
Please tell me if I have misinterpreted. If I have not misinterpreted, then no action required.
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.
That's right.
@@ -1150,4 +1150,20 @@ private static String keyForVariant( final VariantContext variant ) { | |||
return String.format("%s:%d-%d %s, %s", variant.getContig(), variant.getStart(), variant.getEnd(), variant.getReference(), | |||
variant.getAlternateAlleles().stream().map(Allele::getDisplayString).sorted().collect(Collectors.toList())); | |||
} | |||
|
|||
@Test | |||
public void testSachet() throws Exception { |
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.
Can we rename to something about missed duplicates? And if you want to thank Sachet in a comment, that should be okay.
@davidbenjamin Back to you |
@LeeTL1220 This is for Sachet's artifact, ready to review pending the results on Firecloud looking good.