#!/bin/bash -Ceuo pipefail
set -euo pipefail
export JAVA_TOOL_OPTIONS='-Xmx8g -XX:+UseParallelGC'
PicardCommandLine MarkDuplicates \
    I='OP-XX1230_CC-WGSLIB01-gDNA-N.sorted.bam' \
    O='OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bam' \
    M='OP-XX1230_CC-WGSLIB01-gDNA-N.markdup_metrics.txt' \
    TMP_DIR="$PWD" \
    VALIDATION_STRINGENCY=SILENT \
    MAX_RECORDS_IN_RAM=250000 \
    REMOVE_DUPLICATES=false \
    CREATE_INDEX=true \
    TAG_DUPLICATE_SET_MEMBERS=true \
    TAGGING_POLICY=All \
    OPTICAL_DUPLICATE_PIXEL_DISTANCE=2500 \
    READ_NAME_REGEX='^(?:[^:]*:){4}([^:]+):([^:]+):([^:]+)(?::[^:]*)'
# Picard/htsjdk may write <rg>.markdup.bai; we standardize on <rg>.markdup.bam.bai for Nextflow outputs
if [[ -f 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bam.bai' ]]; then
    :
elif [[ -f 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bai' ]]; then
    ln -sf 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bai' 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bam.bai'
else
    echo "MarkDuplicates did not create 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bam.bai' or 'OP-XX1230_CC-WGSLIB01-gDNA-N.markdup.bai'" >&2
    exit 1
fi