SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Splitting a bed file in multiple bed files by track albireo Bioinformatics 3 12-18-2014 11:58 AM
GTF convert to GFF xiongdianguang Bioinformatics 3 12-25-2012 05:16 PM
GFF to GTF, and GTF to GRanges objects lewewoo Bioinformatics 2 04-03-2012 02:52 PM
gff3,gtf to gff parulvk Bioinformatics 2 11-15-2011 11:48 AM
GFF to GTF gen2prot Bioinformatics 9 12-14-2010 10:07 AM

Reply
 
Thread Tools
Old 07-04-2014, 12:00 AM   #1
StevenVanLaere
Junior Member
 
Location: Belgium

Join Date: Jun 2014
Posts: 3
Question Generation of BED/GTF/GFF-files

Dear All,

I am relatively new to RNA-seq data analysis. I did my share of literature but theory and practice are two different things. So I encounter small problems/questions on which I cannot find straight forward answers. Here is one, hope you can help.

I have aligned my reads to the reference genome (hg19) with TopHat2 and now I want to use DeSeq(2) to identify differentially expressed genes. Obviously a genome annotation file (in this case GFF3) is needed. I wonder what the best solution is. Getting the annotation file from e.g. UCSC? Or generate one from the BAM-files I have, thus converting BAM to BED to GFF3 (to GTF)? In the latter case, how do I deal with the fact that I have multiple BAM-files, i.e. one per sample. I expect differential gene expression so I guess that a GFF3 generated from a BAM file of condition 1 will be different from those for condition 2.

Thanks for your help.

Steven
StevenVanLaere is offline   Reply With Quote
Old 07-04-2014, 12:25 PM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,077
Default

Where did you get your hg19 genome from? It would be best to get the GTF file from the same source.

One solution is to get GTF file from iGenomes. The downloads are large but the files (sequences, annotation, indexes are all in sync as far as the names etc goes) and this will save you time down the road by avoiding problems with annotations etc.
GenoMax is online now   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 09:19 AM.


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