PROGRAM correctseq

c	Adjusts sequence to agrre with gene sequence

	CHARACTER*65	line

	WRITE (6, 1000) 'Enter WH input filename:'
1000	FORMAT (a)
	READ (5, 1000) line

	OPEN (unit=1, form='FORMATTED', status='OLD', readonly,
	1	file=line)

	WRITE (6, 1000) 'Enter WH output filename:'
	READ (5, 1000) line

	OPEN (unit=2, form='FORMATTED', status='NEW', file=line,
	1	carriagecontrol='LIST')

	DO WHILE (.TRUE.)
	  READ (1, 1000, end=9) line
	  READ (line(9:11), 1001) nres
1001	  FORMAT (i3)
	  IF (nres .EQ. 59 .OR. nres .EQ. 76) THEN
	    line(4:6) = 'ASN'
	    IF (line(8:8) .EQ. 'D') line(8:8) = 'N'
	    IF (line(12:14) .EQ. 'OD2') line(12:14) = 'ND2'
	  ELSE IF (nres .EQ. 517) THEN
	    line(4:6) = 'ASP'
	    IF (line(8:8) .EQ. 'N') line(8:8) = 'D'
	    IF (line(12:14) .EQ. 'ND2') line(12:14) = 'OD2'
	  endIF
	  WRITE (2, 1000) line
	endDO

9	CLOSE (unit=1)
	CLOSE (unit=2)

	end