In [1]:
from pydna.all import *
pYPK0 =read("pYPK0.gb")
gb=Genbank("bjornjobb@gmail.com")
template = gb.nucleotide("AF298793.1") # pUG6
#template = gb.nucleotide("AJ002680.1") # pFA6a-kanMX4
In [2]:
p678, p666 =parse('''
>678_pYPK0_hygfwd: (77-mer)
ctcacgttaagggattttggtcatgagCACATACGATTTAGGTGACACTATAGAAC
>666_pYPK0_hygrev (70-mer)
catctttgacagcttatcatcgataagctCGACTCACTATAGGGAGACC''', ds=False)
kanmarker =pcr( p678, p666, template)
In [3]:
from Bio.Restriction import SbfI
vect = pYPK0.linearize(SbfI)
In [4]:
asm=Assembly([kanmarker, vect])
asm
Out[4]:
In [5]:
seq = asm.assemble_circular()[0]
pYPK7 = seq.synced(pYPK0)
In [6]:
pYPK7.stamp()
Out[6]:
In [7]:
pYPK7.locus="pYPK7"
In [8]:
pYPK7.write("pYPK7.gb")
In [9]:
from pydna.all import *
reloaded=read("pYPK7.gb")
assert reloaded.cseguid() in reloaded.definition
In [10]: