#!/bin/bash -Ceuo pipefail
analyze \
    --tumor-bam HCC1395_BL.recalibrated.bam \
    --normal-bam normal/reference-NA12878-downsampled-500k.bam \
    --case-id negative_somatic_control_1 \
    --save \
    --save-location . \
    --save-format both \
    --name negative_somatic_control_1 \
    --cores 16 \
    --coverage 100 \
    --confidence-interval 0.95 \
     --microsatellites-list 1500_dropped_panel_with_boosted_msi_regions.tsv --model mi_msi_v0_4_0_200x_attn.model --use-attention 

mimsi_results_to_json.py negative_somatic_control_1_results.txt negative_somatic_control_1_results.json

cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_SOMATIC_VC:MIMSI_ANALYZE":
    mimsi: $(analyze --version 2>&1 | sed 's/.*version - //; s/ .*//' || echo 'unknown')
END_VERSIONS