同意...楼主老让别人给楼主打工,别人也要工作啊...作者: leigh111 时间: 2011-05-09 09:34
open (FH,'gene.txt') || die "$!";
open (IN,'go.txt') || die "$!";
my %seq;
while (<FH>) {
chomp;
if (/^BGIBMGA\d+\-PA/) {
$seq{$&} = $_;
}
}
close FH;
while (<IN>) {
chomp;
if (exists $seq{$_}) {
print "$seq{$_}\n";
}
else {
print "$_\n";
}
}
close IN;作者: kingwmj 时间: 2011-05-09 10:05
#!/usr/bin/perl
use strict;
use warnings;
use 5.010;
say "Enter you GO's name here: ";
chomp(my $go_name = <STDIN>
my $file = "/root/perl/GO";
my @temp;
my %ghash;
open(GO,"<","$file" || die "Can't open target files!\n";
while (<GO> {
@temp = split(" ",$_);
print if /\b$go_name\b/;