<![CDATA[
#!perl -w
# Add to the boulder sequences the reverse complement
use strict;
use SeqFun qw/revcom/;
use Boulder::Stream;
my $debug = grep(/-d/, @ARGV); # look for -d option
shift @ARGV if $debug; # get rid of -d option
my $stream = Boulder::Stream->new( );
my ($stone, $dna, $revseq);
while ($stone = $stream->get( )) {
$dna = $stone->Sequence;
if ( $debug ) {
print "The dna in BoulderProcess is ", length($dna),
" bp long.\n";
} else {
$revseq = revcom($dna);
$stone->insert(REVCOM => $revseq);
$stone->insert(DidRevcom => "OK");
$stream->put($stone)
}
}
]]>