From 5b8592fbd75de93390ca79ec9de3b28b0e4b9c53 Mon Sep 17 00:00:00 2001 From: mpabi Date: Sun, 9 Jun 2024 17:39:24 +0200 Subject: [PATCH] u --- gfp-old.py | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/gfp-old.py b/gfp-old.py index e759c73..999a65e 100644 --- a/gfp-old.py +++ b/gfp-old.py @@ -31,6 +31,8 @@ def load_local_sequence(file_path): record = SeqIO.read(file, "genbank") return record +#%% + # Ścieżki do lokalnych plików gfp_file_path = 'gfp_sequence.gb' plasmid_file_path = 'plasmid_sequence.gb' @@ -65,9 +67,11 @@ else: # Dodanie sekwencji kodującej His6 his6_tag = "CACCACCACCACCACCAC" # Sekwencja kodująca 6x histydynę (His6) +#%% # Sekwencja białka GFP z His6 tagiem na końcu -gfp_seq_with_his6 = gfp_seq + his6_tag +t = gfp_seq_with_his6 = gfp_seq + his6_tag +#%% # Projektowanie starterów forward_primer = gfp_seq[:20] # Pierwsze 20 nukleotydów sekwencji GFP reverse_primer = str(Seq(gfp_seq_with_his6[-20:]).reverse_complement()) # Ostatnie 20 nukleotydów + His6 @@ -75,10 +79,13 @@ reverse_primer = str(Seq(gfp_seq_with_his6[-20:]).reverse_complement()) # Ostat print(f"Forward Primer: {forward_primer}") print(f"Reverse Primer: {reverse_primer}") + +#%% # Miejsca cięcia restryktazy enzymes = RestrictionBatch([EcoRI, BamHI, HindIII]) restriction_sites = enzymes.search(record.seq) +#%% # Utwórz diagram plazmidu diagram = GenomeDiagram.Diagram("pET-28a(+) Plasmid Map") track = diagram.new_track(1, name="Annotated Features", greytrack=True)