Did I find the right examples for you? yes no      Crawl my project      Python Jobs

All Samples(2)  |  Call(1)  |  Derive(0)  |  Import(1)

src/j/c/jcvi-HEAD/apps/restriction.py   jcvi(Download)
import logging
 
from Bio.Restriction.Restriction import AllEnzymes, Analysis
 
from jcvi.formats.fasta import Fasta, SeqRecord, SeqIO
    fw = open(fragfastafile, "w")
    for name, rec in f.iteritems_ordered():
        a = Analysis([enzyme], rec.seq)
        sites = a.full()[enzyme]
        extract(rec, sites, flank, fw)