Paired End RIP-seq analysis with MACS2 reasonable parameters?
Since some recent publications involving MACS2 as a valuable RIP-Seq (RNA not DNA) peak caller. I tried it to use on my data. But I realize that it since its a Chip-Seq peak-caller, has issues with the paired-end calculated insert/fragment sizes. I work with poly-A selected Material (mostly mRNA lacking Introns).
So I wonder if there are any parameters where I can use actual Fragment Sizes not in a genomic setting (but based on transcript annotations)?.
Otherwise, I don't see how MACS2 can be a reasonable paired-end Peak Caller for RIP-seq analysis.
Greetings and Thanks! David
• 2,321 views
•
link
written
by
ddierks91 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Optimal MACS2/3 configuration for paired-end ChIP-seq data: `-f BAM` versus `-f BAMPE`
written by kalavattam •
-
Does the bigwig file of a single-end chip-seq shows 2 peaks for a fragment whereas the paired-end s…
written by Rajendra KC •
-
CUT&Tag data processing: peak caller choice and typical peak size profile
written by nihilior •
-
BAM reads much longer than fastq?
written by Jay •
-
How can I filter and detect circular RNA aligned paired-end reads with star?
written by ddierks91 •
-
Macs2: predicted fragment length less than 2*tag length
written by langya •
-
Summit Extension in ATAC-seq
written by shankasal •
-
MACS2 gives negative estimated fragment length
written by dcheng1 •
-
Differential chip-seq peak calling using macs2
written by lucas.michel •
-
Why estimate fragment size of ChIP-seq with SPP when using Macs2 for peak calling?
written by colin.kern
What do you mean with
has issues with the paired-end calculated insert/fragment sizes.? It has a-f BAMPEoption that simply reads TLEN. Please explain.