Skip to content
Snippets Groups Projects
Commit c3d58e92 authored by nina.marthe_ird.fr's avatar nina.marthe_ird.fr
Browse files

ajout d'une fct pour transcrire les séquences d'adn

parent fe97d997
No related branches found
No related tags found
No related merge requests found
...@@ -497,13 +497,27 @@ def add_feature_sequence(feature,seg_seq): ...@@ -497,13 +497,27 @@ def add_feature_sequence(feature,seg_seq):
feature_sequence="" feature_sequence=""
for segment in feature.segments_list: for segment in feature.segments_list:
if segment==feature.segments_list[0]: 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]: 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: else:
feature_sequence+=get_segment_sequence(seg_seq,segment) feature_sequence+=get_segment_sequence(seg_seq,segment)
feature.sequence=feature_sequence 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 !! # penser à transcrire la séquence codante du gène !!
def get_aa(codon): def get_aa(codon):
match codon[0:2]: match codon[0:2]:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment