[Bioperl-l] retrieval of PRELIMINARY uniprot sequences using Bio::Registry fails
daniel.lang at biologie.uni-freiburg.de
Tue Sep 5 09:57:59 UTC 2006
sorry for the belated response!
I've compiled you a set of 100 PRELIMINARY entries from the latest
uniprot_trembl release. I've tried to reproduce the bug using only these
as input to build an index, but (sadly) all of them can be retrieved
using the latest checkout:-(
Maybe its not connected to these entries after all, but the size or some
other feature of the uniprot distribution?
I now could make it work using the 1.5.1 release.
Originally, I've built the index using flat protocol, when I try bdb and
bioperl-live even more problems occur:
bp_bioflat_index.pl --dbname sw -i bdb -f swiss -l . -c uniprot_sprot.dat
------------- EXCEPTION -------------
MSG: The lineage 'Eukaryota, Metazoa, Chordata, Craniata, Vertebrata,
Euteleostomi, Amphibia, Batrachia, Anura, Mesobatrachia, Pipoidea,
Pipidae, Xenopodinae, Xenopus, Silurana, Xenopus, tropicalis' had two
non-consecutive nodes with the same name. Can't cope!
But I think this is connected to the new changes to taxonomy handling in
I'm unsure wether to submit this separately, but I could also provide an
example of such a swissprot entry that causes this error.
Brian Osborne wrote:
> Bug, presumably in SeqIO/swiss.pm. Can you send me a small file with such a
> PRELIMINARY entry?
> Brian O.
> On 9/1/06 6:11 AM, "Daniel Lang" <daniel.lang at biologie.uni-freiburg.de>
>> when using Bio::Registry (bioperl-live) to fetch uniprot entries from
>> local indexed uniprot *.dats, I had to realize that several entries
>> could not be retrieved despite the fact that they are present in the
>> files! A closer look reveals that they are of status PRELIMINARY:
>> uniprot_trembl.dat:ID Q16EZ1_AEDAE PRELIMINARY; PRT; 222 AA.
>> I don't "grep" PRELIMINARY anywhere in my cvs checkout..
>> I also can't retrieve the sequences from the online database defined as
>> Is this a bug or a feature? If its a feature, how can I bypass it?
>> Thanks in advance,
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
More information about the Bioperl-l