10 MAXFILES=3
11 DIM C$(4000), C1$(4000)
12 tx=160:c=0:gs=0:o2=0:o1=0:q=0
20 INPUT"file 1";F1$:OPEN F1$+".grd" FOR INPUT AS #1
30 INPUT "file 2";F2$:OPEN F2$+".grd" FOR INPUT AS #2
40 input"crit-low";cl
41 input"crit-high";cu
42 OPEN "cmp" FOR OUTPUT AS #3:CH=0
50 PRINT"":IF NOT EOF(2) THEN FOR I%=0 TO TX+1:INPUT#2,C$(I%):NEXT I%:o2=o2+1:ELSE PRINT#3,"chi square =" CH:PRINT "chi square =";CH/o1:print "organisms=";o1/o2,o2,o1/o2+o2:CLOSE:BEEP:SYSTEM
70 IF NOT EOF(1) THEN FOR I%=0 TO TX+1:INPUT#1,C1$(I%):NEXT I%:o1=o1+1:ut=0:ELSE CLOSE 1:OPEN F1$+".grd" FOR INPUT AS #1: PRINT"   "F1$" vs "C$(0);"   ";q/(o1/o2):q=0:GOTO 50
80 T=0:rem IF C$(0)=C1$(0) THEN GOTO 70
81 FOR Z%=0 TO TX+1
82 IF C$(Z%)="" or C$(Z%)="_" OR C1$(Z%)="" or C1$(Z%)="_" THEN GOTO 100
83 if z%=69 or z%=70 or z%=67 or z%=68 or z%=71 or z%=72  or z%=143 or z%=142 then goto 100
84 UT=UT+1:gosub 2000
100 NEXT Z%:IF (UT=0)THEN PRINT"x";:c=0:GOTO 114 
101 C=T/UT:if c >= cu then print "*";:gosub 1000:goto 114
110 IF C >= cl THEN PRINT "!";:gosub 1000:ELSE PRINT ".";
114 CH=CH+ C:q=q+c:GOTO 70
1000 PRINT#3,c","C$(0)","C1$(0)","T","UT","F2$","F1$:gs=gs+1:return
2000 if inSTR(C$(Z%),C1$(Z%)) or C$(Z%)="D" or C1$(Z%)="D" THEN T=T+1:return
2010 if (z%=159 or z%=160 or z%=95) and (c$(z%)<>"0" or c1$(z%)<>"0") then goto 3000 :else return
3000 if abs( val(c$(z%)) - val(c1$(z%)) )  <=.1 then t=t+1
3010 return