get taxa counts per d, works, intermediate. benchmarking