Cuffdiff "No such file or directory" error
0
0
Entering edit mode
6.7 years ago

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

cuffdiff cufflinks RNA-Seq • 1.8k views
ADD COMMENT
0
Entering edit mode

/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 REPLY
0
Entering edit mode

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 REPLY
1
Entering edit mode

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 REPLY
0
Entering edit mode

That worked! Thank you

ADD REPLY

Login before adding your answer.

Traffic: 2291 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6