scrnaseq icon indicating copy to clipboard operation
scrnaseq copied to clipboard

Crashes due to trying to unzip file that is not zipped

Open esrice opened this issue 4 months ago • 0 comments

Description of the bug

I specified an un-gzipped whitelist file with the --barcode_whitelist parameter. In the STAR_ALIGN step, it tries to unzip this file, which causes gzip to crash, which causes the step to fail. This is the offending line of .command.sh:

--soloCBwhitelist <(gzip -cdf 3M-february-2018.txt)

I will try to fix and submit a PR in the next day or two.

Command used and terminal output

$ nextflow run nf-core/scrnaseq
-profile singularity
--input ../samples.csv
--fasta: ../../ref/bGalGal1b_modified.fa
--gtf: ../../ref/bGalGal1b_modified_filtered.gtf
--protocol 10XV3
--aligner star
--outdir out
--barcode_whitelist /mnt/pixstor/data/esrbhb/3M-february-2018.txt
--save_reference

ERROR ~ Error executing process > 'NFCORE_SCRNASEQ:SCRNASEQ:STARSOLO:STAR_ALIGN (D2)'

Caused by: Process NFCORE_SCRNASEQ:SCRNASEQ:STARSOLO:STAR_ALIGN (D2) terminated with an error exit status (104)

Command executed:

STAR
--genomeDir star
--readFilesIn D2_S2_L001_R2_001.fastq.gz D2_S2_L001_R1_001.fastq.gz
--runThreadN 16
--outFileNamePrefix D2.
--soloCBwhitelist <(gzip -cdf 3M-february-2018.txt)
--soloType CB_UMI_Simple
--soloFeatures Gene
--soloUMIlen 12

--sjdbGTFfile bGalGal1b_modified_genes.gtf
--outSAMattrRGline ID:D2 'SM:D2'

--readFilesCommand zcat --runDirPerm All_RWX --outWigType bedGraph --twopassMode Basic --outSAMtype BAM SortedByCoordinate \

if [ -f D2.Unmapped.out.mate1 ]; then mv D2.Unmapped.out.mate1 D2.unmapped_1.fastq gzip D2.unmapped_1.fastq fi if [ -f D2.Unmapped.out.mate2 ]; then mv D2.Unmapped.out.mate2 D2.unmapped_2.fastq gzip D2.unmapped_2.fastq fi

if [ -d D2.Solo.out ]; then # Backslashes still need to be escaped (https://github.com/nextflow-io/nextflow/issues/67) find D2.Solo.out ( -name ".tsv" -o -name ".mtx" ) -exec gzip {} ; fi

cat <<-END_VERSIONS > versions.yml "NFCORE_SCRNASEQ:SCRNASEQ:STARSOLO:STAR_ALIGN": star: $(STAR --version | sed -e "s/STAR_//g") END_VERSIONS

Command exit status: 104

Command output: (empty)

Command error: INFO: Environment variable SINGULARITYENV_TMPDIR is set, but APPTAINERENV_TMPDIR is preferred INFO: Environment variable SINGULARITYENV_NXF_TASK_WORKDIR is set, but APPTAINERENV_NXF_TASK_WORKDIR is preferred INFO: Environment variable SINGULARITYENV_NXF_DEBUG is set, but APPTAINERENV_NXF_DEBUG is preferred gzip: invalid magic

EXITING because of FATAL ERROR: CB whitelist file /dev/fd/63 is empty. SOLUTION: provide non-empty whitelist.

Oct 11 07:18:04 ...... FATAL ERROR, exiting

Work dir: /mnt/pixstor/warrenwc-lab/users/edward/nxf_work/1a/6d24d1b3b8d570f7e134a16d877d51

Tip: you can replicate the issue by changing to the process work dir and entering the command bash .command.run

-- Check '.nextflow.log' file for details ERROR ~ Pipeline failed. Please refer to troubleshooting docs: https://nf-co.re/docs/usage/troubleshooting

-- Check '.nextflow.log' file for details -[nf-core/scrnaseq] Pipeline completed with errors- WARN: Killing running tasks (1)

Relevant files

No response

System information

  • Nextflow version: 24.04.3
  • nf-core/scrnaseq version: v2.7.1-g4171377
  • slurm executor
  • singularity profile
  • CentOS

esrice avatar Oct 11 '24 12:10 esrice