#!/bin/bash -Ceuo pipefail
sentieon driver \
    -t 12 \
    -r Homo_sapiens_assembly38.fasta \
    --interval_padding 100 \
    -i HCC1395_BL.recalibrated.bam \
    --interval xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_chr21.bed \
    --algo TNhaplotyper2 \
    --tumor_sample HCC1395_BL \
    --call_germline_sites \
    --germline_vcf af-only-gnomad.hg38.vcf.gz \
    --pon pon.hotspot_protected.raw.vcf.gz \
    HCC1395_BL.tnseq_tumor_only.vcf.gz \
    --algo OrientationBias --tumor_sample HCC1395_BL HCC1395_BL.tnseq_tumor_only.orientation_data.tsv \
    --algo ContaminationModel --tumor_sample HCC1395_BL --vcf af-only-gnomad.hg38.vcf.gz --tumor_segments HCC1395_BL.tnseq_tumor_only.segments HCC1395_BL.tnseq_tumor_only.contamination_data.tsv

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_CHIP_DETECTION:BAM_VARIANT_CALLING_TUMOR_ONLY_TNSEQ:SENTIEON_TNHAPLOTYPER2":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS