reads are simulated with expression