This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English version of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materials including this page, select Japan from the country navigator on the bottom of this page.


Find restriction enzymes that cut nucleotide sequence


[Enzymes, Sites] = rebasecuts(SeqNT)
rebasecuts(SeqNT, Group)
rebasecuts(SeqNT, [Q, R])
rebasecuts(SeqNT, S)

Input Arguments

SeqNTNucleotide sequence.
GroupCell array of character vectors representing the names of valid restriction enzymes.
Q, RBase positions that limit the search to all sites between base Q and base R.
SBase position that limits the search to all sites after base S.

Output Arguments

EnzymesCell array with the names of restriction enzymes from REBASE®, the Restriction Enzyme Database.
SitesVector of cut sites identified with the base position number before every cut.


[Enzymes, Sites] = rebasecuts(SeqNT) finds all the restriction enzymes that cut SeqNT, a nucleotide sequence.

rebasecuts(SeqNT, Group) limits the search to Group, a list of enzymes.

rebasecuts(SeqNT, [Q, R]) limits the search to those enzymes that cut after the base position specified by Q and before the base position specified by R.

rebasecuts(SeqNT, S) limits the search to those enzymes that cut just after the base position specified by S.

REBASE, the Restriction Enzyme Database, is a collection of information about restriction enzymes and related proteins. For more information about REBASE, see:


  1. Create a nucleotide sequence.

  2. Find all possible enzymes and cleavage sites in the sequence.

     [enzymes, sites] = rebasecuts(seq)
  3. Find where restriction enzymes CfoI and Tru9I cut the sequence.

    [enzymes, sites] = rebasecuts(seq, {'CfoI','Tru9I'})
    enzymes = 
    sites =
  4. Find all possible enzymes that cut after base 7.

    enzymes  = rebasecuts(seq, 7)
    enzymes = 
  5. Find all possible enzymes that cut between bases 11 and 37.

    enzymes  = rebasecuts(seq, [11 37])
    enzymes = 


[1] Roberts, R.J., Vincze, T., Posfai, J., and Macelis, D. (2007). REBASE—enzymes and genes for DNA restriction and modification. Nucl. Acids Res. 35, D269–D270.

[2] Official REBASE Web site:

Introduced before R2006a

Was this topic helpful?