#!/bin/bash -Ceuo pipefail
deconflict_variants.py \
    --chip-vcf Sig_18_Blood.scored.vcf.gz --output-chip Sig_18_Blood.deconflicted_chip.vcf.gz \
    --germline-vcf Sig_18_Blood.germline.harmonized.vcf.gz \
    --somatic-vcf custom_Sig_18_tumor_normal.snv_indel.phased.vep.final.vcf.gz \
    --output-germline Sig_18_Blood.deconflicted_germline.vcf.gz \
    --output-somatic Sig_18_Blood.deconflicted_somatic.vcf.gz \
    --audit Sig_18_Blood.deconfliction_audit.tsv \
    --cosmic-heme cosmic_heme.tsv.gz \


# Chip output only exists when chip VCF was provided
if [ -f Sig_18_Blood.deconflicted_chip.vcf.gz ]; then
    tabix -p vcf Sig_18_Blood.deconflicted_chip.vcf.gz
fi
if [ -s Sig_18_Blood.deconflicted_germline.vcf.gz ]; then
    tabix -p vcf Sig_18_Blood.deconflicted_germline.vcf.gz
else
    touch Sig_18_Blood.deconflicted_germline.vcf.gz.tbi
fi

# Somatic output may be empty if no somatic VCF was provided
if [ -s Sig_18_Blood.deconflicted_somatic.vcf.gz ]; then
    tabix -p vcf Sig_18_Blood.deconflicted_somatic.vcf.gz
else
    touch Sig_18_Blood.deconflicted_somatic.vcf.gz.tbi
fi

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:VCF_DECONFLICT_VARIANTS:DECONFLICT_VARIANTS":
    deconflict_variants: $(deconflict_variants.py --version 2>&1 || echo "1.0.0")
    cyvcf2: $(python -c "import cyvcf2; print(cyvcf2.__version__)")
END_VERSIONS