Study of Bs→ℓ+ℓ−γ decays in covariant quark model