Paired End RIP-seq analysis with MACS2 reasonable parameters?
0
0
Entering edit mode
21 months ago
ddierks91 • 0

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

RNA-Seq RIP-Seq MACS2 MACS Peak Caller • 766 views
ADD COMMENT
0
Entering edit mode

What do you mean with has issues with the paired-end calculated insert/fragment sizes.? It has a -f BAMPE option that simply reads TLEN. Please explain.

ADD REPLY

Login before adding your answer.

Traffic: 1903 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6