#!/bin/bash -Ceuo pipefail
bcftools norm \
-m-both \
-a \
-f Homo_sapiens_assembly38.fasta \
--output-type v \
1029_1BK-T1-TDNA-1_B23KFCWLT4_1_vs_1029_1BK-N1-BDNA-1_B23KGCJLT4_1.snv_indel.phased.vep.final.vcf.gz \
| bcftools norm \
-m-both \
--output-type z \
--output 1029_1BK-T1-TDNA-1_B23KFCWLT4_1_vs_1029_1BK-N1-BDNA-1_B23KGCJLT4_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