Cuffmerge cannot find transcripts.gtf files
0
1
Entering edit mode
7.7 years ago
calubr ▴ 10

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!

RNA-Seq cufflinks cuffmerge • 3.0k views
ADD COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 2788 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