24小时热门版块排行榜    

查看: 921  |  回复: 0

badsunday

新虫 (小有名气)

[求助] windows下bioperl怎么调用clustalw比对序列。

最近LZ对bioperl比较着迷,windows环境下,安装好了biopel和clustalw。clustalw路径为C:\Program Files (x86)\ClustalW2.
代码如下
use Bio::Tools::Run::Alignment::Clustalw;
use Bio::AlignIO;
@params= ('ktuple' => 2, 'matrix' => 'BLOSUM');
$factory= Bio::Tools::Run::Alignment::Clustalw->new(@params);
$inputfilename = '1.fa';
$in = $factory->align($inputfilename);

$out = Bio::AlignIO->new(-file => ">outputfilename",
                         -format => 'pfam');
                        
while ( my $aln = $in->next_aln ) {
  $out->write_aln($aln);
}

运行后出现
windows下bioperl怎么调用clustalw比对序列。
我猜原因就出在bioperl不能找到Clustalw程序。由于在windows环境下,代码开始端添加BEGIN { $ENV{CLUSTALDIR} = 'C:\Program Files (x86)\ClustalW2' }也不行。windos环境下虽然system函数可以调用Clustalw,但终究是不方便。
同样的道理PAML中Codeml怎么调用,Bio::Tools:: Phylo:: PAML::Codeml,PAML中Codeml被我安在了'C:\paml4.7\bin\codeml.exe'
跪求高手赐教。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 badsunday 的主题更新
信息提示
请填处理意见