SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Merging files with Cuffmerge SHeaph Bioinformatics 0 02-02-2013 07:42 AM
transcript strand of cufflink and cuffmerge result ZhigangLi RNA Sequencing 1 01-15-2013 07:10 PM
Create new Reference by merging SNV list with reference genome rdoan Bioinformatics 0 10-12-2012 08:17 AM
Cufflinks merging more than one transcript on bacterial genomes Noa Bioinformatics 2 01-24-2012 01:33 PM
Cufflinks merging more than one transcript on bacterial genomes Noa RNA Sequencing 0 01-24-2012 06:41 AM

Reply
 
Thread Tools
Old 06-27-2013, 12:08 PM   #1
mra
Junior Member
 
Location: USA

Join Date: Feb 2012
Posts: 7
Default Cuffmerge: merging reference and transcript GTFs?

Dear All.

After running cufflinks

Code:
cufflinks -p 10 --frag-len-mean 400 --frag-len-std-dev 100 -N --GTF-guide ${gtf} --frag-bias-correct ${fasta} --multi-read-correct -o cufflinks_all all.accepted_hits.bam
I got a 12M transcripts.gtf file that I eventually merged with a 240M reference GTF file using cuffmerge:

Code:
cuffmerge -p 10 -o cuffmerge_all --ref-gtf ${gtf} --ref-sequence ${fasta} assemblies.txt
(on assemblies.txt there is only the transcripts.gtf produced previously by cufflinks)

Yet, cuffmerge outputs a 142M merged.gtf. Assuming that cuffmerge makes a non-redudant union of the transcripts and reference GTF files, shouldn't merged.gtf weight at least as the reference GTF (240M)?

Thanks in advance
mra 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 11:46 AM.


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