本地Blast的输入和输出问题
在这里学习了本地blast的使用，但是有个问题，本地blast默认的是标准输入和标准输出，由于需要，我希望将序列通过变量输入并且将比对结果输出到另一个变量。 比如 $reference中是一条序列（不是从文件读入） 而将结果输出到变量$blast_result（不是从文件读入）
但是用如下脚本是无法实现的 my $localblast=system("$blastall -p blastp -d $database -i $reference -o $blastresult"); 不知该怎么搞定啊。请大家帮忙，非常感谢！
• 2,942 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
Mapkkk •
•
written
by
Mapkkk •
1 answer
Log in to answer this question.
@mapkkk - I think you should edit your question since most of us can't really read Chinese characters :-)
I think you should use the perl module to do this instead of using 'system'. http://doc.bioperl.org/releases/bioperl-1.2/Bio/Tools/Run/StandAloneBlast.html