转自 https://www.plob.org/article/4603.html 具体请去上面的网页查看. my $DNA="ATGCCCGGT";my $pep=&TranslateDNASeq($DNA); sub TranslateDNASeq{ use Bio::Seq; (my $dna)=@_; my $seqobj=Bio::Seq->new(-seq =>$dna, -alphabet =>'dna'); return
编码 from __future__ import division def get_aa_percentage(protein, aa_list=['A','I','L','M','F','W','Y','V']): protein = protein.upper() protein_length = len(protein) total = 0 for aa in aa_list: aa = aa.upper() aa_count = protein.count(aa) total += a
转载于 Original 2017-06-20 liuhui 生信百科 KEGG 数据库中,把功能相似的蛋白质归为同一组,然后标上 KO 号.通过相似性比对,可以为未知功能的蛋白序列注释上 KO 号.通过KEGG数据库的注释极大的方便我们进行生物学通路的研究,可以直接查看物种某条生物学通路上基因的存在情况. 最简单的方法是看公司给的KEGG注释或者直接下载本物种每个基因的注释结果(比如,植物Phytozome:动植物Ensemble),然后对应到自己的差异基因集里面. 当然如果自己的物种没有KE