Background: Inference of evolutionary trees using the maximum likelihood principle is NP-hard. Therefore, all practical methods rely on heuristics. The topological transformations often used in heuristics are Nearest Neighbor Interchange (NNI), Subtree Prune and Regraft (SPR) and Tree Bisection and Reconnection (TBR). However, these topological transformations often fall easily into local optima, since there are not many trees accessible in one step from any given tree. Another more exhaustive topological transformation is p-Edge Contraction and Refinement (pECR). However, due to its high computation complexity, p-ECR has rarely been used in practice. Results: To make the p-ECR move more efficient, this paper proposes a new method named pECRNJ. The main idea of p-ECRNJ is to use neighbor joining (NJ) to refine the unresolved nodes produced in p-ECR. Conclusion: Experiments with real datasets show that p-ECRNJ can find better trees than the best known maximum likelihood methods so far ...