Files
rnaseq/src/bedtools_genomecov/script.sh
CI 94b5f2024b Build branch multiple_fixes with version multiple_fixes (e14b535)
Build pipeline: viash-hub.rnaseq.multiple-fixes-tnlvp

Source commit: e14b535c3b

Source message: update arguments
2024-09-19 09:48:31 +00:00

26 lines
569 B
Bash

#!/bin/bash
set -eo pipefail
prefix_forward="forward"
prefix_reverse="reverse"
if [ $par_strandedness == 'reverse' ]; then
prefix_forward="reverse"
prefix_reverse="forward"
fi
bedtools genomecov \
-ibam $par_bam \
-bg \
-strand + \
$par_extra_bedtools_args | bedtools sort > $prefix_forward.bedGraph
bedtools genomecov \
-ibam $par_bam \
-bg \
-strand - \
$par_extra_bedtools_args | bedtools sort > $prefix_reverse.bedGraph
mv $prefix_forward.bedGraph $par_bedgraph_forward
mv $prefix_reverse.bedGraph $par_bedgraph_reverse