该问题已被锁定!
2
关注
698
浏览

blastp后,我想提取score的序列,请问使用python应该怎么写这个脚本

为什么被折叠? 0 个回复被折叠
孟浩巍 超级管理员 用户来自于: 北京市
2018-09-21 11:31
问题描述很清楚呀!我给你提供一个思路吧:   1. 建立dict, key是你的蛋白名,value 是一个list,长度为4 对应你的4个物种,里面需要存储你的分数值和蛋白名;   2. 开始循环,获得蛋白名,然后获得这个蛋白名对应的list,根据物种来比较分数,如果分数值大于list中存的分数就替换;   3. 结束循环,打印dict

关于作者

问题动态

发布时间
2018-09-21 11:17
更新时间
2018-09-21 11:31
关注人数
2 人关注

相关问题

如何提取可变剪切位点?
从多个.fa文件中提取以“poptri.”开头的蛋白序列
想要提取指定基因的行信息
怎么准确提取某一基因的ATG起始密码子位点
从bam文件还能提取出来原来比对前的fastaq文件吗?
conda怎么更新blastp到最新版本
xpclr的XPCLR_score的inf值应该怎么处理
sambamba 提取 LITCHI017845的比对情况
用atac-seq数据计算的TSS enrichment score
老铁们,我想问下,如果homer找到了motif1,那我有可能反过来知道哪些peak含有motif1么。

推荐内容

本地nt库blastn时添加staxid,taxid号一列显示N/A
求教一个blastn的问题
已经有不同材料的全基因组二代测序数据,查找是否存在某个已知基因?
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024