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?