作者hgsfhevil (evil)
看板BioMedInfo
标题Re: [程式] 如何利用bioperl去将blastp所得到的结果做分群
时间Sun Aug 17 22:31:04 2008
#! /usr/bin/perl -w
use strict;
use Bio::SearchIO;
my $blast = new Bio::SearchIO(-format => 'blast',-file => $ARGV[0]);
my %Name;
#my $result=$blast->next_result;
while(my $result = $blast->next_result )
{
while(my $sbjct = $result->next_hit)
{
while(my $hsp = $sbjct->next_hsp)
{
$Name{$sbjct->name} = 1 if $hsp->frac_identical >= 0.0;
}
}
print join("\t", sort keys %Name) , "\n";
%Name = ();
}
这是修改好的
他会根据序列相似度去parse
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.138.155.221