function p_OME = OMEprob(TA_GT1, TA_GT2, TA_FN, TA_sN)

% Each TA_GT should be 1x2 matrix

FT1 = TA_GT1(1);   % Family type of cell 1
FT2 = TA_GT2(1);

ST1 = TA_GT1(2);   % Subtype of cell 1
ST2 = TA_GT2(2);

FamDist = GeneDist(FT1, FT2  ,  TA_FN);

if FamDist >= 2
    p_OME = 0;
elseif FamDist == 1

    ST_Dist = GeneDist(ST1, ST2, TA_sN);
    p_OME = 0.2*exp(-0.1*ST_Dist);

else

    ST_Dist = GeneDist(ST1, ST2, TA_sN);
    p_OME = 0.95*exp(-0.1*ST_Dist);

end



end

