#!/bin/bash -Ceuo pipefail
annotateFusionVcf.py \
--vcf COLO829_tumor_100pct_rep1.manta_fusion.somatic_sv.vcf.gz \
--genes protein_coding_genes.bed \
--actionable-bed foresight_panel_for_manta.bed \
--out COLO829_tumor_100pct_rep1.manta_fusion.somatic_sv.annotated.vcf
bgzip COLO829_tumor_100pct_rep1.manta_fusion.somatic_sv.annotated.vcf
tabix -p vcf COLO829_tumor_100pct_rep1.manta_fusion.somatic_sv.annotated.vcf.gz
cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_VARIANT_CALLING_SOMATIC_ALL:BAM_VARIANT_CALLING_SOMATIC_MANTA_FUSION:MANTA_FUSION_ANNOTATE":
python: $( python3 --version 2>&1 | sed 's/Python //' )
annotateFusionVcf: "1.0.0"
END_VERSIONS