Aide problème liste vide

Bonjour, j'ai effectué ce programme, mon problème c'est que dans ma fonction experiment, je n'arrive pas a stocké mais valeur dans ma liste data, cela viens surement de ma fonction hamming qui retourne une difference entre 2 sequence, ce dont j'ai besoin dans la liste (je veux une liste de toutes les differences entre les sequence generées pour ensuite en faire la moyenne avec ma fonction mean). Après beaucoup de manipulation, j'ai du mal a trouver une solution.
En esperant que vous puissiez m'aider

import random

def alphabet(name):
    global alpha 
    if name=="nucleic":
        alpha="acgt"
        print (alpha)
    if name=="protein":
        alpha="ARNDBCEQZGHILKMFPSTWYV"
        print (alpha)
    if name=="iupac_nucleic":
        alpha="ACGTRYSWKMBDHVN"
        print (alpha)
    if name=="iupac_protein":
        alpha="ARNDBCEQZGHILKMFPSTWYV"
        print (alpha)
        
def randseq(num,alpha): 
    seq1=[]
    for i in range(0,num):
        ahbon=random.randint(0,(len(alpha)-1))
        seq1.append(alpha[ahbon])
    seq="".join(seq1)
    return seq
    
def hamming(seq1,seq2):
    count=0
    j=0 
    for i in range(0,len(seq1)-1):
        if (seq1[j] != seq2[j]):
            count=count+1
        j=j+1    
    print (count)
    
def mutate(seq,num_subs):
    seq1=list(seq)
    seq2=list(seq)
    for i in range(num_subs):  
        ahboni=random.randint(0,(len(seq2)-1))
        seq2[ahboni]=randseq(1,alpha)
    return hamming(str(seq1),str(seq2))
    
def experiments(le,su,nb):
    for i in range(nb):
        seq=randseq(le,alpha)
        print (mutate(seq,su))
        data.append(mutate(seq,su))
    return data
    
import math 

def mean(data):
    somme=0
    for i in list(data):
        somme=somme+data[i]
    moyenne=somme/len(data)
    return moyenne

def variance(data):
    k = 0
    for x in data:
        k += 1
        if k == 1:
            M = x
            S = 0
        else:
            Mnext = M + (x - M) / k
            S = S + (x - M)*(x - Mnext)
            M = Mnext
    return (S/(k-1))

def std(data):    
    return math.sqrt(var)

data=[]
alpha="" 
name=str(input("Quel est le type?\n"))
alphabet(name)
le=int(input("Quelle longueur ?\n"))

su=int(input("Combien de substitution dans cette sequence?\n"))
nb=int(input("Combien de fois vous lez vous renouvelez l'experience?\n"))
experiments(le,su,nb)
print ("la liste de distance de hammin est",data)

#STAT
print ("la moyenne est",mean(data))
print ("la variance est",variance(data))
var=variance(data)
print ("le std est",std(data))

Bien a vous !