#!/bin/bash -euo pipefail
preseq \
lc_extrap \
-verbose -seed 1 -seg_len 100000000 -bam \
-pe \
-output 659_bTo-T1-TRNA-1_B23WHTKLT4_1.lc_extrap.txt \
659_bTo-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam

cp .command.err 659_bTo-T1-TRNA-1_B23WHTKLT4_1.attempt1.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