Perl記念日

BLASTの結果から各問い合わせ配列のtop hitだけをとりたい、ってとき。

#!/usr/bin/perl

my $prev_query = '';

while(<>) {
         my($query) = split(/\t/);
	 if($query ne $prev_query) {
	           print $_;
	}
	$prev_query = $query;
}
	   

昨日、BLASTのオプションを探し回って結局ダメだった上記のやりたいことが、これで一発。
これは本当に「呪文」だ。。