die "perl $0 <fq1><OUT1>" unless(@ARGV==2);
use Bio::SeqIO;
use Bio::Seq;
open my $FQ1 ,"zcat $ARGV[0] |" or die "$!";
my$fq1=Bio::SeqIO->new(-fh=>$FQ1,-format=>'fastq');
open my $GZ1 ,"| gzip > $ARGV[1]" or die $!;
my$fqo1=Bio::SeqIO->new(-fh=>$GZ1,-format=>'fasta');
my$i=0;
while ( my $obj1=$fq1->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);
#}
}
$fq1->close();
$GZ1->close();
