Thanks Hans-Rudolf

John




________________________________
 From: Hans-Rudolf Hotz <hrh@fmi.ch>

onductor@r-project.org> 
Sent: Thursday, November 21, 2013 5:38 AM
Subject: Re: [BioC] exon genomic coordinates


Hi John

You can use the BioMart database, which you can access with the biomaRt 
package to get all exons for all transcripts for a given giene, eg:

library(biomaRt)
ensembl = useMart("ensembl")
#assuming you are interested in mouse
mouse.ensembl = useDataset("mmusculus_gene_ensembl",mart=ensembl)

getBM(attributes = 
c("chromosome_name","exon_chrom_start","exon_chrom_end","ensembl_exon_id","ensembl_transcript_id","ensembl_gene_id"), 
filters = 'mgi_symbol', values=c("KIT"),mart=mouse.ensembl)


Hope this helps

Hans-Rudolf


On 11/21/2013 09:14 AM, array chip wrote:
> Hi,
>
>
> Can anyone suggest how to retrieve the genomic coordinates for all exons for a given gene by say gene symbol? For example, how to retrieve the coordinates for all 21 exons for gene KIT?
>
> Thanks
>
> John
>     [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

>
	[[alternative HTML version deleted]]

