Chipseq data peak calling issue
0
0
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

.annotatePeaks.pl <.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
ADD COMMENT

Login before adding your answer.

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

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

Powered by the version 2.3.6