#!/bin/bash -euo pipefail
arriba \
-x aih-tih-sc-b845f7-R1_B23WHYVLT4_1.Aligned.sortedByCoord.out.bam \
-a ref_genome.fa \
-g ref_annot.gtf \
-o aih-tih-sc-b845f7-R1_B23WHYVLT4_1.arriba.fusions.tsv \
-O aih-tih-sc-b845f7-R1_B23WHYVLT4_1.arriba.fusions.discarded.tsv \
-b blacklist_hg38_GRCh38_v2.4.0.tsv.gz \
-k known_fusions_hg38_GRCh38_v2.4.0.tsv.gz \
\
\
-p protein_domains_hg38_GRCh38_v2.4.0.gff3 \
cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:ARRIBA_WORKFLOW:ARRIBA":
arriba: $(arriba -h | grep 'Version:' 2>&1 | sed 's/Version: //')
END_VERSIONS