Macs2: predicted fragment length less than 2*tag length
0
1
Entering edit mode
5.6 years ago
langya ▴ 110

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?

ChIP-Seq macs2 tagsize • 2.7k views
ADD COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 2514 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