Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Megabase 1000 paakhan Introductions 3 12-01-2016 09:10 AM
How to choose a FPKM cut-off wupengpro Bioinformatics 10 05-16-2013 05:19 PM
cut fasta spike1985 General 1 02-14-2012 08:16 AM
Tool for Identifying Restriction Sites at IUPAC SNP sites? brachysclereid Bioinformatics 5 02-09-2011 01:58 AM

Thread Tools
Old 01-18-2014, 06:55 AM   #1
Junior Member
Location: Ann Arbor

Join Date: Jan 2014
Posts: 1
Default Cut Sites Per Megabase

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

kevindoh is offline   Reply With Quote
Old 01-18-2014, 07:20 AM   #2
Senior Member
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372

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..

Last edited by rhinoceros; 01-18-2014 at 07:23 AM.
rhinoceros is offline   Reply With Quote
Old 01-18-2014, 08:04 AM   #3
Senior Member
Location: France

Join Date: Apr 2010
Posts: 143

use Biostar86480 and pipe the output in bedtools/groupby
lindenb is offline   Reply With Quote

batch entrez, programming, restriction enzyme

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:15 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO