Bonjour,
après vérifié manuellement qu’il n’existait pas de caractère de terminaison de la chaine (pas de #), voici ma version :
def count_gc_at(dna):
renvoie une fois pour toutes la longueur de la chaine
totalcar=len(dna)
#initialisation des variables
nb_GC=nb_AT=0
#lancement de la boucle
for i in range (0,len(dna)):
on incrémente que les génes pairés
if dna[i]=="A"or dna[i]=="T": nb_AT+=1
elif dna[i]=="G"or dna[i]=="C": nb_GC+=1
#return "votre code"
# votre code
#calcul du ratio
GC=nb_GC/totalcar
AT=nb_AT/totalcar
#on renvoie les valeurs à la suite du script
return(GC,AT)
Je voudrais juste me permettre de préciser que demander les GC et AT n’est pas identique dans la formulation (selon moi) aux ratios de G+C et A+T dans la chaîne, je veux dire qu’on peut aussi comprendre qu’on cherche le nombre de fois où les bi-gram “AT” et "CG " sont trouvés dans la séquence.