function [occupancy] = motifOccupancy(motifs,loc,rloc)

[G,M]=size(loc);
occupancy=cell(G,M);
tic
for i=1:M
    for g=1:G
    %    fprintf('i=%d, g=%d\n',i,g);
        L =length(motifs{i});
        c1 = length(loc{g,i});
        c2 = length(rloc{g,i});
        occupancy{g,i} = cell(1,c1+c2);
        for r=1:c1
            occupancy{g,i}{r} = [loc{g,i}(r):1:(loc{g,i}(r)+L-1)];
        end
        for r=1:c2
            occupancy{g,i}{r+c1} =  [rloc{g,i}(r):1:(rloc{g,i}(r)+L-1)];
        end        

    end
end
toc