% Dirk Husmeier
% August 2003
% This program plots the topology defining sites,
% reading in the files from TOADS

clf
clear

load topos.in
load topoDefSites.in

Nmax=max(topoDefSites);

set (gca, 'Fontsize', 15);
xTick4= Nmax-mod(Nmax,100);
xTick1= xTick4/4;
xTick2= xTick4/2;
xTick3= 3*xTick4/4;

plot(topoDefSites,topos,'b.','markersize',10);

axis([1 Nmax 0.5 3.5])
set(gca,'YTickMode','manual','YTick',[1 2 3],'XTickMode','manual','XTick',[xTick1 xTick2 xTick3 xTick4])
xlabel('Position in the alignment')
ylabel('Topology')