MJさむらいの実験室

ブラストクラストを実行するとクラスタリングされた遺伝子の名前だけが吐き出される。

各行にはクラスタリングされた遺伝子の名前がリストアップされる。

これらのうち一番左側に書かれている遺伝子がクラスタリングされた最も長い配列。

クラスタリングされた遺伝子のうち最も長い遺伝子名のみを抽出するプログラムを書いてみた。

以下プログラム

#!/usr/bin/perl -w
use strict;
use warnings;
use POSIX;

my $counter = 0;

my $name = <STDIN>;

open (FASTANAME, $name);

while (<FASTANAME>)
{

 my @filename = split( " ", $_);
 my @count = @filename;
if (@count <= 1)
{
 exit;
}
print ">$filename[0]\n";

}
close (FASTANAME);