#!/bin/bash -Ceuo pipefail
/opt/deepvariant/bin/run_deepvariant \
    --ref=Homo_sapiens_assembly38.fasta \
    --reads=HCC1395_BL.recalibrated.bam \
    --output_vcf=germline_control_1.dv_rescue.vcf.gz \
    --output_gvcf=germline_control_1.dv_rescue.g.vcf.gz \
    --model_type=WES \
    --regions=germline_rescue_targets.bed \
     \
    --intermediate_results_dir=tmp \
    --num_shards=4

cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_GERMLINE_VC:BAM_VARIANT_CALLING_SENTIEON_HAPLOTYPER_RF:DV_RESCUE":
    deepvariant: 1.6.1
END_VERSIONS