#!/bin/bash -Ceuo pipefail
bcftools norm \
    -m-both \
    -a \
    -f Homo_sapiens_assembly38.fasta \
    --output-type v \
    1029_ILZ-T1-TDNA-1_B23VCMFLT4_1_vs_1029_ILZ-N1-BDNA-1_A23MHW3LT4_1.snv_indel.phased.vep.final.vcf.gz \
| bcftools norm \
    -m-both \
    --output-type z \
    --output 1029_ILZ-T1-TDNA-1_B23VCMFLT4_1_vs_1029_ILZ-N1-BDNA-1_A23MHW3LT4_1.norm.vcf.gz

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