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

[David.Bauer at SCHERING.DE]

What about this solution:

cutseq foo.fasta -from 2 -to 2 | pasteseq -filter -pos 1 -bs asis:'L' |
descseq -filter -append -desc "A2L"

>foo A2L
MLTSCGLLKIIQRE

Cheers,
David.


emboss-bounces at emboss.open-bio.org schrieb am 06/03/2006 22:29:49:

> 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
>
> My way:  use extractseq to extract two file: one before position 2,
> the other after postion 2. Then creat a fasta file contain 'L'.  After
> that,  I use union to connect these 3 sequence file in to one.
>
> Or write a perl script to do this by change a string's substring.
>
> How If emboss could provide a 'mutate' !
>
> Thank you :)
>
> Best regards!
> --
> Zhiqiang Ye
>
> _______________________________________________
> EMBOSS mailing list
> EMBOSS at emboss.open-bio.org
> http://newportal.open-bio.org/mailman/listinfo/emboss