Text this: Modeling COVID-19 dynamics in the Basque Country: characterizing population immunity profile from 2020 to 2022