#!/bin/bash -euo pipefail
vcf_collect.py \
--fusioninspector 659_1q-T1-TRNA-1_B23MHV2LT4_1.FusionInspector.fusions.abridged.tsv.annotated.coding_effect \
--fusionreport 659_1q-T1-TRNA-1_B23MHV2LT4_1_fusionreport_index.html \
--fusioninspector_gtf 659_1q-T1-TRNA-1_B23MHV2LT4_1.tsv \
--fusionreport_csv 659_1q-T1-TRNA-1_B23MHV2LT4_1.fusions.csv \
--hgnc hgnc_complete_set.txt \
--sample 659_1q-T1-TRNA-1_B23MHV2LT4_1 \
--out 659_1q-T1-TRNA-1_B23MHV2LT4_1.fusion_data.vcf \
--arriba_fusions 659_1q-T1-TRNA-1_B23MHV2LT4_1.arriba.fusions.filtered_pass_only.tsv \
--use_chr_prefix
sort_fusion_vcf_mateid.py \
--input 659_1q-T1-TRNA-1_B23MHV2LT4_1.fusion_data.vcf
bgzip 659_1q-T1-TRNA-1_B23MHV2LT4_1.fusion_data.vcf
cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:FUSIONINSPECTOR_WORKFLOW:VCF_COLLECT":
python: $(python --version | sed 's/Python //g')
HGNC DB retrieval: $(cat HGNC-DB-timestamp.txt)
END_VERSIONS