samtools / htsjdk

A Java API for high-throughput sequencing data (HTS) formats.
http://samtools.github.io/htsjdk/
283 stars 242 forks source link

Mf get read position from ref #1475

Closed fleharty closed 4 years ago

fleharty commented 4 years ago

Description

Adding functionality to optionally include softclips in getAlignmentBlocks

Things to think about before submitting:

tfenne commented 4 years ago

@fleharty I really don't like this - it adds a bunch more public methods to SAMRecord for something that is rather questionable. If you really want to treat soft clips as Ms I feel like the more correct thing to do would be to update the cigar on the record, and then use the existing getReadPositionAtRefPosition() method(s).

codecov-io commented 4 years ago

Codecov Report

Merging #1475 into master will increase coverage by 0.005%. The diff coverage is 78.571%.

@@               Coverage Diff               @@
##              master     #1475       +/-   ##
===============================================
+ Coverage     69.203%   69.207%   +0.005%     
- Complexity      8701      8709        +8     
===============================================
  Files            588       588               
  Lines          34581     34596       +15     
  Branches        5779      5785        +6     
===============================================
+ Hits           23931     23943       +12     
- Misses          8368      8369        +1     
- Partials        2282      2284        +2     
Impacted Files Coverage Δ Complexity Δ
src/main/java/htsjdk/samtools/SAMUtils.java 62.972% <77.778%> (+0.170%) 136.000 <7.000> (+5.000)
src/main/java/htsjdk/samtools/SAMRecord.java 68.148% <80.000%> (-0.175%) 299.000 <5.000> (+2.000) :arrow_down:
...samtools/util/AsyncBlockCompressedInputStream.java 76.000% <0.000%> (+4.000%) 13.000% <0.000%> (+1.000%)
fleharty commented 4 years ago

@tfenne That's an interesting idea. I'll try to accomplish this by modifying the CIGAR. Thanks.

fleharty commented 4 years ago

Closing this PR, I believe there is a better solution.