Biostar Beta. Not for public use.
Question: Binning BED file and counting antisense / sense reads
0
Entering edit mode

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?

ADD COMMENTlink 4.1 years ago Sam • 70 • updated 4.1 years ago Amitm ♦ 1.6k
1
Entering edit mode

hi,

You could skip the step of creating 1bp bins. Just use one line for each 2k window. And then using Bedtools Coverage use parameter _-d. That will give depth at each base anyway. Also you can use the -s_ to force strandedness. Obviously your input BED must have the strand column.

ADD COMMENTlink 4.1 years ago Amitm ♦ 1.6k
Entering edit mode
0

This is great to hear! I'll sort this out today and let you know how it goes so I can accept your answer or ask additional questions.

ADD REPLYlink 4.1 years ago
Sam
• 70

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0