1. #!/usr/bin/perl -w
2. use strict;
3.
4. use ANNOserver;
5. use Data::Dumper;
6.
7. # Reads FASTA protein sequences from STDIN.
8. # Produces sequence ID, role, genome set name, confidence
9. # If "-blast" is specified, uses BLAST option.
10.
11. my $annoObject = ANNOserver->new();
12. my $blastOption = ($ARGV[0] && $ARGV[0] =~ /^-blast/i ? 1 : 0);
13.
14. my $results = $annoObject->assign_function_to_prot(-input => \*STDIN,
15. -kmer => 8,
16. -scoreThreshold => 3,
17. -assignToAll => $blastOption);
18. while (my $data = $results->get_next()) {
19. my ($id, $role, $genomeSet, undef, $hits) = @$data;
20. # Only proceed if a role was found.
21. if ($role) {
22. $genomeSet = '' if ! defined $genomeSet;
23. print join("\t", $id, $role, $genomeSet, $hits) . "\n";
24. }
25. }