Hoohm / dropSeqPipe

A SingleCell RNASeq pre-processing snakemake workflow
Creative Commons Attribution Share Alike 4.0 International
147 stars 47 forks source link

Error in rule repair #116

Open Eigertje opened 2 years ago

Eigertje commented 2 years ago

I'm running DropSeqPipe v5.0 with 5 samples in my samples.csv file. However, the pipeline aborted on 1 sample during "rule repair". After rerunning the pipeline with this sample removed the error occured again on another sample. For the other 3 samples the whole DropSeqPipe pipeline worked fine. Any idea what is happening here? I'm not familiar at all with python and snakemake I'm afraid.

The error: `Error in rule repair: jobid: 82 output: ../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz, ../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz log: ../Results/logs/bbmap/HN4_RNA_repair.txt (check log file(s) for error message) conda-env: /mnt/go/temp/scRNA-seq/dropSeqPipe/.snakemake/conda/65b8583675236a3c6fd8572f06a27964 shell: repair.sh -Xmx20g in=../Results/samples/HN4_RNA/trimmed_R1.fastq.gz in2=../Results/samples/HN4_RNA/trimmed_R2.fastq.gz out1=../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz out2=../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz repair=t threads=4 2> ../Results/logs/bbmap/HN4_RNA_repair.txt (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

Removing output files of failed job repair since they might be corrupted: ../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz, ../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz `

The content of the "/mnt/go/temp/scRNA-seq/Results/logs/bbmap/HN4_RNA_repair.txt" file: `java -ea -Xmx20g -cp /mnt/go/temp/scRNA-seq/dropSeqPipe/.snakemake/conda/65b8583675236a3c6fd8572f06a27964/opt/bbmap-38.22-1/current/ jgi.SplitPairsAndSingles rp -Xmx20g in=../Results/samples/HN4_RNA/trimmed_R1.fastq.gz in2=../Results/samples/HN4_RNA/trimmed_R2.fastq.gz out1=../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz out2=../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz repair=t threads=4 Executing jgi.SplitPairsAndSingles [rp, -Xmx20g, in=../Results/samples/HN4_RNA/trimmed_R1.fastq.gz, in2=../Results/samples/HN4_RNA/trimmed_R2.fastq.gz, out1=../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz, out2=../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz, repair=t, threads=4]

Set threads to 4 Set INTERLEAVED to false Started output stream. /mnt/go/temp/scRNA-seq/dropSeqPipe/.snakemake/conda/65b8583675236a3c6fd8572f06a27964/bin/repair.sh: line 104: 18758 Killed java -ea -Xmx20g -cp /mnt/go/temp/scRNA-seq/dropSeqPipe/.snakemake/conda/65b8583675236a3c6fd8572f06a27964/opt/bbmap-38.22-1/current/ jgi.SplitPairsAndSingles rp -Xmx20g in=../Results/samples/HN4_RNA/trimmed_R1.fastq.gz in2=../Results/samples/HN4_RNA/trimmed_R2.fastq.gz out1=../Results/samples/HN4_RNA/trimmed_repaired_R1.fastq.gz out2=../Results/samples/HN4_RNA/trimmed_repaired_R2.fastq.gz repair=t threads=4 `