Question: Macs2: predicted fragment length less than 2*tag length
0
gravatar for langya
14 months ago by
langya50
langya50 wrote:

I am using

macs2 2.1.1.20160309 to call peaks on narrow peak ChIP-Seq data set that is sequenced with Single Read 101 bp and then trimmed the reads with TRIMMOMATIC and mapped to bowtie2 (with good mappability), I then called Peaks with macs2 with following command: macs2 callpeak -t $1 -c $2 -f BAM -g hs -n $1 -B --SPMR -q 0.05

But then in the macs2 output, there is this warning:

2 predicted fragment length is 93 bps

2 alternative fragment length(s) may be 93,463,534 bps

2 Since the d (93) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem!

Could anyone help me out?

tagsize chip-seq macs2 • 748 views
ADD COMMENTlink modified 12 months ago by Biostar ♦♦ 20 • written 14 months ago by langya50

Could be a poor ChIP efficiency because of a bad antibody, but that is guessing. Can you post a screenshot of the IGV (or whatever browser) for a representative region, like a 10kb window or so, centering on what you would call a peak? How to add images on biostars

ADD REPLYlink modified 14 months ago • written 14 months ago by ATpoint25k
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: 1471 users visited in the last hour