Linkage disequilibrium (LD) mapping is able to localize quantitative trait loci (QTL) within a rather small region (e.g. 2 cM), which is much narrower than linkage analysis (LA, usually 20 cM). The multilocus LD method utilizes haplotype information around putative mutation and takes historical recombination events into account, and thus provides a powerful method for further fine mapping. However, sometimes there are more than one QTLs in the region being studied. In this study, the Bayesian model selection implemented via the Markov chain Monte Carlo (MCMC) method is developed for fine mapping of multiple QTLs using haplotype information in a small region. The method combines LD as well as linkage information. A series of simulation experiments were conducted to investigate the behavior of the method. The results showed that this new multiple QTLs method was more efficient in separating closely linked QTLs than single-marker association studies.