/usr/share/doc/libace-perl/examples/upstream2.pl is in libace-perl 1.92-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | #!/usr/bin/perl
# This example will pull out N base pairs upstream of each predicted
# gene in C. elegans
use lib '..','../blib/arch';
use Ace::Sequence;
use strict vars;
use constant HOST => $ENV{ACEDB_HOST} || 'formaggio.cshl.org';
use constant PORT => $ENV{ACEDB_PORT} || 200005;
$|=1;
my $upstream = shift || die "Usage: upstream.pl <size (bp)>\n";
my $db = Ace->connect(-host=>HOST,-port=>PORT) || die "Connection failure: ",Ace->error;
warn "Fetching all predicted genes, please wait....\n";
#my @genes = $db->fetch('Predicted_Gene' => '*');
my @genes = $db->fetch(Predicted_Gene => '4R79.2');
for my $gene(@genes) {
next unless my $seq = Ace::Sequence->new(-seq=>$gene,
-offset => (-$upstream-1),
-length => ($upstream)
);
print $gene,"\t",$seq->dna,"\n";
}
|