seqret reads and writes (returns) sequences. Here is a sample session with seqret. It is used to extract an entry from a database and write it to a file: % seqret Input sequence: embl:paamir Output sequence [paamir.fasta]: Here seqret is used to display the contents of the sequence on the screen: % seqret Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.fasta]: stdout Here it is used in three different ways to write the result to a file in GCG format. Once by using the qualifier -osformat, once by using the format in the output USA on the command line, and once by specifying the format in the USA at the prompt: % seqret -osf gcg Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.gcg]: % seqret -outseq gcg::paamir.gcg Reads and writes (returns) sequences Input sequence: embl:paamir % seqret Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.fasta]: gcg::paamir.gcg Here seqret is used to produce the reverse-complement of a sequence: % seqret -srev Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.fasta]: Here seqret is used to extract the bases between the positions starting at 5 and ending at 25: % seqret -sbegin 5 -send 25 Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.fasta]: Here seqret is used to extract the bases between the positions starting at 5 and ending at 5 bases before the end of the sequence: % seqret -sbegin 5 -send -5 Reads and writes (returns) sequences Input sequence: embl:paamir Output sequence [paamir.fasta]: Here seqret is used to read all entries in the database tembl that start with "hs" and write them to a file: % seqret Reads and writes (returns) sequences Input sequence(s): embl:hs* Output sequence [hs989235.fasta]: Here seqret is used to read all entries in the database tembl that start with "hs" and writes them to a file. In this example, the specification is all done in the command line. To avoid confusing Unix with the * character, put backslash (\) before it: % seqret embl:hs\* hs989235.fasta Reads and writes (returns) sequences Here seqret is used to read only the entry "hsfau" from the file hs989235.fasta, which contains many entries: % seqret Reads and writes (returns) sequences Input sequence(s): hs989235.fasta:hsfau Output sequence [hsfau.fasta]: Here seqret is used to read all entries in the file hs989235.fasta, but only writes the first one of these entries out to a file: % seqret -firstonly Reads and writes (returns) sequences Input sequence(s): hs989235.fasta Output sequence [hs989235.fasta]: first.fasta Here seqret is used to display the short sequence actgatcgtg in uppercase EMBL format on the screen: % seqret -supper -osf embl asis::actgatcgtg stdout Reads and writes (returns) sequences ID standard; DNA; UNC; 10 BP. SQ Sequence 10 BP; 2 A; 2 C; 3 G; 3 T; 0 other; ACTGATCGTG 10 // To force seqret to both read in and write out features, use the command-line qualifier -feature. seqret does not read in features by default, resulting in slightly faster performance. If you want to read features with your sequence and write them out on output, using -feature will change the default behavior so that any features present in the sequence are used: % seqret -feature Reads and writes (returns) sequences Input sequence(s): em:hs989235 Output sequence [hs989235.fasta]: embl::hs989235.embl Mandatory qualifiers:
Advanced qualifiers:
|