How to create GRanges for TES
1
0
Entering edit mode
@marpieraccioli-23381
Last seen 3.1 years ago
Italy

Hi, I'm interested in analyzing average density of reads from CLIP-seq experiments in the 3'UTR and around TES (transcription end site) of genes. For ChIP-seq analysis, the promoters function of "GenomicFeatures" package computes a GRanges object that spans the promoter region around the transcription start site for the transcripts in a TxDb object.Subsequent upstream 0and downstream arguments define the number of bases upstream and downstream from the transcription start site that make up the promoter region. How to create a GRanges object that spans region around TES (transcription end site) of the genes? Thank you, Marco

GenomicFeatures TxDb CLIP-Seq • 1.2k views
ADD COMMENT
0
Entering edit mode
@michael-lawrence-3846
Last seen 2.9 years ago
United States

I guess you could invert the strand and use promoters(). Another option is flank(gr, width, both=TRUE, start=FALSE) but note that the width is divided equally upstream and downstream (it's symmetric).

ADD COMMENT
0
Entering edit mode

Thank you so much for your quick reply! I tried creating a GRanges object with a list of hg19 TES ( +/- x base pairs upstream/downstream of TES position). Is it a good solution?

ADD REPLY

Login before adding your answer.

Traffic: 240 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6