fasta文件排序

yang7473453

有fasta格式文件如下
>contig5
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
>contig3
GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
>contig1
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT

如何用perl排序成這個樣子
>contig1
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT
>contig3
GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
>contig5
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA

yang7473453

求指點回復 1# yang7473453


   

HardcoreAng

my %sequence;
my %order;
$/= ">";
open IN,"$ARGV[0]" or die;# 輸入文件
while (<IN>){
chomp;
my @array = split /\n/,$_;
my $query = shift @array;
    if (defined $query){
    my @query = split /contig/,$query;
    $order{$query[-1]} = $query;
    $sequence{$query} = join ("",@array);
    }
}
close (IN);

foreach my $m (sort {$a <=> $b} keys %order){
print ">$order{$m}\n$sequence{$order{$m}}\n";
}


測試沒問題

yang7473453

很感謝，很好用回復 3# HardcoreAng