H, W = map(int, input().split())
table = {}
for i in range(H):
    for j, v in enumerate(list(input())):
        if v != '?':
            if i not in table:
                table[i] = {}
            table[i][j] = int(v)

exp1= H
exp2= W
v_chk = [-1] * W
for i, row in table.items():
    h_chk = -1
    for j, v in row.items():
        if exp1 != -1:
            if h_chk == -1:
                h_chk = (j + v) % 2
                exp1 -= 1
            else:
                if h_chk != (j + v) % 2:
                    exp1 = -1
        if exp2 != -1:
            if v_chk[j] == -1:
                v_chk[j] = (i + v) % 2
                exp2 -=1
            else:
                if v_chk[j] != (i + v) % 2:
                    exp2 = -1
    if exp1 == -1 and exp2 == -2:
        break

MOD = 1000000007
ans = 0
if exp1 != -1:
    ans += pow(2, exp1, MOD)
if exp2 != -1:
    if exp1 == -1:
        dup = 0
    elif exp2 == W:
        dup = 2
    else:
        dup = 1
        for j in range(W-1):
            if v_chk[j] != -1 and v_chk[j] == v_chk[j+1]:
                dup = 0
                break
    ans += pow(2, exp2, MOD) - dup
    ans %= MOD
print(ans)