[EMBOSS] Get a partial GenBank of a circular DNAþ
Peter Rice
pmr at ebi.ac.uk
Thu Jun 17 08:38:21 UTC 2010
On 16/06/2010 13:55, Roi Brodo wrote:
>
> After some more reading I think I can do it using union. The problem is that after I create the list (of the two ranges) using yank, union dies on "union terminated: Bad value for '-sequence' and no prompt". Why is that? Shouldn't i use a yank file?
Yes, yank and union is the correct approach.
The output of yank is a list file, so the input to union should be
@filename to read a list of sequence addresses from the file.
If you just give the filename it assumes it is sequences (perhaps a
fasta file of sequences to be joined).
We will add this to our feature requests - it should be possible to make
seqret handle ranges from circular sequences. This will be after the
next release as it requires rewriting the way several library functions
work to allow the circular range.
regards,
Peter Rice
More information about the EMBOSS
mailing list