dbiblast problem
john walshaw (JIC)
john.walshaw at bbsrc.ac.uk
Thu Jun 27 09:04:11 UTC 2002
I have experienced a similar problem, tring to index WU-BLAST-formatted databases with
dbiblast. However I got it to work properly with NCBI-BLAST-formatted databases.
John Walshaw, John Innes Centre, Norwich Research Park,
Colney, Norwich NR4 7UH, UK. +44(0)1603 450827
> -----Original Message-----
> From: David.Bauer at SCHERING.DE [mailto:David.Bauer at SCHERING.DE]
> Sent: 26 June 2002 12:32
> To: emboss at embnet.org
> Subject: dbiblast problem
>
>
> Hi,
>
> I observed the folowing problem when retrieving entries from
> blast databases
> formated with dbiblast:
>
> I get the CORRECT sequence if I specify the ID:
> --------------------------------------------------------------
> ---------------------------------------------------------------
> seqret -auto -stdout cgdb_nt:celsl2a
> >CELSL2A M27263 C.elegans trans-spliced leader 2 (SL2
> RNA-alpha) gene, 5' flank.
> seqret -auto -stdout cgdb_nt:celsl2b
> >CELSL2B M27264 C.elegans trans-spliced leader 2 (SL2
> RNA-beta) gene, 5' flank.
> ##############################################################
> ###########
>
> But I get the WRONG sequences if I specify the ACC:
> --------------------------------------------------------------
> ----------------------------------------------------------------
> seqret -auto -stdout cgdb_nt:M27263
> >CELSL2B M27264 C.elegans trans-spliced leader 2 (SL2
> RNA-beta) gene, 5' flank.
> seqret -auto -stdout cgdb_nt:M27264
> >CELSNTI L15302 C.elegans synaptotagmin I mRNA, complete cds
> and flanking
> regions.
> ##############################################################
> #############
>
> With fastacmd the headers look like this:
> --------------------------------------------------------------
> -------------------------------------------------------------------
> >gb|M27263|CELSL2A C.elegans trans-spliced leader 2 (SL2
> RNA-alpha) gene, 5'
> flank
> >gb|M27264|CELSL2B C.elegans trans-spliced leader 2 (SL2
> RNA-beta) gene, 5'
> flank
> ##############################################################
> ############
>
> Any ideas ?
>
> Thanks, David.
>
>
More information about the EMBOSS
mailing list