今天的目标是随机产生碱基序列，并判断所产生的序列是否符合DNA碱基序列特征。
import random
Nucleotides=['A', 'T', 'C', 'G'] #创建包含DNA四个碱基的列表
DNA_seq=''.join([random.choice(Nucleotides) for i in range(500) ]) #随机产生DNA碱基序列
#函数判断产生的序列是否符合DNA序列
def validateSeq(dna_seq):
tmpseq=dna_seq.upper() #upper函数，将输入的字母全部换成大写字母
for nuc in tmpseq: #遍历产生的序列
if nuc not in Nucleotides: #if函数判断
return False #不是返回错误
return tmpseq #是返回碱基序列
DNAStr=validateSeq(DNA_seq)
Print(DNAStr)
最后
以上就是娇气哈密瓜最近收集整理的关于python 处理生物信息（1）的全部内容，更多相关python内容请搜索靠谱客的其他文章。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复