Hi everyone,
I have a GRO-Seq NGS dataset (single-end). I would like to plot the average coverage around the TSS (+/- 10Kb) for all the genes in the genome, in much the same way as it is done in the plot on the left of the attached figure (which can also be found at http://www.ncbi.nlm.nih.gov/pmc/arti...333/figure/F2/).
I know there are good options for TSS profile plots such as HTSeq of ngsplot, but since they are thought for ChIP-Seq experiments they merge both forward and reverse strands. Due to the nature of the GRO-Seq experiment I need to keep both strands independently.
I would be grateful if someone could point me to the best way to do this type of plot.
Thank you so much in advance!
Xavi.
I have a GRO-Seq NGS dataset (single-end). I would like to plot the average coverage around the TSS (+/- 10Kb) for all the genes in the genome, in much the same way as it is done in the plot on the left of the attached figure (which can also be found at http://www.ncbi.nlm.nih.gov/pmc/arti...333/figure/F2/).
I know there are good options for TSS profile plots such as HTSeq of ngsplot, but since they are thought for ChIP-Seq experiments they merge both forward and reverse strands. Due to the nature of the GRO-Seq experiment I need to keep both strands independently.
I would be grateful if someone could point me to the best way to do this type of plot.
Thank you so much in advance!
Xavi.