#!/bin/bash -euo pipefail
preseq \
lc_extrap \
-verbose -seed 1 -seg_len 100000000 -bam -s 10000000 \
-pe \
-output aih-tih-sc-176637-R1_B23WHYVLT4_1.lc_extrap.txt \
aih-tih-sc-176637-R1_B23WHYVLT4_1.Aligned.sortedByCoord.out.bam

cp .command.err aih-tih-sc-176637-R1_B23WHYVLT4_1.attempt2.command.log

cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:QC_WORKFLOW:PRESEQ_LCEXTRAP":
    preseq: $(echo $(preseq 2>&1) | sed 's/^.*Version: //; s/Usage:.*$//')
END_VERSIONS