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