#!/bin/bash -euo pipefail
inner_distance.py \
-i tih_rna_sample_00090_23H5VFLT4_s32.Aligned.sortedByCoord.out.bam tih_rna_sample_00090_23H5VFLT4_s32.Aligned.sortedByCoord.out.bam.bai \
-r ref_annot.bed \
-o tih_rna_sample_00090_23H5VFLT4_s32 \
\
> stdout.txt
head -n 2 stdout.txt > tih_rna_sample_00090_23H5VFLT4_s32.inner_distance_mean.txt
cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:QC_WORKFLOW:BAM_RSEQC:RSEQC_INNERDISTANCE":
rseqc: $(inner_distance.py --version | sed -e "s/inner_distance.py //g")
END_VERSIONS