Question

Perl写生信脚本，从一坨比对过后的DNA序列中摘出符合度最高的几条

Write a perl script that will parse the BLAST results to return the top hit for each sunflower unigene sequence in each BLAST search. Use an e-value of 1e-6 as your cut-off for calling a sequence a quality hit (ie., exclude sequences with e > 1e-6).

BLAST search（下图）我已经做出来了，但是后面这个脚本我没想法了。

第三列就是符合度。

Answer 1

open(A,"test.txt");
while(my $line=<A>){
    @sss=split("\t",$line); 
    if(! exists $hash{$sss[0]}){
        print $line;
        $hash{$sss[0]}=1;
    }

}