Question: Cuffdiff "No such file or directory" error
0
gravatar for jammydodger123456
11 months ago by
jammydodger12345610 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 • 561 views
ADD COMMENTlink written 11 months ago by jammydodger12345610

/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 11 months ago by genomax51k

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 11 months ago by jammydodger12345610
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 11 months ago • written 11 months ago by genomax51k

That worked! Thank you

ADD REPLYlink written 11 months ago by jammydodger12345610
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: 2093 users visited in the last hour