from pydna.genbank import genbank
lacZ = genbank("NC_000913 REGION: complement(363231..366305)")
from pygenome import saccharomyces_cerevisiae as sg
GAL4 = sg.stdgenes["GAL4"].cds()
PPR1 = sg.stdgenes["PPR1"].cds()
GPD1 = sg.stdgenes["TDH3"].cds()
lacZ.express("sce")
cds | len | cai | gc | sta | stp | n-end | CGA | CGG | CGC | CCG | CTC | GCG | rare |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ATG...TAA | 1025.0 | 0.552 | 0.563 | 1.0 | 0.47 | >30 h | 3 | 7 | 37 | 36 | 9 | 28 | 0.117 |
GAL4.express("sce")
cds | len | cai | gc | sta | stp | n-end | CGA | CGG | CGC | CCG | CTC | GCG | rare |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ATG...TAA | 882.0 | 0.701 | 0.408 | 1.0 | 0.47 | 3 min | 7 | 0 | 3 | 11 | 8 | 6 | 0.04 |
PPR1.express("sce")
cds | len | cai | gc | sta | stp | n-end | CGA | CGG | CGC | CCG | CTC | GCG | rare |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ATG...TAG | 905.0 | 0.718 | 0.386 | 1.0 | 0.23 | 3 min | 7 | 2 | 5 | 4 | 3 | 9 | 0.033 |
x = GPD1.express("sce")
x
cds | len | cai | gc | sta | stp | n-end | CGA | CGG | CGC | CCG | CTC | GCG | rare |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ATG...TAA | 333.0 | 0.814 | 0.47 | 1.0 | 0.47 | >30 h | 0 | 0 | 0 | 0 | 0 | 0 | 0.0 |
x.lol()
[['cds', 'len', 'cai', 'gc', 'sta', 'stp', 'n-end', 'CGA', 'CGG', 'CGC', 'CCG', 'CTC', 'GCG', 'rare'], ['ATG...TAA', 333.0, 0.814, 0.47, 1.0, 0.47, '>30 h', 0, 0, 0, 0, 0, 0, 0.0]]