#!/bin/bash -euo pipefail
read_distribution.py \
    -i BVT_FFPE_TRNA_bst_01_02_A23WKFTLT4_2.Aligned.sortedByCoord.out.bam BVT_FFPE_TRNA_bst_01_02_A23WKFTLT4_2.Aligned.sortedByCoord.out.bam.bai \
    -r ref_annot.bed \
    > BVT_FFPE_TRNA_bst_01_02_A23WKFTLT4_2.read_distribution.txt

cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:QC_WORKFLOW:BAM_RSEQC:RSEQC_READDISTRIBUTION":
    rseqc: $(read_distribution.py --version | sed -e "s/read_distribution.py //g")
END_VERSIONS