Setting Peak Cutoff Using Minimum # of Reads
Is there standard way to set a cutoff for minimum number of reads to consider something an H3K4Me3 peak? The issue we face is that peaks with very few reads often have the greatest fold change in control vs treatment. However, when examining these peaks in IGV sometimes they look real with very few reads in the vicinity, but other times they look like background.
Thank you! Andrew
• 1,453 views
•
link
written
by
Andrew Siladi •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Statistical test to look for peak differentiation between two manhattan plots
written by jfaberha •
-
Cut&Run TF Analysis - Very Low Peak Counts Despite Good Library Quality
written by sk •
-
Is it acceptable to use different FDR cutoffs for peak calling in different ChIP-seq experiments (w…
written by science_lizard •
-
ChIP-Seq experiment, MACS2 bdgpeakcall 0 fold-change, pvalue and qvalue
written by Bastien Hervé
-
Sanity check: DESeq2 analysis of RNAseq data gives atypical plots
written by pemiliano.canton •
-
Peak calling with MACS for paired end DNase-seq data
written by Bioinformatist Newbie
-
how to find differential region between my replicates.
written by maryamforoozani16 •
-
Macs2 callpeak error: IndexError: list index out of range
written by Gary
-
MeDip-seq data interpretation + analysis / MEDIPS, MACS2
written by florian.noack •
-
Manorm For Differential Histone Binding In Chip-Seq Samples
written by Diana
That's a usual problem with differential analyses. Along with fold change, use normalized read count (e.g. RPM) as the cut off, e.g. min RPM must be >5.