Evolutionary algorithm to generate offspring of one or morph two (peptide) sequences into each other.
Generate offspring for a given sequence. This script is designed to work with peptide sequences. However, any sequence can be provided, given a suitable distance matrix is referenced.
python pepevolve.py <parent> --lambd <int> --sigma <float> --matrixfile <file> --skip <str> --seed <int>
parent
: {str} parent string from which you wish to generate childrenlambd
: {int} number of offsprings to generate for the given parentsigma
: {float} width of the gaussian distribution for tuning the distance to the parentmatrixfile
: {str} filename of the distance matrix to useskip
: {str} letters (AA) to skip when samplingseed
: {int} random seed to use when sampling (makes runs reproducible)
python pepevolve.py GLFDIVKKVVGALGSL --lambd 10 --sigma 0.1 --matrixfile grantham.txt --skip CM --seed 42
generated sequences with corresponding distances and sigma values written to the file restult.txt
Evolutionary algorithm to morph two given sequences into each other with the help of an evolutionary algorithm. This script is designed to work with peptide sequences. However, any sequence can be provided, given a suitable distance matrix is referenced.
python pepmorph.py <start> <target> --lambd <int> --sigma <float> --matrixfile <file> --skip <str> --seed <int>
parent
: {str} parent string from which you wish to generate childrentarget
: {str} target string to which the parent should be morphed to (same length as parent!)lambd
: {int} number of offsprings to generate in between the two sequencessigma
: {float} width of the gaussian distribution for tuning the distance to the previous sequencematrixfile
: {str} filename of the distance matrix to useskip
: {str} letters (AA) to skip when samplingseed
: {int} random seed to use when sampling (makes runs reproducible)
python pepmorph.py KLLKLLKKLLKLLK GLFDIVKKVVGALGSL --lambd 10 --sigma 0.1 --matrixfile grantham.txt --skip CM --seed 42
generated sequences with corresponding distances and sigma values written to the file restult.txt