Performance of select from a 3d list – Mathematica slower than Python

I am creating a random 3d data set in Matematica 12.1. Then I am selecting all points that are in a certain range of one axis.

The same I am doing in Python (same computer, Python 3.8.5, numpy 1.19.2)

RESULT:
It seems that Python is able to select much faster (1.7 sec) than Mathematica (5.2 sec). What is the reason for that?
For selection in Mathematica I used the fastest solution, which is by Carl Woll (see here at bottom).

SeedRandom(1);
coordinates = RandomReal(10, {100000000, 3});

selectedCoordinates = 
   Pick(coordinates, 
    Unitize@Clip(coordinates((All, 1)), {6, 7}, {0, 0}), 
    1); // AbsoluteTiming

{5.16326, Null}

Dimensions(coordinates)

{100000000, 3}

Dimensions(selectedCoordinates)

{10003201, 3}

PYTHON CODE:

import time
import numpy as np
 
np.random.seed(1)
coordinates = np.random.random_sample((100000000,3))*10

start = time.time()
selectedCoordinates = coordinates((coordinates(:,0) > 6) & (coordinates(:,0) < 7))
end = time.time()

print(end-start)

print(coordinates.shape)

print(selectedCoordinates.shape)

1.6979997158050537

(100000000, 3)

(9997954, 3)