[Bioperl-l] load_seqdatabase.pl problem
Hilmar Lapp
hlapp at gnf.org
Mon Mar 1 23:16:31 EST 2004
I suspect the database wasn't fresh at the time you loaded the NCBI
taxonomy. Specifically, the offending species probably was already in
there through some sequence load.
Alternatively, the EMBL SeqIO parser fails to parse out the NCBI
taxonID and the species name slightly mismatches with the one in the
NCBI taxonomy database.
What is the result if you issue the following query:
select t.* from taxon t, taxon_name tn
where tn.taxon_id = t.taxon_id
and tn.name = 'Aspergillus oryzae'
and tn.name_class = 'scientific name';
-hilmar
On Monday, March 1, 2004, at 04:57 AM, Eric Bazin wrote:
> Hi.
>
> I try to import embl sequence datas in a biosql schema using the perl
> script "load_seqdatabase.pl". Previously I loaded without any problem
> the taxon schema from the NCBI with load_ncbi_taxonomy.pl. Can anybody
> give
> me a help for this?
>
> This is what I received:
>
> bash#./load_seqdatabase.pl --dbname polymorphix --format embl
> /home/eric/polybase.dat
>
> Could not store AB022832:
> ------------- EXCEPTION -------------
> MSG: Unique key query in Bio::DB::BioSQL::SpeciesAdaptor returned 3
> rows
> instead of 1. Query was [name_class="scie
> ntific name",binomial="Aspergillus oryzae"]
> STACK Bio::DB::BioSQL::BasePersistenceAdaptor::_find_by_unique_key
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/BioSQL/Ba
> sePersistenceAdaptor.pm:960
> STACK Bio::DB::BioSQL::BasePersistenceAdaptor::find_by_unique_key
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/BioSQL/Bas
> ePersistenceAdaptor.pm:856
> STACK Bio::DB::BioSQL::BasePersistenceAdaptor::create
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/BioSQL/BasePersistence
> Adaptor.pm:183
> STACK Bio::DB::Persistent::PersistentObject::create
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/Persistent/PersistentObj
> ect.pm:243
> STACK Bio::DB::BioSQL::BasePersistenceAdaptor::create
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/BioSQL/BasePersistence
> Adaptor.pm:170
> [root at thaler biosql]# :BasePersistenceAdaptor::store
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/BioSQL/BasePersistenceA
> daptor.pm:253
> STACK Bio::DB::Persistent::PersistentObject::store
> /usr/lib/perl5/site_perl/5.8.0/Bio/DB/Persistent/PersistentObje
> ct.pm:270
> STACK (eval) ./load_seqdatabase.pl:517
> STACK toplevel ./load_seqdatabase.pl:500
>
> --------------------------------------
>
> Thanks a lot
>
> Eric bazin
>
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
--
-------------------------------------------------------------
Hilmar Lapp email: lapp at gnf.org
GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
-------------------------------------------------------------
More information about the Bioperl-l
mailing list