#!/bin/bash -euo pipefail
junction_annotation.py \
-i 1173_D6K_T1_RNA_SLD_01_A23T55JLT4_2.Aligned.sortedByCoord.out.bam 1173_D6K_T1_RNA_SLD_01_A23T55JLT4_2.Aligned.sortedByCoord.out.bam.bai \
-r ref_annot.bed \
-o 1173_D6K_T1_RNA_SLD_01_A23T55JLT4_2 \
\
2> >(grep -v 'E::idx_find_and_load' | tee 1173_D6K_T1_RNA_SLD_01_A23T55JLT4_2.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