#!/bin/bash -Ceuo pipefail
set -euo pipefail
nipt_emit_alpha_features.py \
    --vcf HG02855_x_HG02886_FF_5.vcf \
    --bam HG02855_x_HG02886_FF_5.mix.bam \
    --ref Homo_sapiens_assembly38.fasta \
    --out-tsv HG02855_x_HG02886_FF_5.features.tsv \
    --out-vcf HG02855_x_HG02886_FF_5.augmented.vcf \
    --target-tsv HG02855_x_HG02886_FF_5.targets.tsv \
    --indels-only \
    --ff 0.05 \
    --alphas 0.012,0.05,0.10 \
    --window-pad 250 \
    --min-mapq 20 \
    --progress 100 2>&1 | tee HG02855_x_HG02886_FF_5.log