#!/bin/bash -euo pipefail
junction_saturation.py \
    -i A549_FFPE_01_RNA_0001_B23LG7FLT4_3.Aligned.sortedByCoord.out.bam A549_FFPE_01_RNA_0001_B23LG7FLT4_3.Aligned.sortedByCoord.out.bam.bai \
    -r ref_annot.bed \
    -o A549_FFPE_01_RNA_0001_B23LG7FLT4_3 \


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