#!/bin/bash -euo pipefail
junction_annotation.py \
    -i GM24385_0001_RNA_0001_23H5VFLT4_s48.Aligned.sortedByCoord.out.bam GM24385_0001_RNA_0001_23H5VFLT4_s48.Aligned.sortedByCoord.out.bam.bai \
    -r ref_annot.bed \
    -o GM24385_0001_RNA_0001_23H5VFLT4_s48 \
     \
    2> >(grep -v 'E::idx_find_and_load' | tee GM24385_0001_RNA_0001_23H5VFLT4_s48.junction_annotation.log >&2)

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