I am attempting to bin my BED file that contains gene intervals -2k/2k from the TSS. I would like to produce bins that are 1bp each.
After binning I would like to count anti-sense and sense reads independently. The output data i 'm looking for is:
GeneID -2000 -1999 -1998 -1997 ...... 0 ...... 1997 1998 1999 2000
SGIP 0.00 0.00 0.00 1 1 0.00 0.00 9 0.00
Where the first column is the GeneID provided, and the numbers after represent the number of read counts (anti-sense / sense) based on bin interval from TSS.
I'm not quite sure how to do something like this. I at first thought it might be as simple as using
bedtools makewindowsbut i'm unsure if this is appropriate. For read counts, I need a tool that is strand specific and capable of counting anti-sense and sense reads. Which tool would work best for this and what are the input files I require?