SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Obtain the gene names for FPKM generated from cufflinks command (http://seqanswers.com/forums/showthread.php?t=39879)

twotwo 01-13-2014 07:54 AM

Obtain the gene names for FPKM generated from cufflinks command
 
Hi,
I used the cufflinks command to generate the FPKM values. However, the genes.fpkm_tracking, which have the FPKM values, did not have the gene names. The following is my command.
cufflinks -o /mypath/cufflinks2 -g /mypath/hg19/hg19.gtf /mypath/tophat/accepted_hits.bam

I used UCSC hg19 as reference for both tophat and cufflinks.

And the following is part of the output:

tracking_id class_code nearest_ref_id gene_id gene_short_name tss_id locus
CUFF.1 - - CUFF.1 - - chr1:11873-14409
CUFF.2 - - CUFF.2 - - chr1:14361-29370
NR_026820 - - NR_026820 - - chr1:34610-36081
NM_001005484 - - NM_001005484 - - chr1:69090-70008
CUFF.3 - - CUFF.3 - - chr1:323891-328581
CUFF.4 - - CUFF.4 - - chr1:134772-140566
NM_001005224 - - NM_001005224 - - chr1:367658-368597


The gene_id kept missing. Any comments? Thanks!

Zapages 01-13-2014 06:14 PM

hmm maybe loading the files into NCBI IGV and seeing what matches up based on the chromosome location. The chromosome location are indicated.

Steps:

1) Load the your UCSC hg19 into NCBI IGV
2) Load your GTF files and BAM hits files to see the results.

I hope this helps a bit. :)

rboettcher 01-14-2014 01:23 AM

Try bedtools intersect with your Cufflinks GTF and the GTF you used as reference annotation to see if there is any known annotation in the regions identified by Cufflinks (which should not be the case).

FYI: -g means that Cufflinks will perform a RABT assembly of novel transcripts, in case you only want quantification of known genes use -G.


All times are GMT -8. The time now is 12:50 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.