#!/bin/bash -euo pipefail
vcf_collect.py \
--fusioninspector 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.FusionInspector.fusions.abridged.tsv.annotated.coding_effect \
--fusionreport 659_bjQ-T1-TRNA-1_B23MHV2LT4_1_fusionreport_index.html \
--fusioninspector_gtf 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.tsv \
--fusionreport_csv 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.fusions.csv \
--hgnc hgnc_complete_set.txt \
--sample 659_bjQ-T1-TRNA-1_B23MHV2LT4_1 \
--out 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.fusion_data.vcf \
--arriba_fusions 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.arriba.fusions.filtered_pass_only.tsv \
--use_chr_prefix

sort_fusion_vcf_mateid.py \
--input 659_bjQ-T1-TRNA-1_B23MHV2LT4_1.fusion_data.vcf

bgzip 659_bjQ-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