MutateSeq function do not mark saturated sites correctly

[javieriserte]

I think that array named "vector" in the code below is used to mark those sites that has a large number of non- synonymous mutations, however the marked site is a different randomly chosen site.

# Line 713 of tools.c
while(nsyn_persite[p]/Options.runs > print_parameter ) {
    # site p is saturated if code enter the while loop
    r=rndu(); 

    if(r!=1.0)
        p=r*protsize; 
    else
        p=0; 
    # A new site p is chosen randomly.

    vector[p]=1; 
    # The new site p is marked as saturated, not the initial one.

    cont=0; 
    for(i=0;i<protsize;i++) { 
        if(vector[i]!=0){ 
            cont++; 
        } 
    } 
    if(cont==protsize) { 
        average=1.; 
        return(1); 
    } 
}