pCAPs_EcoRV_ScENO2tp = (pYPKa_EcoRV + ScENO2tp).looped() ScRPE1 = pcr( primer[556], primer[557], Dseqrecord(sgd.cds("RPE1"))) pCAPs_AjiI_ScRPE1 = (pYPKa_AjiI + ScRPE1).looped() ScRKI1 = pcr( primer[555], primer[554], Dseqrecord(sgd.cds("RKI1"))) pCAPs_AjiI_ScRKI1 = (pYPKa_AjiI + ScRKI1).looped() CiGXF1 = pcr( primer[553], primer[552], Dseqrecord(gb("AJ937350"))) pCAPs_AjiI_CiGXF1 = (pYPKa_AjiI + CiGXF1).looped() from pYPK0 import pYPK0 as pCAPs_pSU0 pCAPs_pSU0_E_Z, stuffer = pCAPs_pSU0.cut((EcoRV,ZraI)) ypk_primer = parse("/home/bjorn/Dropbox/wikidata/ypk_primers.wiki") #ScPDC1_PacI-ScTKL1-ScRPS19b_PacI A_ScPDC1tp_b = pcr( primer[167], primer[567], pCAPs_ZraI_ScPDC1tp) ScTKL1 = pcr( ypk_primer[8], ypk_primer[9], Dseqrecord(sgd.cds("TKL1"))) C_ScRPS19btp_d = pcr( primer[568], primer[166], pCAPs_EcoRV_ScRPS19btp) circular = Assembly(( A_ScPDC1tp_b, ScTKL1, C_ScRPS19btp_d,
# RPS19atp_ScRKI1tp_TPI1tp # TPI1tp_CiGXF1_ENO2tp from pYPK0_PDC1tp_ScTKL1_RPS19btp import pYPK0_PDC1tp_ScTKL1_RPS19btp as cas1 from pYPK0_RPS19btp_ScRPE1_RPS19atp import pYPK0_RPS19btp_ScRPE1_RPS19atp as cas2 from pYPK0_RPS19atp_ScRKI1_TPI1tp import pYPK0_RPS19atp_ScRKI1_TPI1tp as cas3 from pYPK0_TPI1tp_CiGXF1_ENO2tp import pYPK0_TPI1tp_CiGXF1_ENO2tp as cas4 cas1 = pcr( p167, p647, cas1) cas2 = pcr( p648, p651, cas2) cas3 = pcr( p538, p621, cas3) cas4 = pcr( p419, p166, cas4) pYPK0_E_Z, stuffer = pYPK0.cut((EcoRV, ZraI)) #print [cas1, cas2, cas3, cas4, pYPK0_E_Z] a = Assembly([cas1, cas2, cas3, cas4, pYPK0_E_Z], limit = 61) seq = a.circular_products[0] print seq.small_fig() seq=seq.synced("tcgcgcgtttcggtgatgacggtgaaaacctctg") pYPK0_T=seq ''' # coding: utf-8