SEQanswers (
-   Bioinformatics (
-   -   _peaks.bed (

khb 02-07-2011 12:48 AM

Missing strand _peaks.bed
I've tried to analyze my data with MACS, and want the file _peaks.bed that MACS gives. Does anyone know how to know which strand each peak is on? The file only contains
chrom start stop name score

ffinkernagel 02-07-2011 01:46 AM

The answer to your question is 'mu'.

A chipseq peak has no strandedness... in theory, you should only see reads of one strand on one, and of the other strand on the other side, crossing over wherever your protein bound the chromatin (since fragments starting further down on the same strand would not have a protein bound, and would not be visible after the chip). In practice, it's not that clear cut. Nevertheless, the peak does not have a strand.

khb 02-07-2011 02:58 AM

Of course. Thanks

rossignol 08-14-2011 10:40 AM


I see the point about peaks not having a strand. What should I then submit to HOMER when doing motif analysis? HOMER expects to get a peak file in bed format, with the last column being +/- for the strand of the peak.

All times are GMT -8. The time now is 03:05 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.