HWE.test=function(als, counts=apply(als, 1, sum)){
# als is an nx2 matrix of allele counts
	no.als=2*nrow(als)
	t1=table(counts)
	p=(2*t1[1]+t1[2])/no.als
	cat("p",p,"\n")
	e0=p^2*no.als; e1=2*p*(1-p)*no.als; e2=(1-p)^2*no.als
	es=c(e0, e1, e2)
	chi.sq=sum((es-t1)^2/(es))
	pchisq.new=1-pchisq(chi.sq, 1)
	return(pchisq.new)
	}