#!/bin/bash -euo pipefail
bam_stat.py \
    -i 659_8h-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam 659_8h-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam.bai \
     \
    > 659_8h-T1-TRNA-1_B23MHV2LT4_1.bam_stat.txt

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