#!/bin/bash -Ceuo pipefail
tabix custom_Sig_18_tumor_normal.tnseq.whatshap.chr21.phased.vcf.gz

bcftools concat \
    --output custom_Sig_18_tumor_normal.tnseq.whatshap.phased.vcf.gz \
    --allow-overlaps --output-type z --write-index=tbi \
    --threads 6 \
    custom_Sig_18_tumor_normal.tnseq.whatshap.chr21.phased.vcf.gz

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_VARIANT_CALLING_SOMATIC_ALL:CONCAT_WHATSHAP_TNSEQ_VCFS":
    bcftools: $(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*$//')
END_VERSIONS