1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103
| #depart code
nbA = nbC = nbG = nbT = TotalNb = NbStepsRight = NbStepsUp = 0
#SeqLength = L = I = InitW = 0
#XEndSegment = YEndSegment = Step = 0.00
from turtle import *
def compte_bases(chaine):
nbA = nbC = nbG = nbT = TotalNb = NbStepsRight = NbStepsUp = 0
for base in chaine :
#elif base i "A":
if base in "A":
nbA +=1
elif base == "C":
nbC +=1
elif base == "G":
nbG +=1
elif base == "T":
nbT +=1
TotalNb +=1
NbStepsRight = nbC - nbG
NbStepsUp = nbA - nbT
return(TotalNb, nbA, nbC, nbG, nbT, NbStepsRight, NbStepsUp )
def dessine_un_segment_adn2(x,y):
showturtle()
pendown()
goto(x,y)
#setposition(x,y)
#hideturtle()
#x1=x
#y1=y
#on entre la chaine
chaine_nucleotides = input("écris ta chaine de nucléotides : ")
#on écrit le code
Lg=longueur_chaine_nucleotides = len(chaine_nucleotides)
print("Lg: " + str(Lg))
n=0
pas=8
r=recouvrement=2
nbNC = Nombre_de_mini_chaines_nucléotides = Lg/(pas-r)
print("nbNC: " + str(nbNC))
# On ne veut pas des derniers nucléotides inférieur à pas
nbNC = int(nbNC)
print("nbNC: " + str(nbNC))
for mini_chaines_nucleotides in range (n,nbNC):
a=n*(pas-r)
b=(n*(pas-r)+pas)
mini_chaines_nucleotides = chaine_nucleotides[a:b]
print(mini_chaines_nucleotides)
#on compte
glu = compte_bases(mini_chaines_nucleotides)
print(glu)
#return renvoie un tuple
TotalNb = glu[0]
nbA = glu[1]
nbC = glu[2]
nbG = glu[3]
nbT = glu[4]
NbStepsRight = glu[5]
NbStepsUp = glu[6]
print("total Nb : " + str(TotalNb) )
print("NBA: " + str(nbA))
print("NBC: " + str(nbC))
print("NBG: " + str(nbG))
print("NBT: " + str(nbT))
print("NBStepsRight: " + str(NbStepsRight))
print("NBStepsUP: " + str(NbStepsUp))
#on dessine
dessine_un_segment_adn2(NbStepsRight,NbStepsUp)
n +=1
print(n) |
Partager