Question: How to get gene annotation for a list of genes
0
gravatar for wstla27
8 months 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 • 278 views
ADD COMMENTlink written 8 months ago by wstla270
3
gravatar for RamRS
8 months ago by
RamRS26k
Houston, TX
RamRS26k 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 8 months ago by RamRS26k

Tried biomaRT, worked perfectly! Thanks a lot!

ADD REPLYlink written 8 months 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 8 months ago by RamRS26k
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: 1239 users visited in the last hour