Re: [程式] 如何利用bioperl去將blastp所得到的結果做分群
#! /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
討論串 (同標題文章)
BioMedInfo 近期熱門文章
PTT職涯區 即時熱門文章
