docs now include return values for NB, simulate_experiment functions, and write_reads