RBP AS生信流程(N)RBP表达量提取

use strict;
use warnings;
my %hash=();
open(RF,"RBPgene.txt") or die $!;
while(my $line=<RF>){
    chomp($line);
    $hash{$line}=1;
}
close(RF);
open(RF,"gene expression.txt") or die $!;
open(WF,">RBP expression.txt") or die $!;
while(my $line=<RF>){
    if($.==1){
        print WF $line;
        next;
    }
    my @arr=split(/\t/,$line);
    my @zeroArr=split(/\|/,$arr[0]);
    if(exists $hash{$zeroArr[0]}){
        print WF $line;
        delete($hash{$zeroArr[0]});
    }
}
close(WF);
close(RF);
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容