BioPython ① | 统计蛋白序列中20种氨基酸的的个数和频率

统计蛋白序列中20种氨基酸的的个数和频率

题目

MNAPERQPQPDGGDAPGHEPGGSPQDELDFSILFDYEYLNPNEEEPNAHKVASPPSGPAYPDDVLDYGLKPYSPLASLSGEPPGRFGEPDRVGPQKFLSAAKPAGASGLSPRIEITPSHELIQAVGPLRMRDAGLLVEQPPLAGVAASPRFTLPVPGFEGYREPLCLSPASSGSSASFISDTFSPYTSPCVSPNNGGPDDLCPQFQNIPAHYSPRTSPIMSPRTSLAEDSCLGRHSPVPRPASRSSSPGAKRRHSCAEALVALPPGASPQRSRSPSPQPSSHVAPQDHGSPAGYPPVAGSAVIMDALNSLATDSPCGIPPKMWKTSP

提示和要求：用循环、字典等实现，并将频率输出到名为D:\test\frq.txt的文本文件中，且文件的第一行必须是标题行，包含20列，每一列为对应氨基酸的单字母符号，列间用Tab键’\t’分割。

解答过程

这题用过一个字符串的count()函数会比较简单实现统计频数的操作。

# -*- coding: utf-8 -*-
# @Author  : 张英豪
# @Time    : 2022-11-30
# @Function: 统计蛋白质序列中氨基酸的频率# 蛋白质序列
Protein_Sequence = 'MNAPERQPQPDGGDAPGHEPGGSPQDELDFSILFDYEYLNPNEEEPNAHKVASPPSGPAYPDDVLDYGLKPYSPLASLSGEPPGRFGEPDRVGPQKFLSAAKPAGASGLSPRIEITPSHELIQAVGPLRMRDAGLLVEQPPLAGVAASPRFTLPVPGFEGYREPLCLSPASSGSSASFISDTFSPYTSPCVSPNNGGPDDLCPQFQNIPAHYSPRTSPIMSPRTSLAEDSCLGRHSPVPRPASRSSSPGAKRRHSCAEALVALPPGASPQRSRSPSPQPSSHVAPQDHGSPAGYPPVAGSAVIMDALNSLATDSPCGIPPKMWKTSP'
# 氨基酸的种类
Amino_Acid_type = 'ACDEFGHIKLMNPQRSTVWY'
# 统计氨基酸的频率，保留三位小数
Amino_Acid_freq = {}
for i in range(len(Amino_Acid_type)):Amino_Acid_freq[Amino_Acid_type[i]] = round(Protein_Sequence.count(Amino_Acid_type[i]) / len(Protein_Sequence),3)# 创建并打开一个文本文件
file = open('D:/test/frq.txt', 'w')
for k in Amino_Acid_freq.keys():file.write(str(k) + '\t')
file.write('\n')
for v in Amino_Acid_freq.values():file.write(str(v) + '\t')
# 关闭文件
file.close()

BioPython ① | 统计蛋白序列中20种氨基酸的的个数和频率相关推荐

excel如何快速统计一列中相同数值出现的个数
excel中如何快速统计一列中相同数值出现的个数,数据如下图,A1:A100为账号,有重复,需统计每个人出现了多少次. 1.选中A1:A100的数据范围,选插入-数据透视表,按照默认设置,创建数据透视 ...
【python作业】编写一个函数，由实参传来一个字符串，统计此字符串中字母和数字的个数，在主函数中输入字符串并输出上述的结果。
编写一个函数,由实参传来一个字符串,统计此字符串中字母和数字的个数,在主函数中输入字符串并输出上述的结果. def fun(s):count1 = 0count2 = 0for i in s:if i ...
7-64 统计大写辅音字母 (15 分)英文辅音字母是除A、E、I、O、U以外的字母。本题要求编写程序，统计给定字符串中大写辅音字母的个数。PTA：中M2021春C、Java入门练习第I段
7-64 统计大写辅音字母 (15 分) 英文辅音字母是除A.E.I.O.U以外的字母.本题要求编写程序,统计给定字符串中大写辅音字母的个数. 输入格式: 输入在一行中给出一个不超过80个字符.并以回 ...
英文辅音字母是除A、E、I、O、U以外的字母。本题要求编写程序，统计给定字符串中大写辅音字母的个数。
英文辅音字母是除A.E.I.O.U以外的字母.本题要求编写程序,统计给定字符串中大写辅音字母的个数. 输入格式: 输入在一行中给出一个不超过80个字符.并以回车结束的字符串. 输出格式: 输出在一行中 ...
利用Transformer替代MSA从蛋白序列中学习Contact Map
点击上方"CVer",选择加"星标"置顶重磅干货,第一时间送达本文转载自:GoDesign --背景-- 基于深度学习的蛋白结构预测在近年来取得了不少突破, ...
统计任意字符串中回文字符串的个数
题目来源: 庞果英雄会题目详情: 回文字符串是指从左到右和从右到左相同的字符串,现给定一个仅由小写字母组成的字符串,你可以把它的字母重新排列,以形成不同的回文字符串. 输入:非空仅由小写字母组成的字 ...
统计一个数组中正数和负数的个数
1.遍历数组,将正数重新放在一个数组中,负数放在另一个数组中 1 #include <iostream> 2 using namespace std; 3 int main() 4 { 5 ...
统计csv表格中某一元素的个数（count函数）
我有一个csv表格(data.csv),一共1915列,732行(绝大多数的元素是0,偶尔会有1). 现在我想统计 1 在每一列出现的次数.在整个表中出现的次数. data[data == 1].co ...
EXCEL 快速统计一列中相同数值出现的个数
之前用EXCEL需要统计频次用到的.这篇文章放在笔记本很久了,记录一下. 首先,选中A1:A51的数据范围,选择[插入]-[数据透视表]. 按照默认设置,在新工作表中创建数据透视表,点击[确定]按钮. ...

BioPython ① | 统计蛋白序列中20种氨基酸的的个数和频率

统计蛋白序列中20种氨基酸的的个数和频率

题目

解答过程

BioPython ① | 统计蛋白序列中20种氨基酸的的个数和频率相关推荐

最新文章

热门文章