#!/bin/bash -Ceuo pipefail
merge_vcf_characteristics.py \
--target Sig_18_Blood.GERMLINE.merged.vcf.gz \
--source chip_pon.chip_pon.harmonized.vcf.gz \
--output Sig_18_Blood.CHIP_PON.merged.vcf.gz \
--prefix CHIP_PON \
tabix -p vcf Sig_18_Blood.CHIP_PON.merged.vcf.gz
cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_CHIP_DETECTION:VCF_MERGE_CHARACTERISTICS:MERGE_CHIP_PON":
merge_vcf_characteristics: $(merge_vcf_characteristics.py --version 2>&1 | grep -o 'version.*' | cut -d' ' -f2 || echo "1.0.0")
cyvcf2: $(python -c "import cyvcf2; print(cyvcf2.__version__)")
END_VERSIONS