Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Most distances in the flank_gene_distances column are 0 #235

Open
rejo27 opened this issue Mar 23, 2024 · 0 comments
Open

Most distances in the flank_gene_distances column are 0 #235

rejo27 opened this issue Mar 23, 2024 · 0 comments

Comments

@rejo27
Copy link

rejo27 commented Mar 23, 2024

Dear Pr Guangchuang Yu
I am using the ChIPseeker package to annotate the peaks of ATACseq.
The command is

txdb <- makeTxDbFromGFF("/Users/ltan/Documents/my.gff",format="gff")
peaks <- "/Users/ltan/Downloads/ATAC.peak.bed"
peakAnno <- annotatePeak(peaks, TxDb = txdb,tssRegion = c(-2000, 0), addFlankGeneInfo=T,flankDistance = 5000)

But I found that although there are multiple genes within 5000 bp upstream and downstream of some peaks, most of the flank_gene_distances columns are 0.
such as
image

I don't know why.
Looking forward to your reply

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant