Cut Sites Per Megabase

[kevindoh]

Hey guys,

I can't seem to find the right tool to do this. I was wondering if anyone had advice on how to approach.

I am attempting to get the number of cut sites per megabase for the whole genome (given a recognition site). I'd rather not write all the code from the bottom up. Does anyone have a tool that would do this efficiently? Or a strategy for coding it myself that doesn't involve days of work?

Thanks in advance

K

[rhinoceros]

If it's a fasta file,

grep -v ">" yourFile.fasta | grep -c "ATGCCCGC"

Just replace the seq with whatever..

edit. This tells the total number. You'd obviously have to divide with seq length to get per whatever..

[lindenb]

use Biostar86480 https://github.com/lindenb/jvarkit/wiki/Biostar86480 and pipe the output in bedtools/groupby http://bedtools.readthedocs.org/en/l...s/groupby.html