function [i, j] = reportpairs(F1, F2, F3, F4, cv, mnames, chrid, fid)
%report signifcant hits from a pairwise genome scan
%	using criteria in Sugiyama et al, 1999
% Fi is an mxm array
% cv - is a 1x4 vector of critical values - may be empty
% n is the sample size
% fid is a filepointer to the output file name
% 	these are user specified but see pairpermute.m

[i, j] = find(F1>cv(1) & (F2>cv(2) | (F3>cv(3) & F4>cv(4))));
[temp, ord] = sortrows([chrid(i) chrid(j)]);
if length(i) > 0
   for s = 1:length(i)
      fprintf(fid,'%d\t%d\t%s\t%d\t%d\t%s\t%6.2f\t%6.2f\t%6.2f\t%6.2f\r\n',...
         i(ord(s)), chrid(i(ord(s))), mnames{i(ord(s))},...
         j(ord(s)), chrid(j(ord(s))), mnames{j(ord(s))},...
      	F1(i(ord(s)), j(ord(s))), F2(i(ord(s)), j(ord(s))), F3(i(ord(s)), j(ord(s))), F4(i(ord(s)), j(ord(s))));
   end
end
fprintf(fid,'\r\n');

%fclose(fid)