CHIP-SEQ - multiply aligned reads handling

Do I need to filter out reads that have been aligned more than once before peakcalling with MACS2?

chip-seq

• 1,100 views • link

updated 7.2 years ago by Biostar

• written 7.4 years ago by kozaleo •

https://groups.google.com/forum/#!topic/macs-announcement/ASQ3GsXQC1c

Depends on the data you are trying to process. If your regions of interest are not repeated you can filter out them before peakcalling

• 0 views • link

written 7.4 years ago by Bastien Hervé

0 answers

No answers yet.

More posts like this

Barley, alignment rates
written 1.5 years ago by IrK

Hi everyone, I'm working with single-end reads (150 bp) from barley (Hordeum vulgare) and aligning them against the Morex V3 reference genome. I have a …
Where to check out the alignment results for BWA
written 5.8 years ago by field654 •

I used to use Bowtie2 and it gives the results "Out of xxxxxxxxx reads, yyyyyyyy have been aligned 0 times, zzzzzz have been aligned exactly …
How can I filter RNA-seq reads to remove reads where more than 50% of the read is low quality (Q<20…
written 6.2 years ago by mar77 •

I'm trying to filter out any reads that consist of more than 50% low quality reads (Q<20). I've already performed filtering steps to remove adapter …
Overlap in Peaks of ChIP data for 2 TFs
written 7.3 years ago by yugen30 •

Hi all, I have recently started analyzing ChIP-seq data. I have two datasets from GEO for two different transcription factors in same sample and want …
How to know the mm9 genome size?
written 8.2 years ago by mikysyc2016

When I use MACS2 to do peakcalling, i need to know genome size(gsize), and I use mm9. How I can know that? Thanks!
Uniquly aligned reads by BWA-mem (MAPQ, and XA tag)
written 8.7 years ago by biobiu

I'm trying to filter out reads that were primarily aligned to one genomic region (if a read has secondary alignment I want to keep the …
PeakAnalyzer & PeakSplitter
written 9.2 years ago by Brice222 •

Greetings Please help me... I did peakcalling with MACS2. Output MACS2 includes not wig file. For analysis by PeakAnalyzer in part PeakSplitter of which I …
file "_peaks.bed" on macs2
written 9.2 years ago by Brice222 •

Hi I did PeakCalling with Macs2. I need to file format "_peaks.bed". But in the new version, does not exist the file format "_peaks.bed" . …
Filtering high quality reads from bowtie2
written 11.0 years ago by yulia.gordon •

Hi everyone! I've encountered an interesting discrepancy when trying to filter the results of my paired end CHIP-SEQ. When using bowtie2 in Galaxy interface it …
Filter fastq/sam/bam for reads
written 12.1 years ago by hlsz.laszlo •

Dear All, I'm analyzing a ChIP-seq data, and I having some trouble filtering out "good" reads for us. Briefly, I've got a fastq file, then …