SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
cuffmerge crashes when converting gtf files to sam files swbiggs4 Bioinformatics 20 02-16-2017 09:19 AM
Demultiplexing dual-indexed MiSeq fastq files lynchde Bioinformatics 2 08-18-2013 02:15 PM
Merging files with Cuffmerge SHeaph Bioinformatics 0 02-02-2013 06:42 AM
Can You mix indexed with un-indexed library on MiSeq run ? memento Illumina/Solexa 2 01-22-2013 11:00 AM
bwa indexed files? sbsusers Bioinformatics 2 06-09-2012 06:59 AM

Reply
 
Thread Tools
Old 10-23-2014, 03:24 AM   #1
Marys
Junior Member
 
Location: Warsaw

Join Date: Oct 2014
Posts: 4
Default cuffmerge - from wrong not indexed files?

hi Everyone,
I am new in RNAseq analysis. I have run tophat then cufflinks, and now I wnat to run cuffmerge, but I have found this problem:

# cuffmerge -p 4 -g /path/to/genes.gtf -s /path/to/genome.fa /path/to/assemblies.txt

[Thu Oct 23 13:18:51 2014] Beginning transcriptome assembly merge
-------------------------------------------

[Thu Oct 23 13:18:51 2014] Preparing output location ./merged_asm/
Traceback (most recent call last):
File "/usr/bin/cuffmerge", line 580, in <module>
sys.exit(main())
File "/usr/bin/cuffmerge", line 530, in main
transfrag_list_file = open(args[0], "r")
IOError: [Errno 2] No such file or directory: '/home/michalina/BAM_SAM_tophat_samtools/TS2/assemblies.txt


I have no idea what is going wrong. My cufflinks were wrong? gtf or fa references file are not appropriate? hese are the same reference files that I use for previous analysis (tophat, cufflinks).

many thanks for any commnet or advise
M.
Marys is offline   Reply With Quote
Old 10-23-2014, 03:32 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,077
Default

Have you made a list of the GTF files from your cufflinks assemblies and then put the names in the "assemblies.txt" file? Your error is about this file not being present.

From Cuffmerge help:

Quote:
<assembly_list.txt> Text file "manifest" with a list (one per line) of GTF files that you'd like to merge together into a single GTF file.
GenoMax is offline   Reply With Quote
Old 10-23-2014, 03:48 AM   #3
Marys
Junior Member
 
Location: Warsaw

Join Date: Oct 2014
Posts: 4
Default

Hi Genomax,

no I haven't, how can I do this? could you tell me?

I did something like this:
cat /home/michalina/BAM_SAM_tophat_samtools/TS/transcripts.gtf /home/michalina/BAM_SAM_tophat_samtools/TS2/transcripts.gtf > asseblies.txt

but I was afraid it was wrong.

Honestly, I am new one in all these things...
Marys is offline   Reply With Quote
Old 10-23-2014, 04:03 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,077
Default

Make a file called "assemblies.txt" and put these two lines in it (use a text editor):

Quote:
/home/michalina/BAM_SAM_tophat_samtools/TS/transcripts.gtf
/home/michalina/BAM_SAM_tophat_samtools/TS2/transcripts.gtf
GenoMax is offline   Reply With Quote
Old 10-23-2014, 04:24 AM   #5
Marys
Junior Member
 
Location: Warsaw

Join Date: Oct 2014
Posts: 4
Default

wow! many thanks!!! that was that much simple
Marys is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:56 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO