function rf = plotrecfracs(geno, chrid)
%genome coverage plot
rf = est_rec_fracs(geno);
m = size(rf,1)
temp = [[rf,zeros(m,1)];zeros(1,m+1)];
h = surf(temp)
view(2)
axis([1 m+1 1 m+1])
shading flat
title('Recombination Fractions')
colormap(flipud(jet))
colorbar

if nargin > 1
   n = max(chrid);
   mp = zeros(1,n)
   for i = 1:n
      mp(i) = mean(find(chrid==i));
   end
   xlabel('Chromosome')
   ylabel('Chromosome')
   set(gca,'Ticklength',[0 0])
   set(gca,'XTick',mp)
   set(gca,'XTickLabel',[num2str([1:19]');' X'])
   set(gca,'YTick',mp)
   set(gca,'YTickLabel',[num2str([1:19]');' X'])
end