function run_mus_selection_with_gen(data_file, model_file,his)

if ~exist('his','var')
	his = false;
end;

rng('shuffle');

gen_seq = construct_local_aaseq_draw(model_file, 300000);
model_name = strsplit(model_file,'/');
model_name = model_name{end-2};

data_files = strsplit(ls (data_file));
data_seq ={};
for k=1:length(data_files)-1
	T = readtable(data_files{k},'FileType','text','HeaderLines',3,'delimiter',' ');
	data_seq = [data_seq ; T.(2)(cellfun(@isempty,strfind(T.(2),'*'))) ];
end;

data_seq = unique(data_seq);
data_name = strsplit(data_file,'/');
data_name = strtok(data_name{end},'.');
data_name = strrep(data_name,'*','_aggregate');

if his
     lnq_GD_minimal_alllen(['results/selection/' data_name '_' model_name '_min'],gen_seq,data_seq,his,100000);
else
     lnq_GD_minimal_alllen(['results/selection/' data_name '_' model_name '_min'],gen_seq,data_seq);
end;
