SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cuffcompare output nat Bioinformatics 11 03-11-2015 04:27 AM
no error and no output from cuffcompare hibachings2013 Bioinformatics 1 03-02-2015 03:59 PM
cuffcompare output statistic yasu Bioinformatics 2 01-05-2013 04:06 PM
cufflinks cuffcompare output Mark Bioinformatics 1 07-19-2010 08:23 AM
cuffcompare output and absent transfrags PFS Bioinformatics 0 05-12-2010 08:07 AM

Reply
 
Thread Tools
Old 06-18-2012, 09:01 AM   #1
bruce01
Senior Member
 
Location: .

Join Date: Mar 2011
Posts: 157
Default Cuffcompare - no output, no warnings

Hi,

I am running cuffcompare after making transcripts.gtf and merging those for each of the two conditions I wish to investigate. I am running without a gtf reference (I want to see the difference between running with and without -g ref.gtf makes). I am using cuffcompare 2.0.0:

cuffcompare -o ~/cufflinks/cuffComp/no_gtf/ -V -s /data/ref.gtf
/cufflinks/N_merged.cufflink_no_gtf/transcripts.gtf
/cufflinks/P_merged.cufflink_no_gtf/transcripts.gtf

I get no output at all. Verbose gives me the following:

Warning: Could not connect to update server to verify current version. Please check at the Cufflinks website (http://cufflinks.cbcb.umd.edu).
Prefix for output files:
Processing qfile #1: /cufflinks/N_merged.cufflink_no_gtf/merged.gtf
Loading transcripts from /cufflinks/N_merged.cufflink_no_gtf/merged.gtf..
Processing qfile #2: /cufflinks/P_merged.cufflink_no_gtf/merged.gtf
Loading transcripts from /cufflinks/P_merged.cufflink_no_gtf/merged.gtf..
Tracking transcripts across 2 query files..
Cleaning up..
Done.

This takes ~30 seconds. Tried with transcripts.gtf also. Also tried with another set of gtf files. Running with a ref.gtf gives:

<list of duplicates>
61 duplicate reference transcripts discarded.
..reference annotation loaded
Prefix for output files:
Processing qfile #1: /cufflinks/N_merged.cufflink_gtf/merged.gtf
Loading transcripts from /cufflinks/N_merged.cufflink_gtf/merged.gtf..
2134 redundant cufflinks transfrags discarded.
Found 11538 transcripts with undetermined strand.
Processing qfile #2: /cufflinks/P_merged.cufflink_gtf/merged.gtf
Loading transcripts from /cufflinks/P_merged.cufflink_gtf/merged.gtf..
1145 redundant cufflinks transfrags discarded.
Found 11302 transcripts with undetermined strand.
Tracking transcripts across 2 query files..
Cleaning up..
Done.


Any advice or help would be very much appreciated
bruce01 is offline   Reply With Quote
Old 03-02-2015, 03:59 PM   #2
arkfang
Junior Member
 
Location: Pittsburgh

Join Date: Mar 2015
Posts: 2
Default

Because you gave a output directory while cuffcompare asked for a output prefix. Which means you give a blank prefix to cuffcompare.
So the outputs that are supposed to be output_prefix.combined.gtf are now just .combined.gtf.
And in a lot of linux-like systems, .* file (a file starts with dot) is treated as a hidden file.
Check your output directory and you probably will see them there, hidden.
arkfang 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 09:01 PM.


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