Question: How to get strand information in the bed file of macs1.4 peak calling output?
0
gravatar for harikrishnanc903
12 months ago by
harikrishnanc9030 wrote:

I am new to MACS peak calling. I am not sure how to get the strand information in the bed output file after running macs tool. Using the strand information I need to extract the sequence from the genome. Any suggestion to add the strand information in the MACS output?

chip-seq next-gen • 648 views
ADD COMMENTlink modified 10 months ago by ZZzzzzhong200 • written 12 months ago by harikrishnanc9030

ChIPSeq is not strand specific hence that information is not available in the MACS2 output. Please see @Devon Ryan's answer on this thread Why chip-seq is strand specific?

ADD REPLYlink written 12 months ago by Sej Modha4.3k

Thanks. I have one more doubt. Homer does peak calling using the "findpeak".In homer output, we are getting strand information. Is it because homer uses a different algorithm to call the peak?

ADD REPLYlink written 12 months ago by harikrishnanc9030
1

As said, there is no strand information from the biological side on ChIP-seq data. Ignore the strand, wherever this tools gets it from. Probably only counts reads from one strand or whatever.

ADD REPLYlink written 12 months ago by ATpoint21k
0
gravatar for ATpoint
12 months ago by
ATpoint21k
Germany
ATpoint21k wrote:

There is no strand information, as ChIP-seq binding events are not strand-specific. DNA-binding factors and DNA-associated proteins, such as histones bind/are associated with double-stranded DNA. Why do you need the strand information, so what is your final goal?

ADD COMMENTlink modified 12 months ago • written 12 months ago by ATpoint21k

From macs output(bed file), I have contig start and end position. I will use this information to extract the sequence from the genome. I am concerned about whether we need the stand information in the bed file for doing the same.

ADD REPLYlink written 12 months ago by harikrishnanc9030

You don't. Unless you have a good reason (getting a gene sequence), sequences are always extracted from the reference (= plus) strand.

ADD REPLYlink written 12 months ago by ATpoint21k
0
gravatar for ZZzzzzhong
10 months ago by
ZZzzzzhong200
ZZzzzzhong200 wrote:

Transcriptome usually comes from fixed chains.peak_result is your MACS out put and genes.gtf is your annotation file.

intersectBed -wo -a peak_result -b genes.gtf | awk -v OFS="\t" '{print $1,$2,$3,"*","*",$10}'|uniq > Peaks.bed
fastaFromBed -s -f genome.fa -bed Peaks.bed -fo Peaks.fa
ADD COMMENTlink written 10 months ago by ZZzzzzhong200
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: 1599 users visited in the last hour