From c3d58e9215f20d8a5d9571370f899c499d896c0b Mon Sep 17 00:00:00 2001
From: NMarthe <nina.marthe@ird.fr>
Date: Fri, 20 Oct 2023 11:52:09 +0200
Subject: [PATCH] =?UTF-8?q?ajout=20d'une=20fct=20pour=20transcrire=20les?=
=?UTF-8?q?=20s=C3=A9quences=20d'adn?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
Functions.py | 18 ++++++++++++++++--
1 file changed, 16 insertions(+), 2 deletions(-)
diff --git a/Functions.py b/Functions.py
index 7e35d35..fd28c3d 100644
--- a/Functions.py
+++ b/Functions.py
@@ -497,13 +497,27 @@ def add_feature_sequence(feature,seg_seq):
feature_sequence=""
for segment in feature.segments_list:
if segment==feature.segments_list[0]:
- feature_sequence+=get_segment_sequence(seg_seq,segment)[feature.pos_start:] # vérifier les +/- 1 pour la position
+ feature_sequence+=get_segment_sequence(seg_seq,segment)[feature.pos_start-1:] # revérifier les +/- 1 pour la position, avec de vraies données
elif segment==feature.segments_list[-1]:
- feature_sequence+=get_segment_sequence(seg_seq,segment)[0:feature.pos_stop] # vérifier les +/- 1 pour la position
+ feature_sequence+=get_segment_sequence(seg_seq,segment)[0:feature.pos_stop] # revérifier les +/- 1 pour la position, avec de vraies données
else:
feature_sequence+=get_segment_sequence(seg_seq,segment)
feature.sequence=feature_sequence
+def transcription(dna_sequence):
+ rna_sequence=""
+ for dna_base in dna_sequence:
+ match dna_base:
+ case "A":
+ rna_sequence+="U"
+ case "C":
+ rna_sequence+="G"
+ case "G":
+ rna_sequence+="C"
+ case "T":
+ rna_sequence+="A"
+ return rna_sequence
+
# penser à transcrire la séquence codante du gène !!
def get_aa(codon):
match codon[0:2]:
--
GitLab