poeli/extract_reads_from_fastq_by_desc.pl

## extract_reads_from_fastq_by_desc.pl
#!/usr/bin/env perl

$string=$ARGV[0];
$cnt=0;
$found=1;

while(<STDIN>){
	if( $cnt%4==0){
		if( /$string/i ){
			$found=1;
			print $_;
		}
		else{
			$found=0;
		}
	}
	elsif( $found ){
		print $_;
	}
	$cnt++;
}
	#!/usr/bin/env perl

	$string=$ARGV[0];
	$cnt=0;
	$found=1;

	while(<STDIN>){
	if( $cnt%4==0){
	if( /$string/i ){
	$found=1;
	print $_;
	}
	else{
	$found=0;
	}
	}
	elsif( $found ){
	print $_;
	}
	$cnt++;
	}