#!/bin/bash -Ceuo pipefail
samtools \
merge \
-c \
--threads 15 \
--write-index \
\
\
1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam##idx##1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bai \
1/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 2/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 3/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 4/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 5/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 6/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 7/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam 8/1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bam
# Touch index to ensure timestamp is newer than aligned file
touch 1136_4Z5-T1-TDNA-1_A23MHWWLT4_1.bai
cat <<-END_VERSIONS > versions.yml
"DAQ:SAMTOOLS_MERGE":
samtools: $(echo $(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*$//')
END_VERSIONS