Question: Does strand information still exist after reads merging?
gravatar for Shaway
3.1 years ago by
Shaway10 wrote:

hello everyone! I encountered a question and do not know how to tackle with it. I have some reads(using bedtools bamtobed function) like this

chr10   29334602    29334652    SRR891270.293/2 255 +

chr10   29334606    29334656    SRR891270.293/1 255 -

chr10   61751512    61751562    SRR891270.256/2 255 +

chr10   61751525    61751575    SRR891270.256/1 255 -

chr10   103211892   103211942   SRR891270.153/2 255 +

chr10   103212015   103212065   SRR891270.153/1 255 -

and I want to merge these reads and using the output to call peaks by macs2. I take the fisrt number in + strand reads and second number in -strand reads, the output like below:

chr10   29334602    29334656    SRR891270.293    255

chr10   61751512    61751575    SRR891270.256    255

chr10   103211892   103212065   SRR891270.153    255

but there are no strand information!

in macs2, the 6 column is needed(For BED format, the 6th column of strand information is required by MACS. And please pay attention that the coordinates in BED format is zero-based and half-open (

what can I do to solve this problem?


chip-seq • 788 views
ADD COMMENTlink modified 3.1 years ago by Devon Ryan96k • written 3.1 years ago by Shaway10
gravatar for Devon Ryan
3.1 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

After merging, the concept of strand (actually, it's just orientation) is meaningless. You appear to have a BAM file originally, so just use that as the input to MACS2.

ADD COMMENTlink written 3.1 years ago by Devon Ryan96k

thanks for you answer.

ADD REPLYlink written 3.1 years ago by Shaway10
Please log in to add an answer.


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