Count reads in features