Flux-ratio anomalies from discs and other baryonic structures in the Illustris simulation