Question: macs2 callpeaks genome
gravatar for sophialovechan
2.9 years ago by
sophialovechan40 wrote:

Hi guys, I have a question about MACS2 to call peaks from either chip-sew or ATAC seq. When I used macs2 to call peaks, the program used hs for human genome, which corresponds to hg18. But in later analysis, when I useD homer to find motifs in peaks, I used hg19. Is there any influence since the position of genes in different genome project is different? Thanks. Here is the scripts: /usr/local/macs2/latest/bin/macs2 callpeak --outdir -t * -c * -n * -g hs -B --keep-dup 1 --nomodel --shift -75 --extsize 150 -q 0.01 /home/users/homer/bin/ *.txt hg19 *_discover/ -size -100,100 -mask -S 50

chip-seq • 2.4k views
ADD COMMENTlink modified 2.9 years ago by Sinji2.8k • written 2.9 years ago by sophialovechan40
gravatar for Sinji
2.9 years ago by
UT Southwestern Medical Center
Sinji2.8k wrote:

The difference is minimal at best, at least in the peak calling. I remember reading a post by Istvan once that said that as long as the first number of the effective genome size is the same (ex. 2.7 billion in hg18, 2.8 billion in hg19) then you should be fine.

ADD COMMENTlink written 2.9 years ago by Sinji2.8k
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: 887 users visited in the last hour