We consider within QCD collinear factorization the process p + p → jet + jet + X, where two forward high-p T jets are produced with a large separation in rapidity ∆y (Mueller-Navelet jets). In this case the (calculable) hard part of the reaction receives large higher-order corrections ∼ α n s (∆y) n , which can be accounted for in the BFKL approach with next-to-leading logarithmic accuracy, including contributions ∼ α n s (∆y) n−1 . We calculate several observables related with this process, using the next-to-leading order jet vertices, recently calculated in the approximation of small aperture of the jet cone in the pseudorapidity-azimuthal angle plane. †