Go Back   SEQanswers > Bioinformatics > Bioinformatics

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

Thread Tools
Old 02-02-2012, 04:28 AM   #1
Location: Porto Alegre

Join Date: Jan 2012
Posts: 12
Default cuffmerge problems

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

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?
ccstaats 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 05:01 PM.

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