deepTools
deepTools copied to clipboard
--filterRNAstrand usage for alignmentSieve is missleading
--filterRNAstrand documentation is misleading.
The actual filtration, from what I've tested, actually works as documented for bamCoverage --filterRNAstrand (see bamCoverage --filterRNAstrand.
It indicates this: Selects RNA-seq reads (single-end or paired-end) in the given strand. (Default: %(default)s),
but works like this: Selects RNA-seq reads (single-end or paired-end) originating from genes on the given strand. This option assumes a standard dUTP-based library preparation (that is, --filterRNAstrand=forward keeps minus-strand reads, which originally came from genes on the forward strand using a dUTP-based method)