Chipseq data peak calling issue
Entering edit mode
9 months ago
Adyasha • 0

Hi ,

I'm trying to do analysis of chipseq data. I have 3 samples: Sample1, sample2 and input

I have done QC and then alignment using Bowtie. After that I used samtool to get bam files. Then I have used Picard for duplicate removal. Now I used macs2 for peak calling. I want to see transcription factor binding region.

Then I used

macs2 callpeak -t sample1.bam -c input.bam -f BAM -g genome size --outdir peak calling

But I couldn't get any .bed files from this. I only got .peak, .xlsx, summit.bam.

Then I tried annotating the peaks using this command <.peak(output from macs2)> <hg38>  <output file> 

But this gave me a text file, with some peak ID , starting and end position, but it doesn't give any gene name or ID.

My output doesn't look like this:

enter image description here

Some end columns are missing, and also I think I have very less peaks than expected. Can anyone please help me on this issue?

Thank you

ChIP-seq RNA sequencing • 347 views

Login before adding your answer.

Traffic: 1378 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6