#!/bin/bash -Ceuo pipefail
diff_dv_sentieon.py \
    --dv-vcf germline_control_1.dv_rescue.vcf.gz \
    --sentieon-vcf germline_control_1.haplotyper.unfiltered.vcf.gz \
    --output-vcf germline_control_1.dv_only_sites.vcf.gz

tabix -p vcf germline_control_1.dv_only_sites.vcf.gz

cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_GERMLINE_VC:BAM_VARIANT_CALLING_SENTIEON_HAPLOTYPER_RF:DIFF_DV_SENTIEON":
    python: $(python --version | sed 's/Python //g')
    cyvcf2: $(python -c "import cyvcf2; print(cyvcf2.__version__)")
END_VERSIONS