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

cp .command.err GM24385_0001_RNA_0001_23H5VFLT4_s24.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