getGeneSeq: Get gene sequence using the biomaRt package

Description Usage Arguments Value Note Author(s) Examples

View source: R/getGeneSeq.R

Description

Get gene sequence using the biomaRt package

Usage

1
getGeneSeq(LocationParameters, mart)

Arguments

LocationParameters

c(ensembl_gene_id, distance from the peak to the transcription start site of the gene with the above ensemblID, upstream offset from the peak, downstream offset from the peak, Gene Start, Gene End)

mart

see useMart of bioMaRt package for details

Value

a list with the following items

feature_id

ensemble gene ID

distancetoFeature

distance from the peak to the transcriptionstart site of the gene with the above ensembl gene ID

upstream

upstream offset from the peakStart

downstream

downstream offset from the peakEnd

seq

sequence obtained around the peak with above upstream and downstream offset

Note

internal function not intended to be called directly by users

Author(s)

Lihua Julie Zhu

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
if (interactive())
{
mart <- useMart(biomart="ensembl", dataset="drerio_gene_ensembl")
LocationParameters =c("ENSDARG00000054562",400, 750, 750,40454140,40454935)
getGeneSeq(LocationParameters, mart)

LocationParameters =c("ENSDARG00000054562",752, 750, 750,40454140,40454935)
getGeneSeq(LocationParameters, mart)
 
LocationParameters =c("ENSDARG00000054562",750, 750, 750,40454140,40454935)
getGeneSeq(LocationParameters, mart)
 
 LocationParameters =c("ENSDARG00000054562",-2, 750, 750,40454140,40454935)
 getGeneSeq(LocationParameters, mart)

 LocationParameters =c("ENSDARG00000054562",0, 750, 750,40454140,40454935)
 getGeneSeq(LocationParameters, mart)
 
 LocationParameters =c("ENSDARG00000054562",2, 750, 750,40454140,40454935)
 getGeneSeq(LocationParameters, mart) 
 
 LocationParameters =c("ENSDARG00000054562",1000, 750, 750,40454140,40454935)
 getGeneSeq(LocationParameters, mart)
}

LihuaJulieZhu/ChIPpeakAnno documentation built on Aug. 5, 2020, 12:02 a.m.