all_fits now parallelizes over both genes and regions instead of just over regions per gene