使用python简单实现K核苷酸频率（KNF，k-nucleotide frequencies）或K-mer频率

K核苷酸频率（KNF，k-nucleotide frequencies）或K-mer频率
KNF描述了序列中存在k个核苷酸的所有可能的多核苷酸的频率。如果k=2，则计算的为双核苷酸频率（即AA、AT、AG、AC、……TT），共42=16种；如果k=3，则计算的为双核苷酸频率（即AAA、AAT、AAG、AAC、……TTT），共43=64种；以此类推。
K-mer频率方法与上面所述相同。

from itertools import product  #定义包
def nucleotide_frequency(seq,k):  #seq为输入的多个序列，k为几核苷酸
    X = []
    char = []  
    for i in product('ACGT', repeat = k):  #笛卡尔积（有放回抽样排列）
        char.append(''.join(i))  #去除引号后存放到char里，即把('A,A,A')转变成（AAA）形式
    for s in seq:  #循环调用每一个序列
        frequence = []  #定义空列表，用于存放计算好的频率
        for a in char:  #循环取char里的k核苷酸类型
            number = 0
            for j in range(len(s)-k+1):  #在单个序列中进行循环取k间隔片段
                if s[j:j+k] == a:  #判断k间隔片段是否和char里的k核苷酸类型相等
                    number += 1  
            char_frequence = number/(len(s)-k+1)  #计算频率，此时序列长度=原序列长度-k+1
            frequence.append(char_frequence)
        X.append(frequence)
    return X

#用具体数据进行调用
knf = nucleotide_frequency(data,k)  #data为具体数据，k的值根据需要自己设定