用perl根据ID号在fasta文件中提取序列程序

#!/usr/bin/perl -w

use strict;

die "perl $0 \n" unless(@ARGV == 2);

my %hash;

open IN,$ARGV[0];

while(){

chomp;

$_=~s/[> \t]+//g;

$hash{$_} = 1;

}

close IN;

open IN,$ARGV[1];

$/=">";

;

$/= "\n";

while(){

chomp;

my $name = $_;

$/=">";

my $seq= ;

chomp $seq;

$seq =~s/\n//g;

$/ = "\n";

if(exists($hash{$name})){

print ">$name\t$seq\n";

}

}

close IN;

file:///D|/生物软件/程序文件/perl/用perl根据ID号在fasta文件中提取序列程序.txt[2016/1/7 20:07:29]