SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat v1.1 with GTF files hyjkim Bioinformatics 7 12-17-2012 07:11 AM
Best source for GTF file for use with TopHat/Cufflinks sdarko Bioinformatics 17 12-14-2012 11:48 AM
Tophat/Cufflinks GTF (canine/non-human) M&M RNA Sequencing 0 08-01-2011 03:00 PM
Ensembl gtf to gff3 for tophat genec Bioinformatics 16 05-29-2011 09:15 PM
Acceptable Sp/Sn output from cufflinks and problems with Homo_sapiens.GRCh37.60.gtf nat Bioinformatics 0 12-02-2010 09:58 PM

Reply
 
Thread Tools
Old 12-15-2011, 06:18 PM   #1
shouguogao
Junior Member
 
Location: birminghal al

Join Date: Dec 2011
Posts: 2
Default Homo_sapiens. GRCh37.55.gtf for tophat

I tried to use tophat to generate input files for Altanalyzer

After downloading GRCh37.55.gtf from ensembl, I run

./tophat --GTF Homo_sapiens.GRCh37.55.gtf --output-dir /home/user/tophat_data/NP902 hg19 /home/user/hESC-NP/SRS011902-NP.fq

I got the error message:

Error: reading .... Homo_sapiens.GRCh37.55.gtf .

I tried to follow suggestions in this link http://seqanswers.com/forums/showthread.php?t=8247
But they do not work for my task.

Homo_sapiens.GRCh37.56.gtf and Homo_sapiens.GRCh37.65.gtf do not work either.

Any suggestions?
shouguogao is offline   Reply With Quote
Old 12-15-2011, 08:20 PM   #2
Jon_Keats
Senior Member
 
Location: Phoenix, AZ

Join Date: Mar 2010
Posts: 279
Default

You can't use the hg19 reference genome and ensembl annotations as the chromosomes names don't match (ie. chromosome 1, hg19 = chr1 while ensemble = 1). You can get all the needed files from iGenomes as suggested on the tophat website or create your own bowtie index from the ensemble reference genome file.
Jon_Keats is offline   Reply With Quote
Old 12-16-2011, 02:51 AM   #3
shouguogao
Junior Member
 
Location: birminghal al

Join Date: Dec 2011
Posts: 2
Default

Thanks for your reply

I did use bowtie index in Ensembl GRCh37 from http://tophat.cbcb.umd.edu/igenomes.html as reference genome. I just rename it to hg19. sorry for confusing.
will try to download dna sequence from ensembl and build my own reference genome.
shouguogao is offline   Reply With Quote
Reply

Tags
error, grch37.55.gtf, tophat

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 12:52 AM.


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