local function array(dimension, default_val) local n=dimension local m={}if default_val~=nil then m[1]={__index=function()return default_val end}end for i=2,n do m[i]={__index=function(p, k)local c=setmetatable({},m[i-1])rawset(p,k,c)return c end}end return setmetatable({},m[n])end
----
local MOD = math.floor(10^9)+7
local S = io.read()
local N = #S
local ti = array(1, 0)
local tio = array(1, 0)
tio[0] = 1
for i=1,N do
    local c = S:sub(i,i)
    if c == '?' then
        local oldsum = tio[0] + tio[1] + tio[2] + tio[3] + tio[4] + tio[5] + tio[6] + tio[7] + tio[8] + tio[9] + tio[10] + tio[11] + tio[12]
        ti[0] = (oldsum - tio[4] - tio[8] - tio[12]) % MOD
        ti[1] = (oldsum - tio[3] - tio[8] - tio[12]) % MOD
        ti[2] = (oldsum - tio[3] - tio[7] - tio[12]) % MOD
        ti[3] = (oldsum - tio[3] - tio[7] - tio[11]) % MOD
        ti[4] = (oldsum - tio[2] - tio[7] - tio[11]) % MOD
        ti[5] = (oldsum - tio[2] - tio[6] - tio[11]) % MOD
        ti[6] = (oldsum - tio[2] - tio[6] - tio[10]) % MOD
        ti[7] = (oldsum - tio[1] - tio[6] - tio[10]) % MOD
        ti[8] = (oldsum - tio[1] - tio[5] - tio[10]) % MOD
        ti[9] = (oldsum - tio[1] - tio[5] - tio[9]) % MOD
        ti[10] = (oldsum - tio[0] - tio[5] - tio[9]) % MOD
        ti[11] = (oldsum - tio[0] - tio[4] - tio[9]) % MOD
        ti[12] = (oldsum - tio[0] - tio[4] - tio[8]) % MOD
    else
        c = tonumber(c)
        ti[(0 + c) % 13] = tio[0]
        ti[(10 + c) % 13] = tio[1]
        ti[(20 + c) % 13] = tio[2]
        ti[(30 + c) % 13] = tio[3]
        ti[(40 + c) % 13] = tio[4]
        ti[(50 + c) % 13] = tio[5]
        ti[(60 + c) % 13] = tio[6]
        ti[(70 + c) % 13] = tio[7]
        ti[(80 + c) % 13] = tio[8]
        ti[(90 + c) % 13] = tio[9]
        ti[(100 + c) % 13] = tio[10]
        ti[(110 + c) % 13] = tio[11]
        ti[(120 + c) % 13] = tio[12]
    end
    tio, ti = ti, tio
end
print(tio[5])