command <<< | command <<< | ||||
set -e -o pipefail | set -e -o pipefail | ||||
path=$(dirname ${bam}) | |||||
path=$(dirname ${sep=' ' bam}) | |||||
mv ${sep=' ' bai} $path | mv ${sep=' ' bai} $path | ||||
bedtools multicov -bams ${sep=' ' bam} -bed ${bed_file} > revread_count.txt | bedtools multicov -bams ${sep=' ' bam} -bed ${bed_file} > revread_count.txt | ||||