#!/usr/bin/perl
#---blast output file
my $infile = "blastout";
open(IN, "$infile");
while((my $line =
chomp($line);
if ($line =~ /\
push(@hitnums,split(/\<\/Hit_num\>/,(split(/\
}
if ($line =~ /\
push(@hitdefs,split(/\<\/Hit_def\>/,(split(/\
}
if ($line =~ /\
push(@hitlens,split(/\<\/Hit_len\>/,(split(/\
}
if ($line =~ /\
push(@bitscores,split(/\<\/Hsp_bit-score\>/,(split(/\
}
if ($line =~ /\
push(@scores,split(/\<\/Hsp_score\>/,(split(/\
}
if ($line =~ /\
push(@evalues,split(/\<\/Hsp_evalue\>/,(split(/\
}
if ($line =~ /\
push(@qfroms,split(/\<\/Hsp_query-from\>/,(split(/\
}
if ($line =~ /\
push(@qtos,split(/\<\/Hsp_query-to\>/,(split(/\
}
if ($line =~ /\
push(@hfroms,split(/\<\/Hsp_hit-from\>/,(split(/\
}
if ($line =~ /\
push(@htos,split(/\<\/Hsp_hit-to\>/,(split(/\
}
if ($line =~ /\
push(@qframes,split(/\<\/Hsp_query-frame\>/,(split(/\
}
if ($line =~ /\
push(@hframes,split(/\<\/Hsp_hit-frame\>/,(split(/\
}
if ($line =~ /\
push(@identities,split(/\<\/Hsp_identity\>/,(split(/\
}
if ($line =~ /\
push(@positives,split(/\<\/Hsp_positive\>/,(split(/\
}
if ($line =~ /\
push(@algnlens,split(/\<\/Hsp_align-len\>/,(split(/\
}
}
No comments:
Post a Comment