We calculate angular correlation function between ultra-high energy cosmic rays (UHECR) observed by Yakutsk and AGASA experiments, and most powerful BL Lacertae objects. We find significant correlations which correspond to the probability of statistical fluctuation less than $10^{-4}$, including penatly for selecting the subset of brightest BL Lacs. We conclude that some of BL Lacs are sources of the observed UHECR and present a list of most probable candidates.