#pragma GCC optimize("Ofast","inline")
#include<bits/stdc++.h>
using namespace std;
const int N=2e4,mod=1e9+7;

#define F(a,x) for(int a=0;a<x;a++)

int T,n=1e4,f[N][7][7][7],g[N][7][7][7],cm[N][8],h[3][7],ans,d[4][8][8];
char s[N];

void ad(int&x,int y)
{
    x+=y;
    x-=x>=mod?mod:0;
}

void init()
{
    f[0][0][0][0]=g[0][0][0][0]=1;
    for(int i=1;i<7;i++)cm[0][i]=1;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<7;j++)cm[i][j]=cm[i-1][j]*j%7;
        memcpy(f[i],f[i-1],sizeof f[i-1]);
        memcpy(g[i],g[i-1],sizeof g[i-1]);
        F(a,7)F(b,7)F(c,7)
        {
            if(i&1)ad(f[i][(a+cm[i][1])%7][(b+cm[i][3])%7][(c+cm[i][5])%7],f[i-1][a][b][c]);
            else ad(g[i][(a+cm[i][1])%7][(b+cm[i][3])%7][(c+cm[i][5])%7],g[i-1][a][b][c]);
        }
    }
}

int find(int x,int y,int wx,int wy,int k)
{
    int a=-1,b=-1,c=-1;
    if(0==x)a=wx;if(0==y)a=wy;
    if(1==x)b=wx;if(1==y)b=wy;
    if(2==x)c=wx;if(2==y)c=wy;
    int s=0;
    for(int i=0;i<7;i++)ad(s,g[k][a==-1?i:a][b==-1?i:b][c==-1?i:c]);
    return s;
}

void fd(int a,int b,int c,int w)
{
    ad(ans,mod-1ll*w*h[a][b]%mod);
    if(b!=c)
    ad(ans,mod-1ll*w*h[a][c]%mod);
}

void solve(int k,int*p,int w,int s)
{
    ad(ans,1ll*w*s%mod);
    for(int i=0;i<7;i++)p[i]=(7-p[i])%7;
    int A[3]={p[1],p[6]},B[3]={p[3],p[4]},C[3]={p[5],p[2]};
    fd(0,p[1],p[6],w);
    F(j,2)if(j==0||B[1]!=B[0])
    {
        int val=h[1][B[j]];
        ad(val,mod-d[0][A[0]][B[j]]);
        if(A[1]!=A[0])ad(val,mod-d[0][A[1]][B[j]]);
        ad(ans,mod-1ll*w*val%mod);
    }
    for(int t=0;t<2;t++)if(t==0||C[1]!=C[0])
    {
        int val=h[2][C[t]];
        F(i,2)if(i==0||A[1]!=A[0])ad(val,mod-d[1][A[i]][C[t]]);
        F(j,2)if(j==0||B[1]!=B[0])ad(val,mod-d[2][B[j]][C[t]]);
        F(i,2)if(i==0||A[1]!=A[0])
        F(j,2)if(j==0||B[1]!=B[0])ad(val,g[k][A[i]][B[j]][C[t]]);
        ad(ans,mod-1ll*val*w%mod);
    }
}

void sol()
{
    cin>>s;n=strlen(s);
    for(int i=0;i<n/2;i++)swap(s[i],s[n-i-1]);
    int r[8];ans=0;
    for(int i=1;i<7;i++)r[i]=1;
    for(int i=n-1;i;i--)
    {
        if(s[i]=='1')
        {
            int s=0;
            memset(h,0,sizeof h);
            memset(d,0,sizeof d);
            F(a,7)F(b,7)F(c,7)ad(s,g[i-1][a][b][c]),ad(h[0][a],g[i-1][a][b][c]),ad(h[1][b],g[i-1][a][b][c]),ad(h[2][c],g[i-1][a][b][c]);
            F(a,7)F(b,7)
            {
                d[0][a][b]=find(0,1,a,b,i-1);
                d[1][a][b]=find(0,2,a,b,i-1);
                d[2][a][b]=find(1,2,a,b,i-1);
            }
            F(a,7)F(b,7)F(c,7)
            {
                int p[8];
                for(int k=1;k<7;k++)p[k]=r[k];
                p[1]=(p[1]+a)%7;p[6]=(p[6]+7-a)%7;
                p[3]=(p[3]+b)%7;p[4]=(p[4]+7-b)%7;
                p[5]=(p[5]+c)%7;p[2]=(p[2]+7-c)%7;
                solve(i-1,p,f[i-1][a][b][c],s);
            }
            for(int j=1;j<7;j++)r[j]=(r[j]+cm[i][j])%7;
        }
    }
    int f=0;
    for(int i=1;i<7;i++)f|=(r[i]==0);
    if(!f&&s[0]=='1')ans++;
    cout<<ans%mod<<'\n';
    for(int i=0;i<n;i++)s[i]=0;
}

int main()
{
    init();
    cin>>T;
    while(T--)sol();
}