Are multi mapped reads used in ChIP-sea peaking calling such as MACS2?
I wonder whether multi mapped reads are used in ChIP-sea peaking calling such as MACS2? If so, what are the suggestions on QC?
chip-seq
peaking
calling
reads
macs2
multi-mapped
• 1,342 views
•
link
updated
by
Friederike
•
written
by
chao.xing •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Run MACS2 on merged BAM file Error
written by cthangav •I'm trying to run MACS2 on a merged ATAC-seq BAM file (because I have two replicates. The Pipeline I'm using involves PECA, which doesn't allow …
-
low FRiP(Fraction of Reads in Peaks) score in ATAC-seq
written by yoonho •Hi. I'm doing ATAC-seq analysis of colon tissue. I analyzed 1)QC -> 2)Mapping -> 3)Post alignment processing(remove mt reads, duplicated reads, multi-mapped reads) -> 4)Peak …
-
Chip-seq peak calling best way usage of parameters
written by luckysardar171 •Hello, currently working on chip-seq analysis calling peaks using MACS2, need suggestions how to set parameters according to my data i have H3K4me1 H3K4me3 and …
-
Providing R1,R2 to macs2
written by AspireI have paired-end ChIP-Seq data that resulted in a low amount of uniquely mapped reads. When mapping the reads PE, only a small percentage maps …
-
Cleaning phantom peaks
written by chualec •From the ENCODE3 QC suggestions, phantom peaks give cross-correlation values that may help identify if your chip-seq worked or not, and also gives numbers based …
-
RNA Seq analysis of rice samples
written by alim.hcu •Dear all I have RNA-seq data from rice crop which is model and well annotated crop. When i mapped the reads, i found that most …
-
How to do the peakcalling for fly spike-in ChIPseq data
written by lhaiyan3Dear all, We have the fly spike-in ChIP-seq data, I use the following tab to calculate the correlation factor, File/Sample name DMSO_input Induced_Input DMSO_TF1 Induced_TF1 …
-
peaks overlap: ATAC & CHIP seq
written by sophialovechan •Hi guys. I have two peak files. One is from ATAC seq and the other is from H3K27ac chip seq. I am trying to merge …
-
How can I increase the number of peaks using macs2 peak caller?
written by GK1610I have a set of chip seq bam files. I have ~ 24 million paired-end reads in these bam files I used macs2 for narrow …
-
samtools flagstat result is different from samtools view result
written by Louis Kok •Hi All. I wonder if anyone of you have inconsistent mapping result generated by `samtools flagstat` and `samtools view` command? My samtools flagstat output is …
Yes, afaik only non-properly paired reads are removed. Hence, a MAPQ filter should be applied before peak calling, e.g. MAPQ > 10 or 20.
If the multi-mapped reads are part of your BAM or BED file that you feed into MACS2, it will use them.