| 查看: 1426 | 回复: 0 | |||
[交流]
高通量数据做本地blast时，如何将97%以上的比对结果全部显示出来？
|
|
各位有懂的吗？之前在另一个区求助，知晓的烦请点一下链接 http://muchong.com/bbs/viewthread.php?tid=7834451 各位虫友，我有一批16S的高通量数据，然后我现在想做一个本地blast，把我手上的一部分序列（isolate.fasta）与高通量数据进行blast，然后把所有高于97%的结果全部显示出来，请大家帮忙。 我是这样做的，首先用格式化本地数据库（即16S高通量数据515F-907R.format.fq.fasta），用的是formatdb命令 formatdb -i 515F-907R.format.fq.fasta -p F -t My database -o T 然后生成三个扩展名分别是：.pin、.psq、.phr的文件，然后就开始进行本地blast，用blastall命令 blastall -p blastn -i isolate.fasta -d 515F-907R.format.fq.fasta -o myresult.blastout -T T -e 1e-10 -m 9 这样运行出来会有一个所有isolate.fasta的序列与库序列比对的结果，但是每个比对结果只显示前250个，有些可能比到第250个相似性也是99%以上，就说说还有许多库序列与query的相似性高于97%，但是这部分就不显示出来了。 请问一下大家，可以用什么程序解决这个问题，可不可以设置一个97%的阈值，比对高于这个值得全部显示出来，低于97%的不显示。刚刚接触生物信息方面的工作，很多不太清楚，谢谢大家帮忙！ |
回复此楼» 猜你喜欢
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有180人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
-
氢氧化钾活化操作
已经有18人回复
-
散金
已经有68人回复
-
生物质超级电容器
已经有0人回复
-
ALK5 高选择性小分子 LY-3200882：可放大合成路线与质控要点解析
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
请教，转录组测序与基因组测序结果分析
已经有4人回复
- 有木有一个基因组疾病的数据库？
已经有5人回复
- 测序完之后，怎么看测序结果？请指教
已经有6人回复
- RACE拼接后扩增全长的测序结果不一样
已经有5人回复
- 测转录组遇到的一个问题 求高手解答
已经有4人回复
- 基因组测序结果问题
已经有7人回复
- 各位，我800bp的基因测序后有3个点突变，后续能做表达用吗？
已经有7人回复
- 我扩了某个物种的某个基因，测序结果比对后，我惊呆了
已经有9人回复
- Illumina高通量RNA测序数据分析中不同标准化方法的综合比较
已经有6人回复
- 鉴定一个基因的特定功能需要从哪些方面开展研究才能确定？
已经有10人回复
- 基因测序对比之后，该怎样判断下一步怎么做。
已经有3人回复
- 植物转录组测序结果
已经有9人回复
- 急！急！急！对转录组测序结果进行验证，如何选定基因
已经有9人回复
- 转录组数据分析免费软件
已经有11人回复
- 测序结果的疑问
已经有7人回复
- 【求助/交流】玉米基因组库
已经有7人回复
- 【求助/交流】primer-blast
已经有3人回复
- 生物信息学专题－生物版，医学版和信息科学版共同创建
已经有87人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30