#!/bin/bash -Ceuo pipefail
INDEX=`find -L ./ -name "*.amb" | sed 's/\.amb$//'`

export bwt_max_mem=112G

sentieon bwa mem \
    -K 100000000 -Y -R '@RG\tID:B23KGCJLT4.1029_0I-N1-BDNA-1_B23KGCJLT4_1.L004\tPU:L004\tSM:1029_0I-N1-BDNA-1_B23KGCJLT4_1\tLB:1029_0I-N1-BDNA-1_B23KGCJLT4_1\tDS:s3://natera-platform-sandbox/pipeline-resources/ngi-igenomes/igenomes/Homo_sapiens/GATK/GRCh38/Sequence/WholeGenomeFasta/Homo_sapiens_assembly38.fasta\tPL:ILLUMINA' \
    -t 31 \
    $INDEX \
    1029_0I-N1-BDNA-1_B23KGCJLT4_1_R1.fastq.gz 1029_0I-N1-BDNA-1_B23KGCJLT4_1_R2.fastq.gz \
    | sentieon util sort \
        -r Homo_sapiens_assembly38.fasta \
        -t 31 \
        -o 1029_0I-N1-BDNA-1_B23KGCJLT4_1.bam \
        --sam2bam \
        -

# When outputting CRAM, sentieon util sort may create a spurious .bai index.
# Remove it since CRAM uses .crai indices.
if [ "bam" = "cram" ] && [ -f "1029_0I-N1-BDNA-1_B23KGCJLT4_1.cram.bai" ]; then
    rm -f 1029_0I-N1-BDNA-1_B23KGCJLT4_1.cram.bai
fi

cat <<-END_VERSIONS > versions.yml
"DAQ:SENTIEON_BWAMEM":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS