File Info

Filename
.command.sh
Full Path
s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/94/5282eae97aa5e33706f1a15ca271e6/.command.sh
Size
1.2 KB
Task
NFCORE_RNAFUSION:RNAFUSION:CTATSPLICING_STARTOCANCERINTRONS (659_b4x-T1-TRNA-1_B23MHV2LT4_1)
Attempt

Content

View in new tab 
#!/bin/bash -euo pipefail
python -u /usr/local/src/CTAT-SPLICING/STAR_to_cancer_introns.py \
    --SJ_tab_file 659_b4x-T1-TRNA-1_B23MHV2LT4_1.SJ.out.tab \
    --chimJ_file 659_b4x-T1-TRNA-1_B23MHV2LT4_1.Chimeric.out.junction \
    --bam_file 659_b4x-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam \
    --output_prefix 659_b4x-T1-TRNA-1_B23MHV2LT4_1 \
    --ctat_genome_lib ctat_genome_lib_build_dir \
    --vis --sample_name 659_b4x-T1-TRNA-1_B23MHV2LT4_1

# Sort output files to ensure consistent ordering across runs while preserving headers
if [ -f 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns ]; then
    # Extract header (first line) and sort data lines by fourth column (uniq_mapped) descending, then by first column (intron)
    head -n 1 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns > 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.tmp
    tail -n +2 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns | LC_ALL=C sort -k4,4nr -k1,1 >> 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.tmp
    mv 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.tmp 659_b4x-T1-TRNA-1_B23MHV2LT4_1.cancer.introns
fi

cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:CTATSPLICING_STARTOCANCERINTRONS":
    ctat-splicing: 0.0.3
END_VERSIONS