In this answer @RussellBorogove discussed the frequency of dockings of each space shuttle with the Mir space station. I've summarized those frequencies, total number of missions flown, and year of first and last flight below.
Not counting the Enterprise and the Challenger (lost early in its expected life) there are four shuttles that span the period over which Mir docking missions were carried out.
Of the 9 dockings, 7 were made by the Atlantis. I've done a rough calculation tossing nine balls into four equal bins repeatedly, and then calculated roughly how many times the maximum number of balls in a given bin would be $N_{max}$. The probability that 7 of 9 balls would end up in any one bin is less than 0.5%. So I think the high frequency of Atlantis-Mir docking is not likely to have been just a statistical fluke, although it's not impossible.
Were there any reasons why NASA would have used Atlantis for such a large fraction of all the Shuttle-Mir dockings?
Enterprise: 0 / 0 (1977-2012, model only)
Columbia: 0 / 28 (1981-2003, lost)
Challenger: 0 / 10 (1983-1986, lost)
Discovery: 1 / 34 (1984-2011)
Atlantis: 7 / 33 (1985-2011)
Endeavour: 1 / 25 (1992-2011)
With 1,000,000 tries:
Nmax %
---- ----
3 37.2
4 43.2
5 15.6
6 3.5
7 0.5
8 <0.1
import numpy as np
from collections import Counter
n_tries = 100000 # I used a million for the question
# shamelessly inefficient way to do this
bob = np.random.randint(4, size=9*n_tries).reshape(9, -1)
qq = [max(Counter(thing).values()) for thing in bob.T]
a,b = np.histogram(qq, bins=range(11))
for A, B in zip(a, b[:-1]):
print B, float(A)/a.sum()