Version 2.2.0
=== Notable changes in this release === 1. coverageBed will optionally only count features in BED file A (e.g. sequencing reads) that overlap with the intervals/windows in BED file B. This has been requested several times recently and facilitates CHiP-Seq and RNA-Seq experiments. 2. intersectBed can now require a minimum __reciprocal__ overlap between intervals in BED A and BED B. For example, previously, if one used -f 0.90, it required that a feature in B overlap 90% of the feature in A for the "hit" to be reported. If one adds the -r (reciprocal) option, the hit must also cover 90% of the feature in B. This helps to exclude overlaps between say small features in A and large features in B: A ========== B ********************************************************** -f 0.50 (Reported), whereas -f 0.50 -r (Not reported) 3. The score field has been changed to be a string. While this deviates from the UCSC definition, it allows one to track much more meaningful information about a feature/interval. For example, score could now be: 7.31E-05 (a p-value) 0.334577 (mean enrichment) 2:2.2:40:2 (several values encoded in a string) 4. closestBed now, by default, reports __all__ intervals in B that overlap equally with an interval in A. Previously, it merely reported the first such feature that appeared in B. Here's a cartoon explaining the difference. **Prior behavior** A ============== B.1 ++++++++++++++ B.2 ++++++++++++++ B.3 +++++++++ ----------------------------------------- Result = B.1 ++++++++++++++ **Current behavior** A ============== B.1 ++++++++++++++ B.2 ++++++++++++++ B.3 +++++++++ ----------------------------------------- Result = B.1 ++++++++++++++ B.2 ++++++++++++++ Using the -t option, one can also choose to report either the first or the last entry in B in the event of a tie. 5. Several other minor changes to the algorithms have been made to increase speed a bit.
Showing
- README 3 additions, 3 deletionsREADME
- RELEASE_HISTORY 54 additions, 0 deletionsRELEASE_HISTORY
- src/closestBed/closestBed.cpp 43 additions, 13 deletionssrc/closestBed/closestBed.cpp
- src/closestBed/closestBed.h 2 additions, 1 deletionsrc/closestBed/closestBed.h
- src/closestBed/closestMain.cpp 23 additions, 3 deletionssrc/closestBed/closestMain.cpp
- src/coverageBed/computeDistance.py 20 additions, 0 deletionssrc/coverageBed/computeDistance.py
- src/coverageBed/coverageBed.cpp 5 additions, 3 deletionssrc/coverageBed/coverageBed.cpp
- src/coverageBed/coverageBed.h 3 additions, 1 deletionsrc/coverageBed/coverageBed.h
- src/coverageBed/coverageMain.cpp 6 additions, 1 deletionsrc/coverageBed/coverageMain.cpp
- src/genomeCoverageBed/genomeCoverageMain.cpp 1 addition, 1 deletionsrc/genomeCoverageBed/genomeCoverageMain.cpp
- src/intersectBed/a.bed 1 addition, 0 deletionssrc/intersectBed/a.bed
- src/intersectBed/b.bed 1 addition, 0 deletionssrc/intersectBed/b.bed
- src/intersectBed/intersectBed.cpp 46 additions, 10 deletionssrc/intersectBed/intersectBed.cpp
- src/intersectBed/intersectBed.h 2 additions, 1 deletionsrc/intersectBed/intersectBed.h
- src/intersectBed/intersectMain.cpp 11 additions, 1 deletionsrc/intersectBed/intersectMain.cpp
- src/utils/bedFile/bedFile.cpp 19 additions, 10 deletionssrc/utils/bedFile/bedFile.cpp
- src/utils/bedFile/bedFile.h 1 addition, 1 deletionsrc/utils/bedFile/bedFile.h
- src/utils/version/version.h 1 addition, 1 deletionsrc/utils/version/version.h
Loading
Please register or sign in to comment