Numerical inverse Laplace transform

When using the Inverse Laplace Transform at a specific point, what is the default method for numerical integration if the function is monotonically decreasing with a singularity only at the origin?

You may try an experiment:

results = 
 InverseLaplaceTransform[1/(Sqrt[s + 1] + s), s, 5.3, 
    Method -> #] & /@ {Automatic, "Stehfest", "Piessens", "Talbot", 
   "Durbin", "Crump", "Papoulis", "Weeks"}
results - results[[1]]

It would seem that for this function the chosen method was Stehfest. From this

results - InverseLaplaceTransform[1/(Sqrt[s + 1] + s), s, 53/10]

it seems that the best method here is Talbot, but the difference from Stehfest is around MachinePrecision.