ENH: Add n_jobs to report

Still need to parallelize MRI.
list(chain(*htmls)) -> sum(htmls, [])

and get rid of chain. Same thing below.