function d=distH(F) % distH.m % % usage: d=distH(F) % % F is a 4x4 frequency array from comparing two sequences % (e.g., by using compseq.m). Output is Hammingdistance % between sequences. % % 7/3//03 A=F-diag(diag(F)); d=sum(sum(A))/sum(sum(F));