首页
问答
文章
专栏
大咖
更多
话题
帮助
请输入关键字进行搜索
查看更多 "
" 的搜索结果
登录
blastp后,我想提取score的序列,请问使用python应该怎么写这个脚本
关注问题
回答问题
该问题已被锁定!
2
关注
785
浏览
blastp后,我想提取score的序列,请问使用python应该怎么写这个脚本
blast
系统是centos,python版本是2.7,第一列是我进行比对的主体,第二列是比对的另外四个物种,最后一列是蛋白比对的得分,现在我想要将spaf每个蛋白对应四个物种每个物种的蛋白最高得分那一行分别提取出来,请问使用python应该怎么编写这个脚本
阅读全文
收起全文
关注问题
回答问题
邀请回答
好问题
0
评论
收藏
举报
分享
复制链接
新浪微博
腾讯空间
微信扫一扫
展开
收起
0
评论
为什么被折叠?
0 个回复被折叠
1
回答
热门排序
最新排序
热门排序
只看楼主
孟浩巍
超级管理员
用户来自于: 北京市
2018-09-21 11:31
问题描述很清楚呀!我给你提供一个思路吧: 1. 建立dict, key是你的蛋白名,value 是一个list,长度为4 对应你的4个物种,里面需要存储你的分数值和蛋白名; 2. 开始循环,获得蛋白名,然后获得这个蛋白名对应的list,根据物种来比较分数,如果分数值大于list中存的分数就替换; 3. 结束循环,打印dict
阅读全文
收起全文
赞同
3
0
评论
分享
复制链接
新浪微博
腾讯空间
微信扫一扫
0
评论
关于作者
生信撰
注册会员
这家伙很懒,还没有设置简介
1
回答
0
文章
1
问题
问题动态
发布时间
2018-09-21 11:17
更新时间
2018-09-21 11:31
关注人数
2 人关注
相关问题
如何提取可变剪切位点?
958 浏览
2 关注
2 回答
0 评论
从bam文件还能提取出来原来比对前的fastaq文件吗?
1314 浏览
2 关注
2 回答
0 评论
sambamba 提取 LITCHI017845的比对情况
742 浏览
2 关注
1 回答
0 评论
我想从蠕虫库寄生虫 (wormbase.org)上下载全部的fa文件,请问怎么样才能批量的下载,同时下载的速度可以稍微快一点
1122 浏览
2 关注
3 回答
0 评论
SRR1924430的sra数据与hdf5提取的fastq不一致
573 浏览
1 关注
0 回答
0 评论
用atac-seq数据计算的TSS enrichment score
1316 浏览
2 关注
2 回答
0 评论
老铁们,我想问下,如果homer找到了motif1,那我有可能反过来知道哪些peak含有motif1么。
647 浏览
0 关注
1 回答
0 评论
conda怎么更新blastp到最新版本
1491 浏览
2 关注
2 回答
0 评论
NIPT算出了z-score后想要基于z>3的值对Z值进行校正该用什么方法
1221 浏览
2 关注
1 回答
0 评论
根据GFF和fasta等文件提取某一基因的ATG位点信息,比如具体的位置?
642 浏览
2 关注
1 回答
0 评论
推荐内容
问
已经有不同材料的全基因组二代测序数据,查找是否存在某个已知基因?
937 浏览
2 关注
2 回答
0 评论
问
本地nt库blastn时添加staxid,taxid号一列显示N/A
787 浏览
2 关注
1 回答
0 评论
问
求教一个blastn的问题
856 浏览
2 关注
1 回答
0 评论
All Rights Reserved Powered BY
WeCenter V4.1.0
© 2024
关于我们
社区规范
你的浏览器版本过低,可能导致网站部分内容不能正常使用!
为了能正常使用网站功能,请使用以下浏览器
Chrome
Firefox
Safari
IE 10+