[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