[genV, genD, genJ] = load_Mus_TCR_genomic_segments();
seq = repmat('A',1,150);

cfg = coder.config('mex');
cfg.DynamicMemoryAllocation = 'AllVariableSizeArrays';

h=pwd();

%cd scripts
model = Mus_V_DJ_cuts_corrinserts_dinuc_ntbias_model(genV, genD, genJ);
model.max_assignments = 5000;
model.call_codegen(genV, genD, genJ);

cd(h)
