[EMBOSS] Help to build a motif for fzzpro

Guy Bottu gbottu at vub.ac.be
Fri Oct 19 09:55:22 UTC 2012


	Dear Aengus,

I have an idea of how to do it. You must of course complete the 
motif/pattern as much as possible because with just {DE}(4) you will 
find much to much unless you search only a few sequences. You must run 
the program fuzzpro with parameter -rformat=listfile. You will obtain as 
output an EMBOSS list file. You can then run fuzzpro again with as 
sequence input list::xxx (xxx the name of your file) and as pattern 
input @yyy where yyy contains :

 > basic_1
<[HKR].
 > basic_2
<x-[HKR].
 > basic_3
<x-x-[HKR].
 > basic_4
<x-x-x-[HKR].

Indeed, the input will contain only sequence fragments matching the 
pattern and the basic amino acid hence should be in one of the first 4 
positions.

	Regards,
	Guy Bottu,
	ex-collaborator of the Belgian EMBnet Node



More information about the EMBOSS mailing list