#!/bin/bash -euo pipefail
draw_fusions.R \
--fusions=1173_G2E_T1_RNA_SLD_01_A23T55JLT4_1.FusionInspector.fusions.tsv \
--alignments=1173_G2E_T1_RNA_SLD_01_A23T55JLT4_1.Aligned.sortedByCoord.out.bam \
--output=1173_G2E_T1_RNA_SLD_01_A23T55JLT4_1_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