Macs2: predicted fragment length less than 2*tag length
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?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
I do not understand how to get total peaks from MACS2 shown to me. What mistake(s) am I making here?
written by Ronin •
-
MACS2: broad and broad no-model output peaks are exactly the same. Is something wrong?
written by piyushjo
-
Finding Fragment Size chip-seq paired end reads
written by luckysardar171 •
-
MACS2- struct.error: unpack requires a string argument of length 4
written by 200dumplings •
-
Peak calling using macs2?
written by star
-
MACS2 MODEL VS NO MODEL FOR CHIRP SEQ
written by prathyushakonda •
-
macs2 warnings unknown sequencing problem
written by GR
-
how to set a proper d parameter in MACS
written by Lila M
-
MACS2 Nomodel doesnt take the fragment size
written by GK1610
-
criteria for choosing proper parameters for MASC2 -- chip-seq ?
written by jimmy_zeng
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