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