取转录本fasta最长的当作基因fasta

#!/usr/bin/env perl
use warnings;
use strict;
use Bio::SeqIO;

die "perl $0 <fasta> > <outfile>\n" if(@ARGV != 1);

my @len = `fastalength $ARGV[0]`;
my %xsh = map { chomp; my($length,$name)=split /\s+/; $name=>$length } @len;

my %hash;
open FA, $ARGV[0] or die $!;
while(<FA>)
{
	chomp;
	next if($_ !~ /^>/);
	if(/gene:([^ ]+) transcript:([^ ]+)/)
	{
		if(exists $xsh{$2})
		{
			push @{$hash{$1}}, "$2 $xsh{$2}";
		}
	}
}

my %trans;
open OUT, ">longest_transcript.list" or die $!;
foreach my $g(keys %hash)
{
	my $tmp = 0;
	my $lt;
	foreach my $t(@{$hash{$g}})
	{
		my @x = split / /, $t;
		($x[1] > $tmp) ? $lt = $x[0] : next;
	}
	$trans{$lt} = 0;
	print OUT "$lt\n";
}

$/ = "\n>";
open FB, $ARGV[0] or die $!;
while(<FB>)
{
	chomp;
	s/^>//;
	my @tmp = split /\n/;
	my @t = s