[EMBOSS] Using seqret to fetch from .nal index databases

Audra Johnson john8376 at uidaho.edu
Fri Jul 29 19:08:50 UTC 2005


Apologies for the length, but I want to be thorough.  I'm doing blast  
searches and then trying to fetch the sequences from the our genembl  
database using seqret.  For example:

blastall -p tblastn /gcgdata_10.3/gcgblast/genembl -i  
dp00061_disordered_115_168.fasta

Gives me results of:

GB_PR:HUMRPA70KD        2e-08   412     573     1       54      54
GB_PR:BC018126  2e-08   386     547     1       54      54
GB_PAT:AX335048 2e-08   412     573     1       54      54
GB_PAT:AR175924 2e-08   412     573     1       54      54
GB_RO:BC019119  0.003   399     584     1       53      62

I've tried using a seqret just for the database name I'm giving  
blastall, and specifically saying the genembl.nal file:

$ seqret
Reads and writes (returns) sequences
Input sequence(s): /gcgdata_10.3/gcgblast/genembl.nal:HUMRPA70KD
Error: Unable to read sequence '/gcgdata_10.3/gcgblast/ 
genembl.nal:HUMRPA70KD'
Input sequence(s): /gcgdata_10.3/gcgblast/genembl:HUMRPA70KD
Error: failed to open filename '/gcgdata_10.3/gcgblast/genembl'
Error: Unable to read sequence '/gcgdata_10.3/gcgblast/ 
genembl:HUMRPA70KD'
Died: seqret terminated: Bad value for '-sequence' and no more retries

But neither works.  (I've omitted the beginning prefix GB_PR: and  
similar prefixes, but I've tried that way and it doesn't work,  
either.)  Is there any way to get seqret functioning with these  
databases?

-- Audra Johnson, University of Idaho



More information about the EMBOSS mailing list