Question: How to get gene annotation for a list of genes
0
gravatar for wstla27
4 weeks ago by
wstla270
wstla270 wrote:

I have a list of genes with their ensembl ids. I want to get their positions on chromosomes or functional annotations, so that I can use bedtools to align them with some bam files and get their gene expressions (counts). Is there anyway that I can do this in bulk instead of one by one?

(I looked at SnpEff but I don’t think that’s what I want. )

Many thanks!

gene • 104 views
ADD COMMENTlink written 4 weeks ago by wstla270
3
gravatar for RamRS
4 weeks ago by
RamRS23k
Houston, TX
RamRS23k wrote:

I'm assuming these are human genes.

You can use various R packages to get these data dimensions - org.Hs.eg.db, biomaRt, TxDb (I think), etc. search for "get sequence from gene symbol" (say) to get a more concrete lead on the package to try.

ADD COMMENTlink written 4 weeks ago by RamRS23k

Tried biomaRT, worked perfectly! Thanks a lot!

ADD REPLYlink written 29 days ago by wstla270

Glad it helped. If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they all work.

Upvote|Bookmark|Accept

ADD REPLYlink written 29 days ago by RamRS23k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1076 users visited in the last hour