function [geneIndex] = getGeneIndex(id, locusIDs)
N= length(locusIDs);
geneIndex=0;
for i=1:N
    if strcmp(locusIDs{i}, id)
        geneIndex=i;
        break;
    end
end