SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problems with cuffmerge Nnogueira Bioinformatics 23 06-25-2016 10:38 PM
Cuffmerge Error ercfrtz Bioinformatics 12 05-19-2014 10:02 PM
cuffmerge failed wall_y Bioinformatics 4 07-12-2012 07:25 AM
Questions about cuffmerge Jolin Bioinformatics 0 10-04-2011 01:46 AM
cuffmerge fail papori Bioinformatics 0 07-31-2011 04:01 PM

Reply
 
Thread Tools
Old 02-02-2012, 05:28 AM   #1
ccstaats
Member
 
Location: Porto Alegre

Join Date: Jan 2012
Posts: 12
Default cuffmerge problems

Hello
I'm founding some problems with cuffmerge
I've a fasta file of 27 supercontigs from broad institute, which I edited to this format
>Supercontig_1.1

I've downloaded also the gtf file from Broad
As an example, the first lane is
Supercontig_1.1 CNB1_FINAL_CALLGENES_3 start_codon 5909 5911 . - 0 gene_id "CNBG_0001"; transcript_id "CNBG_0001T0";

I ran bowtie-build, tophat and cufflinks on my fastq files and everything is fine
However, when I run cuffmerge, I got the message:
You are using Cufflinks v1.3.0, which is the most recent release.
Warning: couldn't find fasta record for '>Supercontig_1.1'!
Warning: couldn't find fasta record for '>Supercontig_1.10'!
Warning: couldn't find fasta record for '>Supercontig_1.11'!

Moreover, I got a merged gtf file with 2 entries for the same sequence. For example:

>Supercontig_1.1 Cufflinks exon 13250 13293 . + . gene_id "XLOC_000001"; transcript_id "TCONS_00000001"; exon_number "1"; gene_name "CNBG_0004"; oId "CNBG_0004T0"; nearest_ref "CNBG_0004T0"; class_code "="; tss_id "TSS1";

Supercontig_1.1 Cufflinks exon 13233 13293 . + . gene_id "XLOC_006239" transcript_id "TCONS_00006365" exon_number "1" oId "CUFF.148.1" class_code "u" tss_id "TSS6299"

I noted that in the first case, a ">" commes before the Supercontig. So I replace every entry in gtf adding or reoving the ">" before the "Supercontig"
Nothing worked.
Any clues?
Thanks
Charley
ccstaats 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:04 PM.


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