Question: Cuffmerge cannot find transcripts.gtf files
0
gravatar for calubr
2.7 years ago by
calubr0
calubr0 wrote:

Hi Everyone, I'm new to RNAseq and I'm trying to run cuffmerge but I'm encountering the following error message:

[Mon Aug 29 12:09:13 2016] Beginning transcriptome assembly merge

[Mon Aug 29 12:09:13 2016] Preparing output location ./merged_asm/

Traceback (most recent call last):

File "/software/cufflinks/2.2.1/bin/cuffmerge", line 580, in <module>

  sys.exit(main())

File "/software/cufflinks/2.2.1/bin/cuffmerge", line 538, in main

  gtf_input_files = test_input_files(transfrag_list_file)

File "/software/cufflinks/2.2.1/bin/cuffmerge", line 268, in test_input_files

  g = open(line,"r")

IOError: [Errno 2] No such file or directory: './path/to/1_clout/transcripts.gtf'

My cuffmerge command is:

cuffmerge -p 8 -g /path/to/gff-file.gff -s /path/to/genome.fa assemblies.txt

and assemblies.txt is a text file that looks like:

./path/to/1_clout/transcripts.gtf

./path/to/2_clout/transcripts.gtf

./path/to/3_clout/transcripts.gtf

./path/to/4_clout/transcripts.gtf

./path/to/5_clout/transcripts.gtf

./path/to/6_clout/transcripts.gtf

It seems to me that cuffmerge can't find the .gtf files made by cufflinks. Is my syntax wrong? I used a reference gff file when running tophat and cufflinks and I'm wondering if my transcripts.gtf files are now wrong?

Thanks for your help!

cuffmerge rna-seq cufflinks • 1.2k views
ADD COMMENTlink modified 2.7 years ago by genomax67k • written 2.7 years ago by calubr0

If you are proving full directory paths for files then you would not use the period (.) at the beginning of those paths (like you have in your example for assemblies.txt file above).

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by genomax67k
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: 924 users visited in the last hour