We performed 2D PIC simulations of Kelvin Helmholtz instability (KHI) with symmetric and asymmetric density and temperature profiles along the flow shear with a northward interplanetary magnetic field. The Magnetic Flux Transport method, field topology and magnetic field minimums are used to identify the reconnection X-lines. We start to observe the reconnection signatures such as magnetic field and flow reversals at the vortex edges in the nonlinear phase of the KHI when the vortices are rolling up. The number of reconnection regions increases at the turbulence phase. The signatures eventually decrease and finally disappear at very turbulent stages of KHI developments. Our results qualitatively agree with MMS observations of reconnection signatures at KHI along the magnetospheric flanks.