PROGRAM homocore

c	Calculates numbers of pariwise sequence identities

	CHARACTER*1	line(127,6)
	INTEGER		homo(6,6)

	OPEN (unit=1, form='FORMATTED', status='OLD', readonly,
	1	file='core.seq')

	DO n = 1,6
	  READ (1, 1000) (line(i,n), i=1,127)
1000	  FORMAT (127a1)
	endDO
	CLOSE (unit=1)

	DO i = 1,6
	  DO j = 1,6
	    homo(i,j) = 0
	    DO n = 1,127
	      IF (line(n,i) .EQ. line(n,j)) homo(i,j) = homo(i,j)+1
	    endDO
	  endDO
	endDO

	WRITE (6, 1001) homo
1001	FORMAT (6i6)
	end