function [start_position_array] = find_genJ_start_positions(sequence_data_file, genJ, num_nuc_cmp);

[data, header] = load_data(sequence_data_file);

Sequence = data.Sequence;

clear data header

max_genJ_seql = 0;
for genJ_in = 1:length(genJ)
    if length(genJ(genJ_in).seq)>max_genJ_seql
        max_genJ_seql = length(genJ(genJ_in).seq);
    end
end

start_position_array = zeros(length(genJ), max_genJ_seql);

for genJ_in = 1:length(genJ)
    current_genJ_seq = genJ(genJ_in).seq;
    disp(genJ_in);
    for start_pos_in = num_nuc_cmp:length(genJ(genJ_in).seq)
        for i = 1:length(Sequence)
            if strcmp(current_genJ_seq((start_pos_in - num_nuc_cmp + 1) : start_pos_in), Sequence{i}((end-num_nuc_cmp +1):end))
                start_position_array(genJ_in, start_pos_in) = start_position_array(genJ_in, start_pos_in) + 1;
            end
        end
    end
end




end