问题是这样的：有很多很多序列，几百条，想大致了解一下这些序列分别是什么样的微生物，如果一条一条去blast，那是相当的累。想找一个工具告诉我每条序列blast结果的前几条的名称是什么即可，不需要其它信息。

在网上找了一下，没找到合适的软件或工具，虽然有些关于批量blast的教程之类的，比如这个，但是给出的结果及其繁琐，很多不需要的信息。

后来发现Biopython可以很简单就进行批量Blast。只需先安装Python和Biopython，Python和Biopython的下载地址分别为：
http://www.python.org/download/
http://www.biopython.org/wiki/Download

继续阅读 »

前段时间做了一个克隆文库，有点大，拿到上百条序列，有了序列后，想选择一种合适的酶进行酶切，进而将克隆的结果与T-RFLP的结果联系起来。试了几个软件和一些在线酶切分析工具，好像都是可以用多种酶去切一个序列，而不能同时切多个序列。要是序列不多，一个个切切还可以，可是这上百条序列就需要几百次甚至更多的重复操作，不但能把手累的抽筋，眼都能被累的抽筋，实在让人无法忍受。于是决定自己做个小工具。

首先用了几个小时速成了一下Perl语言及Bioperl，发现这玩意处理字符串真的是太方便了！只需要下面这几行代码就可以自动或半自动同时查找多个序列酶切位点了(由于T-RFLP结果只与第一个酶切位点有关，因此只找每个序列的第一个切点)。
继续阅读 »