I am trying to find closest genes to my chip-seq peaks with bedtools closest. You may know that if peak and gene have overlap the distance would be zero. So it means that this computation depends on the length of peaks. Since the MAX and MIN length of my peaks are in far range, I was wondering which would be more precise? Computing distance from total length of peaks or summits of peaks to genes? I know that the summit is center of peaks but I am not sure that summit is exact location of TF binding sites.
Thanks for any discussion.
2 answers
Summits of peaks would be preferred. Actually, the resolution of your analysis would be perfect if it can be controlled within 100-500bp. Do not hope it can be very accurate. if you want to make sure the accurate location for a specific gene, you need to do further precise experiments.
•
written
by
Shicheng Guo
It depends from what you are looking for.
The help page for bedtools will probably answer your question (or part of it)
Log in to answer this question.
-
Question about constructing a consensus peak file for diffTF
written by lilithslily •
-
Differential binding analysis between different TFs
written by Alexandra •
-
Summit Extension in Cut&tag analysis
written by black_blue •
-
analyze enrichment of SNP in ChIP-seq peaks
written by BenHu •
-
Using MEMEChip for finding transcription factor binding sites motifs
written by may6343 •
-
Is it possible to control the coordinate of the peak that HOMER takes into account to annotate?
written by nanoide
-
MEME-chip: selecting the input sequences to be equal length
written by guixien •
-
Tool: ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization
written by Guangchuang Yu
-
Combining RNA seq and H3K27ac ChIP seq data
written by thomasterTW •
-
Tools To Define Cis-Regulatory Modules And Associate Them To A Nearby Gene?
written by Ahdf-Lell-Kocks •
