import numpy as np

T = int(input())
MOD = 10 ** 9 + 7

M = np.array(
    [
        [1, 0, 0, 0, 1, 1],
        [0, 1, 0, 1, 0, 1],
        [0, 0, 1, 1, 1, 0],
        [3, 0, 0, 0, 0, 0],
        [0, 3, 0, 0, 0, 0],
        [0, 0, 3, 0, 0, 0],
    ]
)

doubling = [M]
for i in range(40):
    doubling.append(doubling[-1] @ doubling[-1])
    doubling[-1] = np.mod(doubling[-1], MOD)

fraction = pow(3, MOD - 2, MOD)

for _ in range(T):
    N = int(input())
    V = np.array([1, 0, 0, 3, 0, 0])

    for i in range(40)[::-1]:
        if N >> i & 1:
            V = doubling[i] @ V
            V = np.mod(V, MOD)

    print(V[3] * pow(fraction, N + 1, MOD) % MOD)