Question: Macs2: predicted fragment length less than 2*tag length
gravatar for langya
7 months ago by
langya40 wrote:

I am using

macs2 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 • 455 views
ADD COMMENTlink modified 5 months ago by Biostar ♦♦ 20 • written 7 months ago by langya40

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 7 months ago • written 7 months ago by ATpoint15k
Please log in to add an answer.


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