Hi!
I am new to bioinformatics and am doing it for the first time for my master's thesis. I've performed CUT&RUN between two conditions with H3K27ac and IgG as a control. I crosslinked the histone mod to DNA with PFA which is known to reduce signal.
Everything was looking fine between the two conditions, nothing strange during alignment (Bowtie2) or clean-up (remove MAPQ <10, chrM/Un/N), both conditions show very similar alignment rates, number of reads dropped, and size distributions. Even after peak calling (SEACR) the FRiP scores are very similar. However, the IGV tracks are almost completely peakless for one of the conditions, and correspondingly DiffBind showed 11 differential peaks in this condition and >5000 in the other.
What I can't seem to figure out is why the data processing doesn't give any indication of a failed experiment? Or does it, and I've just been looking in the wrong places?
Let me know if there is any relevant information missing. Thanks!
0 answers
No answers yet.
Log in to answer this question.
-
Pseudogene - scarce info
written by mf810 •
-
Help with DESeq2 output
written by randlkc •
-
DeSeq2 analysis - very low variability among replicates
written by srhic •
-
DEseq2 design matrix with 3 factors
written by kand3e •
-
How to assemble viral genomes when my data contains host DNA as well
written by GBC_Zonatos •
-
Most recent tools for circadian gene expression analysis (in plants)
written by Raito92
-
Tutorial: InsideDNA: Analysing transcription factors with ChIP-Sequencing data: a full bioinformati…
written by shilparaopradeep •
-
Forum: Visualization of ChIP-seq data using Heatmaps (Updated: 06/10/16)
written by Sinji
-
Single-cell RNA-seq coverage track not corresponding to htseq-count output or IGV gene annotations?
written by shanasabri •
-
Cross-species RNA-Seq analysis using edgeR with multiple factors?
written by dan.shea •
Are these tracks scaled or normalized? I would say in the lower one you do not see anything as the scale is too large for the present signal. Is this the BAM files or bigwigs?
Thank you for the reply! The tracks are scaled. These are bedgraph files, the output from https://github.com/Henikoff/Cut-and-Run/blob/master/spike_in_calibration.csh. I agree that the scale is too large, but I set it with regards to the tracks for condition 1.