None available.
sub run
{
my $self = shift if (defined(caller(1)));
my $source_id = shift;
my $species_id = shift;
my $files = shift;
my $release_file = shift;
my $verbose = shift;
my $file = @{$files}[0];
my @xrefs;
my $ag_io = $self->get_filehandle($file);
if ( !defined $ag_io ) {
print STDERR "Could not open $file\n";
return 1;
}
my $probe;
while ( $_ = $ag_io->getline() ) {
chomp;
my $xref;
$_ =~ s/\015//g;
if(/^>(.+)/){
$probe = $1;
}
else{
my $sequence = $_;
$sequence =~ s/\n//g;
$xref->{ACCESSION} = $probe;
$xref->{LABEL} = $probe;
$xref->{SEQUENCE} = $sequence;
$xref->{SOURCE_ID} = $source_id;
$xref->{SPECIES_ID} = $species_id;
$xref->{SEQUENCE_TYPE} = 'dna';
$xref->{STATUS} = 'experimental';
push @xrefs, $xref;
}
}
$ag_io->close();
XrefParser::BaseParser->upload_xref_object_graphs(\@xrefs);
print scalar(@xrefs) . " Agilent xrefs succesfully parsed\n" if($verbose);
return 0;
}
1; } |