I think this has drifted slightly off topic. I am not looking to justify single run + error = accurate solution.
If the solution, to the problem on which you are working, is to be compared to a limiting stress then would it not be a reasonable approach to use the error estimation to estimate an upper bound stress value to ensure that the acceptance criteria is still satisfied. In suggesting this I am making the assumption that convergence has been demonstrated, mesh refinement is adequate and that the 'error' is 'low'. The solution will never be 100% accurate. If you have a number of stress categories to which stress limits apply then it seems logical, as a final operation, to quantify the maximum predicted stress. This would be particularly relevant where the predicted stresses are close to the stress limits.
I am also considering other software, such as Ansys, that can report upper bound values.