cuffmerge warnings Warning: couldn't find fasta record for 'chr1_GL456211_random'!

I got a cuffmerge warning like this

Warning: couldn't find fasta record for 'chr1_GL456211_random'!

Warning: couldn't find fasta record for 'chr1_GL456221_random'!

Warning: couldn't find fasta record for 'chrUn_JH584304'!

Warning: couldn't find fasta record for 'chrX_GL456233_random'!

After search, I noticed that this is a common warning due to the inconsistent names.

It is a warning, not error. And it still gives me the output files.

My questions are:

1 is the output file incomplete or not?

2 If I do "$grep -v random genes.gtf > genes.grep.gtf" and then use this gtf file for cuffmerge, is it OK to do it?

3 How to fix these warnings

Thanks so much~~

rna-seq

• 5,084 views • link

updated 10.3 years ago by Devon Ryan

• written 10.3 years ago by jolin0701-dy

1 answer

It's complete. Were it not, then those warnings would be errors instead.
That will get rid of most of the warnings. The better solution would be to just download the fasta file that includes the unplaced contigs from UCSC.
See above

• 0 views • link

written 10.3 years ago by Devon Ryan

More posts like this

Cuffdiff - Error in running cuffdiff by using the merged transcriptome assembly along with the BAM …
written 6.9 years ago by evashan96 •

hello everyone, I ran cuffdiff using following command: cuffdiff -o ./diff_out -b ./genome.fa -p 7 -L c1,c2 -u ./merged.gtf \ ./accepted_hits_S1.bam,./accepted_hits_S2.bam,./accepted_hits_S3.bam ./accepted_hits_S4.bam,./accepted_hits_S5.bam,./accepted_hits_S6.bam,./accepted_hits_S7.bam,./accepted_hits_S8.bam,./accepted_hits_S9.bam,./accepted_hits_S10.bam,./accepted_hits_11.bam,./accepted_hits_S12.bam And the output …
Warning with cuffmerge: couldn't find fasta record
written 9.0 years ago by irving.jgalo23 •

Hi im trying to run cuffmerge with my gtf anotation using the next command cuffmerge -p 8 -g TAIR_gff/TAIR10_GFF3_genes.gff -s TAIR10_chr_all_new.fa assemblies.txt But I get …
Issue with gffread during converting gtf file to fasta
written 9.2 years ago by seta

Hi all, I used cufflinks package for genome-guided transcriptome assembly, followed by using gffread for converting the produced gtf file to fasta file, but some …
Cufflinks segmentation fault
written 9.2 years ago by murphy.charlesj

I' am running cufflinks on some RNA-seq: cufflinks -o . -b Mus_musculus.GRCm38.dna.primary_assembly.fa -G Mus_musculus.GRCm38.79.gtf test.40m.bam I get the following error: You are using Cufflinks v2.2.1, …
Could not find fasta record: CuffCompare & Cuffmerge
written 10.2 years ago by EpiExplorer

Hi , I am using the genes.gtf and genome.fa files from iGenomes for Homo sapiens Ensemble GRCh37, for cuffcompare and cuffmerge. But I get the …
error in cuffmerge
written 10.7 years ago by zizigolu

Hi, I downloaded gtf and genome.fa for yeast already but after trying many fastq files again I get the same error May you tell me …
Cuffdiff terminated after GffObj::getSpliced() error
written 10.7 years ago by nalandaatmi •

Dear All, I am currently analyzing mouse RNASeq samples. I used Tophat2.1.1 version, cufflinks2.2.1 version for my analysis and I downloaded the latest version of …
How could I use the Ensembl ncRNA fasta file to annotate my RNA-seq data
written 11.8 years ago by i19870503 •

I have use tophat and cufflinks get the out put, but I do not understand how to use ncRNA fasta file (from Ensembl database) annotate …
Cuffmerge Error : Couldn'T Find Fasta Record
written 13.9 years ago by zahidawani86 •

<p>Hi There, I am newbie to RNAseq analysis with basic knowledge of linux and perl. I am using tutorial from "Differential gene and transcript expression …
Tophat Protocol Error Occured At Step 4
written 14.2 years ago by ashis.csedu •

<p>Hi all, I was doing differential gene and transcript expression analysis by following the Tophat protocol (<a href='http://www.nature.com/nprot/journal/v7/n3/full/nprot.2012.016.html'>published in Nature protocol in March 2012</a>) on …