虽然有多种多样的办法来解析BLAST,即时单单对于Perl/BioPerl来说,你可以找到好多已有的脚本,如:blast_parsing.pl。但要想灵活的控制筛选条件及输出结果,还需要自己撰写脚本或者修改别人的脚本。所有这些的前提是你对BLAST结果中的各个项目有所了解,同时要知道如何提取出自己想要的项目。
此文的目的就是授渔而非授鱼。下面是几个使用BioPerl来解析BLAST的资料,只要你仔细阅读(主要是前两个)并进行简单的测试,就可以把BLAST玩弄于股掌之间了。
Parsing BLAST HSPs:简单快速的入门知道与脚本
NCBI-BLAST parsing problems:详细的解释与清晰的Method列表
Parsing BLAST and FASTA reports with Search and SearchIO:更多方法
”Parsing BLAST and FASTA reports with Search and SearchIO:更多方法“中的链接如果打不开的话,可以尝试一下下面几个:
http://classes.soe.ucsc.edu/bme060/Winter07/bptutorial.html#iii.4.2_parsing_blast_and_fasta_reports_with_search_and_searchio
http://www.cs.huji.ac.il/course/2003/bioskill/BioPerl/bptutorial.html#III_4_3_Parsing_BLAST_and_FASTA_
http://classes.soe.ucsc.edu/bme060/Winter07/bptutorial.html#iii.4.2_parsing_blast_and_fasta_reports_with_search_and_searchio