Question: MACS2 won't build model - too few paired peaks, but I'm not using paired data
0
gravatar for bertb
23 months ago by
bertb10
bertb10 wrote:

Hello,

I am trying to run a macs2 peak analysis on bam files generated from unpaired reads. Format = AUTO, with no additional options invoked. The warning I receive is "Too few paired peaks (0) so I can not build the model! Broader your MFOLD range parameter may erase this error..."

Re running the script with --nomodel and a suggested --extsize provides strange data that does not line up with what I can see on IGV.

I suspect the issue is that macs is trying to build a model with paired reads when I haven't provided paired data - is this possible? Should I be adjusting my usage parameters?

Thanks,

chip-seq • 1.4k views
ADD COMMENTlink modified 23 months ago • written 23 months ago by bertb10

This has nothing to do with paired-end sequencing. The paired peaks are based on the shifting model that MACS builds based on tags aligned to the plus and minus strand. Please read the MACS paper. Can you post a screenshot of the genome browser? Maybe the data simply do not contain peaks (be it for a biological or technical reason). What was the target of the ChIP (which protein)?

ADD REPLYlink written 23 months ago by ATpoint40k

Thanks for the fast response - I'll try and get a screenshot up. Appreciate the help. (Originally posted as answer in error).

ADD REPLYlink written 23 months ago by bertb10

I'm not able to post a screenshot, but when I load the bam and bedgraph files in the genome browser, the bedgraph files basically look like a compressed version of the bam files. The background is higher, and peaks that appear more than 5-fold above background in the bam files are less than 2-fold enriched in the bedgraph files.

Thanks in advance for any suggestions to help generate a more accurate peakset.

ADD REPLYlink modified 23 months ago • written 23 months ago by bertb10

See this: How to add images to a Biostars post

ADD REPLYlink written 23 months ago by genomax91k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1049 users visited in the last hour