Although isogeometric analysis has shown great potential in achieving highly accurate numerical solutions of partial differential equations, how to efficiently implement the method is one of the challenges that makes it more competitive in practical simulations. In this paper, an integration of isogeometric analysis and a moving mesh method is proposed, providing a competitive approach to resolve the efficiency issue. Focusing on the Poisson equation, the implementation of the algorithm is presented in detail, including the numerical discretization of the governing equation using isogeometric analysis, and a mesh redistribution technique developed via harmonic maps. It is found that the isogeometric analysis brings attractive features in the realization of moving mesh method, such as it provides an accurate expression for moving direction of mesh nodes, and allows for more choices to construct monitor functions. Through a series of numerical experiments, the efficiency and effectiveness of the proposed method are successfully verified. Moreover, the potential of the method towards the practical applications is also well presented with the simulation of a helium atom in Kohn-Sham density functional theory.