How to obtain sense and antisense read counts from Gro-Seq data
I am looking to obtain read-counts for sense and antisense strands from a GRO-Seq bedgraph file -2000/+2000 from the TSS for a set of genes in 1bp windows. Currently I am using HOMER with the -strand +/- -hist 1 -ghist option, is this appropriate?
Would htseq-count or bedtools coverage be better for this? I'm trying to generate a txt file with the following structure:
- 1st column contain the gene ID
- columns 2 to 4001 contain the number of Sense or Antisense reads respectively per nucleotide surrounding the TSS (+/- 2000bp)
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Ideas to split peaks into subpeaks - beyond PeakSplitter
written by Sergio Martínez Cuesta
-
Forum: Visualization of ChIP-seq data using Heatmaps (Updated: 06/10/16)
written by Sinji
-
Binning BED file and counting antisense / sense reads
written by Sam
-
Convergent Transcription Heatmap
written by Sam
-
Same data, very different heatmaps
written by Sam
-
Overlapping regions in a BedGraph file
written by Sam
-
Antisense transcription - how to detect it?
written by John
