Fixed logic relating to "split" coverage.
1. Fixed bug in BamAncillary::getBamBlocks where the switch statement was misbehaving. Changed to if/else 2. Added countSplitHits to bedFile class to proper count coverage from "blocked" alignments. Previously, coverage was counted separately for __each__ block yet what we actually want is the base coverage to be counted for each block, but we only want the blocks to be treated as __one__ alignment or feature. 3. Updated coverageBed to use getBamBlocks.
Showing
- src/coverageBed/coverageBed.cpp 6 additions, 10 deletionssrc/coverageBed/coverageBed.cpp
- src/genomeCoverageBed/genomeCoverageBed.cpp 4 additions, 1 deletionsrc/genomeCoverageBed/genomeCoverageBed.cpp
- src/intersectBed/intersectBed.cpp 1 addition, 0 deletionssrc/intersectBed/intersectBed.cpp
- src/utils/BamTools/BamAncillary.cpp 19 additions, 27 deletionssrc/utils/BamTools/BamAncillary.cpp
- src/utils/bedFile/bedFile.cpp 56 additions, 0 deletionssrc/utils/bedFile/bedFile.cpp
- src/utils/bedFile/bedFile.h 8 additions, 0 deletionssrc/utils/bedFile/bedFile.h
Loading
Please register or sign in to comment