text: user8543: I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between them and which ones are different. What is the best way to do this? which package should I use?
  Many thanks
  
  user4766: 
  For sharing, you can test the proportion of overlapping peaks using bedtools intersect or the find.overlap function in GenomicRanges packages in R.
  What makes more sense is perhaps to merge all the peaks together using bedtools merge, and then count the number of reads/fragments in each peak for each sample using featureCounts. The output of this is a matrix, which is very similar to gene expression results. Then you can do some differential tests on this data.
  
  
original_text: 8543: <p>I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between them and which ones are different. What is the best way to do this? which package should I use?</p>
  
  <p>Many thanks</p>
   4766: <ol>
  <li>For sharing, you can test the proportion of overlapping peaks using <code>bedtools intersect</code> or the <code>find.overlap</code> function in <code>GenomicRanges</code> packages in R.</li>
  <li>What makes more sense is perhaps to merge all the peaks together using <code>bedtools merge</code>, and then count the number of reads/fragments in each peak for each sample using <code>featureCounts</code>. The output of this is a matrix, which is very similar to gene expression results. Then you can do some differential tests on this data.</li>
  </ol>
  