I did a same for Wood (I forget exact name of software, need to check it). Ugopen application, for all cylindrical faces output xyz for center point and ijk for tool axis. Later I did a same for riveting machine.
I remember, it was WoodWop software for woodworking CNC Macine. Exact your problem, I recognize all cylindrical holes and output xyz for center and ijk vector for tool axis.. See a video.
Or - you can insert few point-to-point CAM operations with defined tool axis for each point and output as CLS, each line will be GOTO/X,Y,Z,I,J,K as you need..