You still can do a lot. There is no perfect sampler, jitter always exists. Moreover, traditional systems include quantization which prevents (theoretically) the satisfaction of Shannon-Nyquist-Kotelnikov-Raabe... conditions.
If your system includes low-pass filtering, or/and the phenomenon you look at is very weak above 166 Hz, the variations are likely to be small enough for your purposes. I would propose to look at an harmonic behavior with the following three approaches:
- Suppose the periods are even, then perform Fourier analysis,
- Resample in a simple fashion at even samples (linear interpolation), then perform Fourier analysis,
- Use a Fourier tool adapted to irregular sampling, like the Lomb-Scargle periodogram.
If you detect no meaningful difference in the features you extract, then you are fine (for the time being). You even can, routinely, use these methods in parallel to detect a weird behavior when they differ.
If you see meaningful difference then you might want to invest into irregular sampling processing methods. Luckily, in theory, when you have irregular samples, and a sufficient precision in the timing, you can extract information beyond Nyquist, at the cost of more involved methods. There is a handful of methods for Fourier or time-frequency methods, FIR results, etc.