Question: MACS2 transcription start site analysis of .BAM files
0
gravatar for carlee.bettler
13 months ago by
carlee.bettler10 wrote:

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!

tss chip-seq macs2 • 340 views
ADD COMMENTlink modified 13 months ago by ATpoint38k • written 13 months ago by carlee.bettler10
2
gravatar for ATpoint
13 months ago by
ATpoint38k
Germany
ATpoint38k wrote:

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 COMMENTlink modified 13 months ago • written 13 months ago by ATpoint38k

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

ADD REPLYlink written 13 months ago by carlee.bettler10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1285 users visited in the last hour