Why can you solve LASSO with gradient descent if LASSO is not differentiable

We know that LASSO is

min ||Ax - b||_2² + lambda||beta||_1

where

the L1 regularization is not differentiable due to the absolute operator

so I am really confused to how we use gradient descent to solve LASSO

Thanks