[EMBOSS] Seqret fetch by URL
Steve Taylor
stephen.taylor at imm.ox.ac.uk
Wed Jun 24 10:04:17 UTC 2009
Hi,
I am trying to fetch a sequence using seqret. In my EMBOSS defaults I have the following entry set up. I am using EMBOSS 6.0.1 both on Linux and Solaris.
DB NCBI_prot [ method: url format: fasta type: P
url: "http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi?db=protein&qty=1&c_start=1&list_uids=%s&uids=&dopt=fasta&dispmax=5&sendto=t&fr
om=begin&to=end"
comment: "NCBI protein - FASTA seqs only"
]
If I try:
seqret NCBI_prot:ACA69081
I get:
Reads and writes (returns) sequences
Error: Unable to read sequence 'NCBI_prot:ACA69081'
Died: seqret terminated: Bad value for '-sequence' and no prompt
However, if I try
http://www.ncbi.nlm.nih.gov/sviewer/viewer.fcgi?db=protein&qty=1&c_start=1&list_uids=ACA69081&uids=&dopt=fasta&dispmax=5&sendto=t&from=begin&to=end
I get
>gi|169751563|gb|ACA69081.1| glycoside hydrolase family 1 [Yersinia pseudotuberculosis YPIII]
MSYQQLPKDFLWGGAVAAHQVEGGWDKGGKGVSIADVLSGGSHGVDRVMTDGVLEGYRYPNHEAVDFYSH
YKEDIALFAEMGFKCFRTSIAWTRIFPHGDEQQPNEAGLQFYDDMFDELLKYGIEPVITLSHFEMPWHLV
KEYGGWKNRKVVDFFVKFSEVVMARYKSKVKYWMTFNEINNQRNWKYPLFGYCCSGVVFTEQENPEETLY
QVLHHQFVASAKVVKLGHAINPEFKIGCMVAMVPLYPFSCHPDDMMYSVEAMRERYLFGDVHMRGYYPSY
ILQEWARRGFNIHMEEGDLETLRDGCADYMGLSYYMSNAVSAINPGSGNSLSGFEGSVPNPHVKASDWGW
QIDPVGLRYSLSVLYERYQKPLFIVENGFGAIDKVAADGMVHDDYRIAYLKAHIEQMKKAVFEDGVDLMG
YTPWGCIDCVSFTTGEYSKRYGFIYVDKNDDGTGTMARSRKLSFDWYKKVIASNGEVL
Any ideas why or is there a way I can get more verbose information about what is going on?
Thanks,
Steve
More information about the EMBOSS
mailing list