Hi,
the distribution of parameters obtained from the non maximum-lhood solutions cannot be used for CI inference, as it is not representing evolutionary variance, but due to some variance in the optimization procedure.
The procedure you describe above is a parametric bootstrap procedure, but you could also use a non-parametric (block) bootstrap procedure (probably better).
L