A mathematical model for tick-borne infections: a numerical study