This vector is identical to pMEC1126 (pYPK0_TEF1_PsXYL1_TDH3tp) but for a N272D mutation in the XYL1 reading frame.
In [1]:
    
from pydna.all import *
    
In [2]:
    
pYPK0_TEF1_PsXYL1_TDH3 =read("pYPK0_TEF1_PsXYL1_TDH3.gb")
    
In [3]:
    
len(pYPK0_TEF1_PsXYL1_TDH3)
    
    Out[3]:
In [4]:
    
# NBVAL_IGNORE_OUTPUT
from IPython.display import IFrame
IFrame('http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2988607', width='100%', height=250)
    
    Out[4]:
In [5]:
    
pYPK0_TEF1_PsXYL1_TDH3.list_features()
    
    Out[5]:
In [6]:
    
xyl1 = pYPK0_TEF1_PsXYL1_TDH3.extract_feature(5)
    
In [7]:
    
wt_protein_seq = str(xyl1.seq.translate())
    
In [8]:
    
print(wt_protein_seq)
    
    
In [9]:
    
print(wt_protein_seq[271])
    
    
In [10]:
    
pYPK0_TEF1_PsXYL1_TDH3.features[6].location.start
    
    Out[10]:
In [11]:
    
codon = pYPK0_TEF1_PsXYL1_TDH3[1163+271*3:1163+271*3+3]
    
In [12]:
    
print(codon.seq)
print(codon.seq.translate())
    
    
In [13]:
    
candidate = (pYPK0_TEF1_PsXYL1_TDH3[:1163+271*3]+
             Dseqrecord("GAC")+
             pYPK0_TEF1_PsXYL1_TDH3[1163+271*3+3:]).looped()
    
In [14]:
    
new_feature = pYPK0_TEF1_PsXYL1_TDH3.features[8]
new_feature.qualifiers = {"label":"PsXYL1_N272D"}
candidate.features.append(new_feature)
    
In [15]:
    
candidate
    
    Out[15]:
In [16]:
    
candidate.cseguid()
    
    Out[16]:
In [17]:
    
candidate.id = "pMEC1127"
#candidate.description = "pYPK0_TEF1_PsXYL1_N272D_TDH3tp"
candidate.stamp()
    
    Out[17]:
In [18]:
    
candidate.write("pYPK0_TEF1_PsXYL1_N272D_TDH3.gb")
    
    
In [19]:
    
reloaded = read("pYPK0_TEF1_PsXYL1_N272D_TDH3.gb")
    
In [20]:
    
reloaded.cseguid()
    
    Out[20]: