perl中如何輸出指定ID對(duì)應(yīng)的雙端fastq文件

這篇文章主要介紹“perl中如何輸出指定ID對(duì)應(yīng)的雙端fastq文件”的相關(guān)知識(shí)，小編通過實(shí)際案例向大家展示操作過程，操作方法簡單快捷，實(shí)用性強(qiáng)，希望這篇“perl中如何輸出指定ID對(duì)應(yīng)的雙端fastq文件”文章能幫助大家解決問題。

給一個(gè)fastqID列表，輸出指定ID對(duì)應(yīng)的雙端fastq文件：

perl get_fq_by_id.pl <id><fq1><fq2><OUT1><OUT2>

給一個(gè)fastqID列表，輸出指定ID對(duì)應(yīng)的雙端fastq文件：

die "perl $0 <id><fq1><fq2><OUT1><OUT2>" unless(@ARGV==5);
use Bio::SeqIO;
use Bio::Seq;

open my $FQ1 ,"zcat $ARGV[1] |" or die "$!";
open my $FQ2 ,"zcat $ARGV[2]|" or die "$!";
my$fq1=Bio::SeqIO->new(-fh=>$FQ1,-format=>'fastq');
my$fq2=Bio::SeqIO->new(-fh=>$FQ2,-format=>'fastq');

open my $GZ1 ,"| gzip >$ARGV[3]" or die $!;
open my $GZ2 ,"| gzip >$ARGV[4]" or die $!;
my$fqo1=Bio::SeqIO->new(-fh=>$GZ1,-format=>'fastq');
my$fqo2=Bio::SeqIO->new(-fh=>$GZ2,-format=>'fastq');
my%keep=(
);
open IN ,"$ARGV[0]" or die "$!";
while(<IN>){
chomp;
next if /^#/;
my@tmp=split(/\s+/);
$keep{$tmp[0]}=1;
}
close(IN);
my$i=0;
while ( my $obj1=$fq1->next_seq() and my $obj2=$fq2->next_seq() ) {
my$id1=$obj1->id;
my$id2=$obj2->id;
#print "$id1,$id2\n";
#die;
if( exists $keep{$id1} or exists $keep{$id2}){
$fqo1->write_seq($obj1);
$fqo2->write_seq($obj2);
}
}
$fq1->close();
$fq2->close();
$GZ1->close();
$GZ2->close();

關(guān)于“perl中如何輸出指定ID對(duì)應(yīng)的雙端fastq文件”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識(shí)，可以關(guān)注億速云行業(yè)資訊頻道，小編每天都會(huì)為大家更新不同的知識(shí)點(diǎn)。