Hello everyone,
I am studying bioinformatics and as an amateur I'd like to ask you a few questions. Right now I am focusing in a chip-seq experiment from Encode.I've downloaded the raw data from both biological replicates,I've done quality check via the Fastqc program and I've finished the mapping procedure against the reference genome (Mus Musculus).The output of the mapping was a sam file that I converted to bam via the samtools program in order to utilize it to the peak calling algorithm.Here is my issue!As I read the macs algorithm has the option to use as a parameter a control sample together with the treatment sample for finding more reliable peaks.Where can I find the control data in Encode?The chip-seq experiment I am talking about can be seen from the below url: https://www.encodeproject.org/experiments/ENCSR985ZTV/ At first, I used the alignments in bam format in the field of processed data.Then, I noticed that in the summary of the experiment there is an option named "Controls" and so here I am!This is the link for the control data: https://www.encodeproject.org/experiments/ENCSR601HSS/ Here I can see only one replicate.Could you please help me with this issue because I am really confused!
Thank you in advance for your help :)
2 answers
ENCODE had a penchant for using the same control sample for multiple things. In this case, there really is only a single control sample for both CTCF IPs, so just use that for peak calling both IP BAM files.
Thank you Devon for your response.Could you please tell me the exact name (accession code) of the control sample that I should use as control?
It's the one you linked to (ENCSR601HSS).
Yes of course but which exact file?The processed unfiltered bam file (ENCFF040OST) or the filtered one(ENCFF574XMX)?
Since you aligned the IP samples yourself, get the fastq file and process it the same way.
Ok now it is clear!Thank you very much indeed!
Should I use the same parameters that I used at the first mapping (treatment samples) or more "relaxed"? Thank you
Use the same settings.
Thank you very much Devon!
Log in to answer this question.
-
Question regarding MACS2 approach of peak calling
written by rkc5 •
-
Peak Calling on a single gene?
written by Maycon •
-
Sample size estimation for differential ChIP-Seq
written by Russ
-
ENCODE: Why are there more than one bed narrowPeak file for each DNase-seq experiment (but only one…
written by edostyles •
-
Chip-seq IgG as control or Input DNA or usage of both
written by luckysardar171 •
-
What does "Fold change over controls" from Encode ChIP-Seq experiments mean?
written by marquimedes21 •
-
Gene expression values from ENCODE RNA-seq experiments
written by fabrizio.frasca •
-
Choosing the right RNA-seq experiment: polyA mRNA vs polyA depleted
written by fabrizio.frasca •
-
False Positive Peaks In Chip-Seq - How To Locate Them And How Do They Arise?
written by Tky
-
Is It Normal To Have Very Low Mapping For Treatment File When Compared To Control File In Chip-Seq?
written by Jordan
