In the current work, we consider the inverse conductivity problem of recovering inclusion with one measurement. First, we use conformal mapping techniques for determining the location of the anomaly and estimating its size. We then get a good initial guess for quasi-Newton type method. The inverse problem is treated from the shape optimization point of view. We give a rigorous proof for the existence of the derivative of the state function and of shape functionals. We consider both least squares fitting and Kohn and Vogelius functionals. For the numerical implementation, we use a parameterization of shapes coupled with a boundary element method. Several numerical examples indicate the superiority of the Kohn and Vogelius functional over least squares fitting.