Can we show that $\quad Pr(X_1>X_2) = Pr(X_2>X_1)$ ?
No, we can't. Counterexamples are simple to construct.
e.g. consider a bivariate distribution on $(0,1)^2$, and split the unit square up into 4x4 subregions, only some of which have any probability. Within sub-regions that have probability, make the density in each sub-region uniform (the same number need to be occupied in every "row" and "column"), such that they all have the same density. If you choose the "occupied" subregions intelligently, you can get the probabilities you ask about to be plainly unequal, as here:
$\qquad$ ![enter image description here](https://i.stack.imgur.com/kMVDV.png)
Clearly the margins are both U(0,1), yet the probabilities you mention are 3/4 and 1/4 by inspection (aside from single points, which have probability 0, the entirety of the 3 upper blocks is above the diagonal and the remaining block below it).
(The points I make in my demonstrations here could be made rigorous, but I think the intuition will be better conveyed if I stick to pictures and general reasoning.)
Just as plainly, this example applies to other continuous marginal distributions with montonic increasing cdf $F$, simply by applying a monotonic transformation $t(X) = F^{-1}(X)$ to these variates (which doesn't change the probabilities involved). [The existence of such simple counterexamples doesn't only apply to that case, it's just that the above argument generalizes to it especially easily.]
e.g. here's the above example transformed so that the margins are both standard normal:
![enter image description here](https://i.stack.imgur.com/7jXqB.png)
[plot code adapted from Ken Kleinman's function here]
By using a similar trick, you can make the probabilities arbitrarily far apart (with the restriction that they stay in (0,1) ).
For example, I could as easily split into a 100x100 subgrid and put 99 of the occupied subcells above the diagonal, and only one below it, for a 99%-1% split.
Here's a different bivariate distribution with uniform margins that I then transform to exponential:
![enter image description here](https://i.stack.imgur.com/dkAEe.png)
If you're happy to consider degenerate bivariate distributions, it's even easier to come up with examples. Choose $X_1$ uniform on (0,1) and let $X_2 = X_1+a \mod 1$, choosing $a$ to be any convenient small number.
![enter image description here](https://i.stack.imgur.com/gbLIy.png)
You may find it instructive to investigate copulas, which are multivariate distributions with uniform margins.
Best Answer
First things first. There needs to be greater information given as this does not have a universally correct answer. Different types of distributions have to be looked at with different types of procedures.
But just to show that yes this is possible, we assume that each of the variables that you have mentioned are normally distributed but the parameters of the normal distributions are different from each other for any given pair.
Now we take n samples each of these variables. Then calculate the correlation coefficients for each pair of the variables. If we cannot reject the hypothesis of these correlation coefficients being zero, we hypothesize that the variables are independent of each other. So we have a set of variables which are independent from each other, but they have different probability distributions.