Exercise: seqret

unix % seqret
Reads and writes (returns) a sequence
Input sequence: tembl:xlrhodop
Output sequence [xlrhodop.fasta]:

unix % more xlrhodop.fasta

>XLRHODOP L07770 Xenopus laevis	rhodopsin mRNA,	complete cds.
ggtagaacagcttcagttgggatcacaggcttctagggatcctttgggcaaaaaagaaac
acagaaggcattctttctatacaagaaaggactttatagagctgctaccatgaacggaac

$\vdots$

Now let's retrieve the sequence using its the accession number:

unix % seqret
Reads and writes (returns) a sequence
Input sequence: tembl:L07770
Output sequence [xlrhodop.fasta]: xlrhodop2.fasta

unix % more xlrhodop2.fasta

>XLRHODOP L07770 Xenopus laevis	rhodopsin mRNA,	complete cds.
ggtagaacagcttcagttgggatcacaggcttctagggatcctttgggcaaaaaagaaac
acagaaggcattctttctatacaagaaaggactttatagagctgctaccatgaacggaac

$\vdots$

You could also run this example entirely from the command line:

unix % seqret tembl:xlrhodop -outseq xlrhodop.fasta

By default, seqret writes the sequence in fasta format. You can also tell it to use a different output format:

unix % seqret tembl:L07770 -outseq gcg::xlrhodop.gcg

unix % more xlrhodop.gcg

!!NA_SEQUENCE 1.0

Xenopus	laevis rhodopsin mRNA, complete	cds.

XLRHODOP  Length: 1684	Type: N	 Check:	9453 ..

    1 ggtagaacag cttcagttgg gatcacaggc ttctagggat cctttgggca

   51 aaaaagaaac acagaaggca ttctttctat acaagaaagg actttataga

$\vdots$

A list of the various formats that EMBOSS understands is given at
Usa/formats.html