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

fix 4648 #4670

Merged
merged 1 commit into from
Apr 19, 2018
Merged

fix 4648 #4670

merged 1 commit into from
Apr 19, 2018

Conversation

SHuang-Broad
Copy link
Contributor

fix #4648

The problems are in clipping the alignments when de-overlapping the alignments, particularly computing the new ref and read spans.

Also adds a check on arguments used for constructing alignments, which helped finding errors in test data (now corrected)

@codecov-io
Copy link

codecov-io commented Apr 18, 2018

Codecov Report

Merging #4670 into master will increase coverage by 0.008%.
The diff coverage is 100%.

@@               Coverage Diff               @@
##              master     #4670       +/-   ##
===============================================
+ Coverage     79.835%   79.843%   +0.008%     
- Complexity     17328     17333        +5     
===============================================
  Files           1074      1074               
  Lines          62907     62921       +14     
  Branches       10181     10182        +1     
===============================================
+ Hits           50222     50238       +16     
  Misses          8704      8704               
+ Partials        3981      3979        -2
Impacted Files Coverage Δ Complexity Δ
.../discovery/alignment/ContigAlignmentsModifier.java 87.195% <100%> (+1.661%) 45 <0> (+4) ⬆️
...park/sv/discovery/alignment/AlignmentInterval.java 93.04% <100%> (+0.237%) 84 <1> (+1) ⬆️
...lugin/DefaultGATKReadFilterArgumentCollection.java 100% <0%> (ø) 4% <0%> (ø) ⬇️

Copy link
Member

@cwhelan cwhelan left a comment

Choose a reason for hiding this comment

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

This looks fine, thanks for fixing.

newTigStart = input.startInAssembledContig;
newTigEnd = input.endInAssembledContig - clipLengthOnRead;
newTigStart = originalContigStart;
newTigEnd = Math.min(originalContigEnd - clipLengthOnRead,
Copy link
Member

Choose a reason for hiding this comment

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

As discussed in person maybe add a couple of notes in comments to highlight why this is necessary?

* fixed bug in ref. span and read span when clipping alignment
* add argument check in AlignmentInterval
* correct invalid test data
@SHuang-Broad SHuang-Broad merged commit 3a8d9e3 into master Apr 19, 2018
@SHuang-Broad SHuang-Broad deleted the sh-sv-issues-4648 branch April 19, 2018 20:54
cwhelan pushed a commit to cwhelan/gatk-linked-reads that referenced this pull request May 25, 2018
* fixed bug in ref. span and read span when clipping alignment
* add argument check in AlignmentInterval
* correct invalid test data
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Error in StructuralVariationDiscoveryPipelineSpark: invalid interval in CpxVariantInterpreter
3 participants