+ sigma++;
+ uint * occ = new uint[sigma+1];
+ for(uint i=0;i<=sigma;i++) occ[i] = 0;
+ for(uint i=0;i<n;i++) occ[am->map(seq[i])+1]++;
+ for(uint i=1;i<sigma;i++) occ[i] += occ[i-1];
+ uint * pos = new uint[n];
+ for(uint i=0;i<n;i++) pos[i] = 0;
+ for(uint i=0;i<n;i++) pos[occ[am->map(seq[i])]++]=i;
+ bitmaps = new static_bitsequence*[sigma];
+ uint * bm = new uint[uint_len(n,1)];
+ uint pp=0;