Question: peak calling for medip seq analysis using macs
gravatar for subhoshree.ghose
10 weeks ago by
subhoshree.ghose0 wrote:

Dear all,

I am trying to use macs1.4 for peak calling with my alignment file which is in .sam format. This is a medip sequencing data. However I am facing an issue. 1. While I am using the trimmed and aligned file (.sam format) I am getting 258939 no of peaks.

  1. Afterwards I have converted the sam file to a bam file sorted them, filtered unmapped reads and again did run run macs on trimmed,aligned,sorted and filtered reads file (.sam format) and agin getting same no. of peaks i.e 258939

this is confusing. since the no of peaks are exactly same even after the filtering unmapped reads. Hope to get some response..

Thanking in advance,

next-gen • 169 views
ADD COMMENTlink written 10 weeks ago by subhoshree.ghose0
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: 880 users visited in the last hour