MACS2 transcription start site analysis of .BAM files
1
0
Entering edit mode
2.2 years ago

Hi! I'm new to using MACS2 and I would love some help :)

I have a .BAM file from an AF9 ChIP-Seq done with mouse (mm9 genome) cells. I'm trying to find binding sites -500 to +1kb from transcription start sites but I'm not sure how to do this.

Do I need to provide a second file with TSS locations for mm9? How do I even set the parameters for my command to have a start and end point? Is MACS2 even the best software for this task (maybe Homer would better)?

Thank you so much!

ChIP-Seq MACS2 TSS • 611 views
ADD COMMENT
2
Entering edit mode
2.2 years ago
ATpoint 55k

Run macs as usual to obtain a peak file (narrowPeak format). Make a BED file with your TSS coordinates. Use bedtools intersect to extract peaks overlapping the TSS coordinates. By the way, if you have no legacy burden change to mm10, mm9 is very old and there is no reason to use it unless your project relies on legacy / previous data on mm9.

ADD COMMENT
0
Entering edit mode

Thank you! Unfortunately I had to duplicate some research to get my .bam file that used mm9 :/

ADD REPLY

Login before adding your answer.

Traffic: 2169 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