#!/bin/bash -Ceuo pipefail
sentieon driver \
    -t 12 \
    -r Homo_sapiens_assembly38.fasta \
     \
    -i reference-NA12878-downsampled-500k.bam -i HCC1395_tumor.recalibrated.bam \
    --interval xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity_padded100bp.bed \
    --algo TNhaplotyper2 \
    --tumor_sample HCC1395_tumor --normal_sample reference-NA12878 \
    --call_germline_sites \
    --germline_vcf af-only-gnomad.hg38.vcf.gz \
    --pon pon_tnseq_42_curated_v4.vcf.gz \
    positive_somatic_control_1.tnseq.vcf.gz \
    --algo OrientationBias --tumor_sample HCC1395_tumor positive_somatic_control_1.tnseq.orientation_data.tsv \
    --algo ContaminationModel --tumor_sample HCC1395_tumor --normal_sample reference-NA12878 --vcf af-only-gnomad.hg38.vcf.gz --tumor_segments positive_somatic_control_1.tnseq.segments positive_somatic_control_1.tnseq.contamination_data.tsv

cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_SOMATIC_VC:BAM_VARIANT_CALLING_SOMATIC_TNSEQ:SENTIEON_TNHAPLOTYPER2":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS