[EMBOSS] can't access databases indexed by dbifasta
    Marcus Claesson 
    m.claesson at student.ucc.ie
       
    Tue Feb 17 10:58:29 UTC 2004
    
    
  
Hello,
I have a silly little problem indexing databases in Emboss-2.8.0. After
running dbifasta and adding DB entries in emboss.default I can only
access the database when being in the same directory as the fasta file.
Here is what I did:
[blast_db]$ uname -a
Linux neo.ucc.ie 2.4.9-e.35enterprise #1 SMP Tue Dec 23 00:06:16 EST
2003 i686 unknown
[blast_db]$ pwd
/var/data/blast_db
[blast_db]$ ll ecoli.nt
-rw-r--r--    1 marcus   bioinfo   4763013 Jan 15 01:38 ecoli.nt
[blast_db]$ dbifasta
Index a fasta database
    simple : >ID
     idacc : >ID ACC
     gcgid : >db:ID
  gcgidacc : >db:ID ACC
      dbid : >db ID
      ncbi : | formats
ID line format [idacc]: 
Database directory [.]: /var/data/blast_db
Wildcard database filename [*.dat]: ecoli.nt
Database name: ecoli.nt
Release number [0.0]: 
Index date [00/00/00]:
[blast_db]$ ll entrynam.idx division.lkp acnum.*
-rw-rw-r--    1 marcus   bioinfo       300 Feb 17 10:39 acnum.hit
-rw-rw-r--    1 marcus   bioinfo       300 Feb 17 10:39 acnum.trg
-rw-rw-r--    1 marcus   bioinfo       330 Feb 17 10:39 division.lkp
-rw-rw-r--    1 marcus   bioinfo       300 Feb 17 10:39 entrynam.idx
Added these lines in /usr/local/EMBOSS-2.8.0/emboss/emboss.default:
DB ecoli.nt [
   type: "N"
   format: "fasta"
   method: "emblcd"
   dir: "/var/data/blast_db/"
]
[blast_db]$ showdb
Displays information on the currently available databases
# Name        Type ID  Qry All Comment
# ====        ==== ==  === === =======
ecoli.nt      N    OK  OK  OK  -
[blast_db]$ cd ~
[marcus]$ seqret ecoli.nt
Reads and writes (returns) sequences
Error: failed to open filename 'ecoli.nt'
Error: Unable to read sequence 'ecoli.nt'
Died: seqret terminated: Bad value for '-sequence' and no prompt
But it works when I'm the same directory as ecoli.nt:
[blast_db]$  seqret ecoli.nt
Reads and writes (returns) sequences
Output sequence [ae000111.fasta]: 
etc...
Clearly it must be possible to access ecoli.nt from other directories?
Extremly grateful for any help on this!
Regards,
Marcus
    
    
More information about the EMBOSS
mailing list