#!/bin/bash -Ceuo pipefail
if [ "1" -ne 10 ]; then
    printf '[manta]\nminSomaticScore = 1\n' > manta_override.ini
fi

configManta.py \
    --tumorBam   FFPE_ref_01_gDNA_pool_01_A23LGTTLT4_1.recalibrated.bam \
    --normalBam  FFPE_GM24385_pool_gDNA_01_B23LG35LT4_1.recalibrated.bam \
    --reference  Homo_sapiens_assembly38.fasta \
    --runDir     manta \
    --fusionGeneBed foresight_panel_for_manta.bed \
    --maxFusionBinPartners 20 \
    --config manta_override.ini \
    --exome

python manta/runWorkflow.py -m local -j 14

mv manta/results/variants/candidateSmallIndels.vcf.gz     seracare_100pct_rep1.manta_fusion.candidate_small_indels.vcf.gz
mv manta/results/variants/candidateSmallIndels.vcf.gz.tbi seracare_100pct_rep1.manta_fusion.candidate_small_indels.vcf.gz.tbi
mv manta/results/variants/candidateSV.vcf.gz              seracare_100pct_rep1.manta_fusion.candidate_sv.vcf.gz
mv manta/results/variants/candidateSV.vcf.gz.tbi          seracare_100pct_rep1.manta_fusion.candidate_sv.vcf.gz.tbi
mv manta/results/variants/somaticSV.vcf.gz                seracare_100pct_rep1.manta_fusion.somatic_sv.vcf.gz
mv manta/results/variants/somaticSV.vcf.gz.tbi            seracare_100pct_rep1.manta_fusion.somatic_sv.vcf.gz.tbi

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_VARIANT_CALLING_SOMATIC_ALL:BAM_VARIANT_CALLING_SOMATIC_MANTA_FUSION:MANTA_FUSION_CALL":
    manta: $( configManta.py --version )
END_VERSIONS