Different position for same gene using annotatePeak from ChIPseeker
0
0
Entering edit mode
6 weeks ago
QX ▴ 60

Hi all,

I try to annotate peak with ChIPseeker, but the outcome show the start sites of a gene are differ for the same gene:

                      chr     start       end  SYMBOL geneStart   geneEnd
1:26233263-26235629     1  26233263  26235629   CEP85  26234200  26277687
1:26470149-26473811     1  26470149  26473811   HMGN2  26472440  26476642
1:27913903-27915369     1  27913903  27915369    RPA2  27892073  27914513
1:35748222-35748658     1  35748222  35748658   CLSPN  35746965  35748774
1:35768251-35770632     1  35768251  35770632   CLSPN  35736405  35769928
1:35771985-35772203     1  35771985  35772203   CLSPN  35732112  35769978
1:37573334-37573553     1  37573334  37573553    GNL2  37566817  37570893
1:37593320-37593538     1  37593320  37593538    GNL2  37590779  37594092
1:37595191-37596612     1  37595191  37596612    GNL2  37593485  37595935
1:75725766-75726947     1  75725766  75726947   ACADM  75726849  75763640
1:75741781-75743312     1  75741781  75743312   ACADM  75744318  75763483
1:93874059-93874743     1  93874059  93874743 DNTTIP2  93869780  93873369
1:93878084-93880000     1  93878084  93880000 DNTTIP2  93876759  93879156

for e.g: CLSPN have genestart at: 35746965, 35736405 , 35732112

I use the code from the tutorial:

peakAnno <- annotatePeak(macsPeaks_GR,
                        tssRegion=c(-3000, 3000),
                        TxDb=TxDb.Hsapiens.UCSC.hg38.knownGene, 
                        annoDb="org.Hs.eg.db")

Can anyone help me explain this?

ChIPseeker ATACseq • 140 views
ADD COMMENT

Login before adding your answer.

Traffic: 1372 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