#!/usr/bin/perl
use strict;
use warnings;
my $annotation = shift; #annotation文件
my $target_gene = shift; #目的基因
my $name = shift; #命名
open IN,$annotation;
open INID,$target_gene;
open(OUT, ">/data/00/user/user159/wm/OGgene/dataAnalysis/Go_anno/$name.GO.gene.") || die "$!";
my %keep_id=();
while(my$line=<INID>){
chomp $line;
$keep_id{$line}=1;
}
close(INID);
while(my$line=<IN>){
chomp $line;
my @tmp=split(/\s/, $line);
if (exists ($keep_id{$tmp[0]})){
print OUT$line."\n";
}
}
close(IN);
close(OUT);