Question: Cuffdiff "No such file or directory" error
gravatar for jammydodger123456
4 months ago by
jammydodger1234560 wrote:

I am trying to use the job script;

   cuffdiff -p 10 -L Egg,Metacercaria PRJEBWB9PAULAnnotations.gff3 EggB1/tophat_out/accepted_hits.bam,EggB2/tophat_out/accepted_hits.bam MetB1/tophat_out/accepted_hits.bam,MetB2/tophat_out/accepted_hits.bam,MetB3/tophat_out/accepted_hits.bam

However, I get the error;

  Error: cuffdiff requires at least 2 SAM files


/EggB1/tophat_out/accepted_hits.bam,/EggB2/tophat_out/accepted_hits.bam: No such file or directory

I have checked the paths to these files and everything appears to be as it should for my script to work.

Any suggestions?

Thanks, Duncan

rna-seq cuffdiff cufflinks • 330 views
ADD COMMENTlink written 4 months ago by jammydodger1234560


Why is cuffdiff looking for these under the / directory? Can you provide full paths for all the files?

ADD REPLYlink written 4 months ago by genomax39k

EggB1 is one directory within my working directory. Here's my script;

#$ -N EvM
#$ -o /mnt/scratch/users/40100371/Read_Files/Fhep/EggvMet.$JOB_ID

# send email when job begins end or aborts
#$ -M -m bea

#$ -pe smp-verbose 10

#load cufflinks module environment
module load apps/cufflinks

# set working directory to scratch space project folder
#$ -wd /mnt/scratch/users/40100371/Read_Files/Fhep

#cuffdiff [options]* <transcripts.gtf> <sample1_replicate1.sam,sample1_replicateN.sam> <sample2_replicate1.sam,sample2_replicateN.sam>
cuffdiff -p 10 -L Egg,Metacercaria PRJEBWB9PAULAnnotations.gff3
/EggB1/tophat_out/accepted_hits.bam,/EggB2/tophat_out/accepted_hits.bam /MetB1/tophat_out/accepted_hits.bam,/MetB2/tophat_out/accepted_hits.bam,/MetB3/tophat_out/accepted_hits.bam

I tried to run this without the / before each rep, but the same error message comes back

ADD REPLYlink written 4 months ago by jammydodger1234560

I assume these directories are in /mnt/scratch/users/40100371/Read_Files/Fhep? Can you try to explicitly say /mnt/scratch/users/40100371/Read_Files/Fhep/EggB1/tophat_out/accepted_hits.bam etc?

ADD REPLYlink modified 4 months ago • written 4 months ago by genomax39k

That worked! Thank you

ADD REPLYlink written 3 months ago by jammydodger1234560
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 986 users visited in the last hour