#!/bin/bash -euo pipefail
draw_fusions.R \
--fusions=HD789_0002_RNA_0008_23H5VFLT4_s19.FusionInspector.fusions.tsv \
--alignments=HD789_0002_RNA_0008_23H5VFLT4_s19.Aligned.sortedByCoord.out.bam \
--output=HD789_0002_RNA_0008_23H5VFLT4_s19_combined_fusions_arriba_visualisation.pdf \
--annotation=ref_annot.gtf \
--cytobands=cytobands_hg38_GRCh38_v2.4.0.tsv \
--proteinDomains=protein_domains_hg38_GRCh38_v2.4.0.gff3 \
cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:FUSIONINSPECTOR_WORKFLOW:ARRIBA_VISUALISATION":
arriba: $(arriba -h | grep 'Version:' 2>&1 | sed 's/Version: //')
END_VERSIONS