提取Blast2go blast结果中的一部分

问题

有两个样本,两者之间有些序列是重复的,所以其中一个过滤掉了重复的序列后才进行 Blast2go 分析,,另外一个则全体。 如今须要把过滤掉那部分序列的也给合并过来。服务器

思路

能够根据序列 id,直接 grep xxxx *.xml,从另外一样本那里得到 XML 文件,而后拷贝到新的目录下,再导入到相应 b2g 文件中。code

# 例如
grep -P '\bTSA_246786\b' *.xml

注意:xml

  • 在脚本里要转义:\\b
  • Blast2go 本地 BLAST 的结果格式为 XML2

虽然有点慢,但的确可行,毕竟 BLAST 只会更慢,咱没有服务器,只有破电脑一台。ast

相关文章
相关标签/搜索