The most general model-independent analysis of the rare $B$ decay, $Bsll$, is presented. There are ten independent local four-Fermi interactions which may contribute to this process. The branching ratio, the forward-backward asymmetry, and the double differential rate are written as functions of the Wilson coefficients of the ten operators. We also study the correlation between the branching ratio and the forward-backward asymmetry by changing each coefficient. This procedure tells us which types of operator contribute to the process, and it will be very useful to pin down new physics systematically, once we have the experimental data with high statistics and the deviation from the Standard Model is found.