dbiblast

Qingjuan Gu qingjuan at wistar.upenn.edu
Thu Jul 25 16:17:08 UTC 2002


>The answer is that you should use 'chromosome' as the wildcard 
>database name. Secondly, the author of dbiblast thinks there 
>might be some problems with the 2.4.1 version and this will 
>be corrected in 2.5.0 which is out soon. Thirdly, do not use 
>the very latest formatdb (NCBI have not yet released the 
>details of their index format so EMBOSS cannot index them). 

I did not use 'formatdb' for blast database, I just download from ncbi ftpsite "ftp://ftp.ncbi.nih.gov/blast/db/FormattedDatabases/" download and use gunzip and untar command to get these database. 

I follow you instruction, it failed also. Here is the interactive result: 
# ls -l 
total 434903 
-rwxr--r--   1 nobody   blast     76626843 Jul  4 22:09 patnt.nhr 
-rwxr--r--   1 nobody   blast      8066892 Jul  4 22:19 patnt.nin 
-rwxr--r--   1 nobody   blast      5377840 Jul  4 22:10 patnt.nnd 
-rwxr--r--   1 nobody   blast        21052 Jul  4 22:10 patnt.nni 
-rwxr--r--   1 nobody   blast    142820076 Jul  4 22:13 patnt.nsd 
-rwxr--r--   1 nobody   blast      2926655 Jul  4 22:19 patnt.nsi 
-rwxr--r--   1 nobody   blast     92506630 Jul  4 22:19 patnt.nsq 
-rwxr--r--   1 nobody   blast       119501 Jul  4 22:19 patnt.ntd 
-rwxr--r--   1 nobody   blast        17904 Jul  4 22:19 patnt.nti 
-rwxr--r--   1 nobody   blast    116819853 Jun 28 09:30 patnt.tar.gz 
# 
# dbiblast 
Index a BLAST database 
Database name: patnt 
Database directory [.]: 
Wildcard database filename [patnt]: patnt 
Release number [0.0]: 
Index date [00/00/00]: 
         N : nucleic 
         P : protein 
         ? : unknown 
Sequence type [unknown]: N 
         1 : wublast and setdb/pressdb 
         2 : formatdb 
         0 : unknown 
Blast index version [unknown]: 2 

   EMBOSS An error in dbiblast.c at line 333: 
No files selected 
# 
if i use patnt.nin in 'wildcard databse filename', that works without complain. but when I use seqret or other command to search the sequence, there are always error message: 

# ls -l 
total 435242 
-rwxr-xr-x   1 root     system         416 Jul 25 10:58 acnum.hit 
-rwxr-xr-x   1 root     system         742 Jul 25 10:58 acnum.trg 
-rwxr-xr-x   1 root     system         322 Jul 25 10:56 division.lkp 
-rwxr-xr-x   1 root     system      338334 Jul 25 10:58 entrynam.idx 
-rwxr--r--   1 nobody   blast     76626843 Jul  4 22:09 patnt.nhr 
-rwxr--r--   1 nobody   blast      8066892 Jul  4 22:19 patnt.nin 
-rwxr--r--   1 nobody   blast      5377840 Jul  4 22:10 patnt.nnd 
-rwxr--r--   1 nobody   blast        21052 Jul  4 22:10 patnt.nni 
-rwxr--r--   1 nobody   blast    142820076 Jul  4 22:13 patnt.nsd 
-rwxr--r--   1 nobody   blast      2926655 Jul  4 22:19 patnt.nsi 
-rwxr--r--   1 nobody   blast     92506630 Jul  4 22:19 patnt.nsq 
-rwxr--r--   1 nobody   blast       119501 Jul  4 22:19 patnt.ntd 
-rwxr--r--   1 nobody   blast        17904 Jul  4 22:19 patnt.nti 
-rwxr--r--   1 nobody   blast    116819853 Jun 28 09:30 patnt.tar.gz 
# seqret 
Reads and writes (returns) sequences 
Input sequence(s): patnt:E00007.1 
Error: BLAST Entry failed 
Error: Unable to read sequence 'patnt:E00007.1' 
Input sequence(s): patnt:E00007 
Error: BLAST Entry failed 
Error: Unable to read sequence 'patnt:E00007' 
# seqret 
Reads and writes (returns) sequences 
Input sequence(s): patnt:2168318 
Error: BLAST Entry failed 
Error: Unable to read sequence 'patnt:2168318' 
Input sequence(s): patnt:2168314 
Error: BLAST Entry failed 
Error: Unable to read sequence 'patnt:2168314' 
# 

In my emboss.default:
DB patnt [
        type: N
        method: blast
        format: ncbi
        dir: /db4/users/qgu/test/test_dbi
        file: patnt.*
        indexdir: /db4/users/qgu/test/test_dbi
        release: "1.0"
        comment: "patnt formatted database --BLAST NCBI"
]

Also, my fasta and genbank files works fine, only blast database can not work well. 

Your help will be highly appreciated. 

Grace



More information about the EMBOSS mailing list