function [v,x,count,mutates] = NEW_find_mutate(mu,t,N)

 x = rand(t,N);
 VEC = ones(1,N);
 VEC0 = VEC;
 skip = [1 2 3];
 string = 'actg';

 v(1,:) = VEC0;
 count = 0;
 for i =1:t
  ii = i + 1;
  j = find(x(i,:)>1-mu);
  v(ii,:) = v(ii-1,:);
  for k = 1:length(j)
   a = randperm(3);
   v(ii,j(k)) = mod( v(ii,j(k)) + skip(a(1)),4);    
  end

  count = count + length(j); 
 end

 mutates=0;
 for i =1:N
  if v(ii,i) ~= VEC(i)
    mutates=mutates+1;
  end
 end
  
% mutates = sum(abs(v(ii,:)-VEC0));