Multi-interval subfactors and modularity of representations in conformal field theory