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