#!/bin/bash -Ceuo pipefail
gt_correction.py \
    --input-vcf germline_control_2.haplotyper_rf.rf_filtered.vcf.gz \
    --output-vcf germline_control_2.haplotyper_rf.gt_corrected.vcf.gz \
    --gt-model gt_correction_model_v2.joblib \
    --repeatmasker-bed rf_repeatmasker.bed.gz \
    --threshold 0.5 \


tabix -p vcf germline_control_2.haplotyper_rf.gt_corrected.vcf.gz

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