I want to put on a QGIS map some historical sites to check if they are aligned. If an alignement exists, the points should stay on a great circle (is this correct?). I want to visually verify this using a map projection where orthodromic lines are represented as straight lines. In a second step I would draw two boundary lines and calculate bearings from a determinate point.
For example, here is a snapshot of a QGIS map with the points of "The St Michael-Apollo 'Axis'" as described by Lucien Richer, 1977.
Map boundaries: (-11E, 32N) (36E, 52N). EPSG 3785 projection (popular visualization / Mercator) .
What would be the most appropriate map projection to get great circles as straight lines?
Do I need a custom projection? If yes, what are the parameters to use in QGIS?
Any suggestions would be helpful.