[EMBOSS] Does emboss have a handy way for mutate a protein sequence?

Peter Rice pmr at ebi.ac.uk
Tue Mar 7 13:55:39 UTC 2006


Zhiqiang Ye wrote:

> Dear all,
> 
>       Does emboss have a handy way for mutate a protein sequence by
> the specified way?
> For example, I have a sequence foo.fasta
> 
> 
>>foo
> 
> MATSCGLLKIIQRE
> 
> It has a mutant called  'A2L'. Is there any way to do this operation
> to output(with an option to check the foo.fasta has 'A' at position
> 2):
> 
>>foo A2L
> 
> MLTSCGLLKIIQRE

EMBOSS has several programs to change sequences. None does exactly what you ask.

You could look at:

biosed (does what you ask for longer replacements, but will change all 'A's to 
'L's.)

We could extend biosed to specify the position of the pattern ... is that what 
you need?

regards,

Peter










More information about the EMBOSS mailing list