Time-dependent probability distribution for number of infection in a stochastic SIS model: case study COVID-19