Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
cuffdiff: use merged.gtf from cuffmerge or combined.gtf from cuffcompare? turnersd Bioinformatics 21 10-02-2014 04:41 AM
Cuffcompare: difference between "nearest_ref" and "oID" wariobrega Bioinformatics 0 05-07-2014 03:49 AM
Compare "-GTF" and "-GTF-guide" of cufflink options jiyan Bioinformatics 1 01-23-2014 09:42 AM
Inquery: not-stranded transcripts ("." in strand column) in cufflinks output gtf file sunfuhui Bioinformatics 0 11-15-2012 03:43 PM
Filtering "merge.gtf" (from cuffmerge) by redundance Ramon Vidal Bioinformatics 3 02-13-2012 07:47 AM

Thread Tools
Old 10-17-2014, 08:36 AM   #1
Junior Member
Location: Rockville MD

Join Date: Oct 2014
Posts: 1
Default Cuffcompare merged.gtf to self not all "="

I have been working on comparing the results of 20 merged.gtf files to see the differences between each pair of files. When doing this work, I used the following type of command:
cuffcompare -p 4 -r merged_1.gtf -o ./compare/merged2_vs_1 merged_2.gtf

This should use merged_1.gtf as the "reference" and compare merged_2.gtf to it. However, I noticed that if merged_1.gtf and merged_2.gtf are switched the results are different. Furthermore, if I compare merged_1.gtf to itself ( merged_1.gtf), not all of the class codes are "=". I actually get something like this:

"Class_code" #_occurrences_in_.tmap_file
"=" 215622
"u" 11040
"p" 1172
"c" 1051
"i" 1707

How can this be the case? if I compare something to itself, I should expect all class codes to be "=".

Any thoughts?
travis_jensen123 is offline   Reply With Quote


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 02:59 PM.

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