Question: Cuffdiff "No such file or directory" error
0
gravatar for jammydodger123456
9 days 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

And;

/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 • 194 views
ADD COMMENTlink written 9 days ago by jammydodger1234560

/EggB1/tophat_out/accepted_hits.bam

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

ADD REPLYlink written 9 days ago by genomax32k

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

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

# send email when job begins end or aborts
#$ -M dwells01@qub.ac.uk -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 9 days ago by jammydodger1234560
1

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 9 days ago • written 9 days ago by genomax32k

That worked! Thank you

ADD REPLYlink written 5 days ago by jammydodger1234560
Please log in to add an answer.

Help
Access

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