Without a reprex or experience with pegas or haploFreq, there isn't much to go on here.
The docs suggest (haploFreq: Haplotype Frequencies With a Covariate in pegas: Population and Evolutionary Genetics Analysis System) the number of sequences in Gull need to match the number of values of fac, Gull.region, but your error message suggests these numbers don't match.