Flip the heatmaps of ChIP-seq data to match TSS of +/- strand in deeptools

Question

[Sho]

Hello everyone,

I'm now making heatmaps from ChIP-seq data to see the tendency of where transcription factors(TFs) bind. I have a list of promoter regions(10bp to 2000bp) including chr, start and end. (start < end)

When I plotted ChIP-seq data into these regions, the pattern of heatmaps looked symmetrical because + strand promoter's TSS is on the right side of the regions, and - strand promoter's TSS is on the other side.

So I'd like to flip the - strand promoter regions such that TSS of all the regions can be seen on the right side of the regions. I can separate the regions into + and - strand regions, but I don't know how I can merge and plot them.

Is there anyway I can do this?

Thanks for your help.

Answer 1

[Devon Ryan]

Just include the strand in the normal position (column 6) and this will be done automatically by deepTools.

[Sho]

It worked! Thank you!