#include<bits/stdc++.h>
using namespace std;
int dx[]={-1,0,1,0},dy[]={0,1,0,-1};
int n,m;
int X,Y,Z;
struct Array{
    int a[1000005];
    int* operator [](int i){return a+(i-1)*m;}
}ans;
vector<pair<int,int>> vec;
int flag,flag1,flag2;
void print(){
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)Z-=!ans[i][j];
    while(Z--)ans[vec.back().first][vec.back().second]=0,vec.pop_back();
    puts("Yes");
    if(flag1)for(int i=1;i<=n-i+1;i++)for(int j=1;j<=m;j++)swap(ans[i][j],ans[n-i+1][j]);
    if(flag2)for(int i=1;i<=n;i++)for(int j=1;j<=m-j+1;j++)swap(ans[i][j],ans[i][m-j+1]);
    if(flag){
        for(int i=1;i<=m;i++,puts(""))
            for(int j=1;j<=n;j++)putchar(ans[j][i]?'#':'.');
    }else{
        for(int i=1;i<=n;i++,puts(""))
            for(int j=1;j<=m;j++)putchar(ans[i][j]?'#':'.');
    }
}
void solve(){
    scanf("%d%d%d%d%d",&n,&m,&X,&Y,&Z);
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)ans[i][j]=0;
    flag=0;flag1=flag2=0;
    vec.clear();
    if(n==1){
        if(Y<=Z){
            for(int j=m;j>Y;j--)vec.emplace_back(1,j),ans[1][j]=1;
            print();
            return;
        }
        if(m-Y+1<=Z){
            for(int j=1;j<Y;j++)vec.emplace_back(1,j),ans[1][j]=1;
            print();
            return;
        }
    }
    if(n==2&&m-2+(Y==1||Y==m)<=Z){
        for(int i=1;i<=m;i++)vec.emplace_back(3-X,i),ans[3-X][i]=1;
        if(Y>1)vec.emplace_back(X,1),ans[X][1]=1;
        if(Y<m)vec.emplace_back(X,m),ans[X][m]=1;
        print();
        return;
    }
    if(n>=3&&m>=3&&((n-1)*(m-1)+1)/2-1+(X==1||X==n)+(Y==1||Y==m)<=Z){
        flag1=flag2=0;
        if(X>n-X+1)flag1=1,X=n-X+1;
        if(Y>m-Y+1)flag2=1,Y=m-Y+1;
        {
            vector<pair<int,int>> vv;
            for(int j=1;j<m;j++)vv.emplace_back(1,j),ans[1][j]=1;
            for(int i=1;i<n;i++)vv.emplace_back(i,m),ans[i][m]=1;
            for(int j=m;j>1;j--)vv.emplace_back(n,j),ans[n][j]=1;
            for(int i=n;i>1;i--)vv.emplace_back(i,1),ans[i][1]=1;
            int all=0;
            if(X==1||Y==1||X==n||Y==m)ans[X][Y]=0,all++;
            if(X==1&&Y==1)ans[1][2]=0,all++;
            int now=0;
            for(auto p:vv)
                if(!ans[p.first][p.second])now++;
                else if(now==all)vec.push_back(p);
            for(auto p:vv)
                if(!ans[p.first][p.second])now--;
                else if(now==all)vec.push_back(p);
        }
        if(n&1){
            int p=Y+(Y==1);
            for(int i=3;i<n;i+=2){
                for(int j=2;j<p;j++)vec.emplace_back(i,j),ans[i][j]=1;
                for(int j=m-1;j>p;j--)vec.emplace_back(i,j),ans[i][j]=1;
            }
        }else{
            int p=Y+(Y==1);
            for(int i=3;i<n-1;i+=2){
                for(int j=2;j<p;j++)vec.emplace_back(i,j),ans[i][j]=1;
                for(int j=m-1;j>p;j--)vec.emplace_back(i,j),ans[i][j]=1;
            }
            for(int j=3;j<m-1;j+=2)vec.emplace_back(n-1,j),ans[n-1][j]=1;
        }
        print();
        return;
    }
    swap(n,m);swap(X,Y),flag=1;
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)ans[i][j]=0;
    if(n==1){
        if(Y<=Z){
            for(int j=m;j>Y;j--)vec.emplace_back(1,j),ans[1][j]=1;
            print();
            return;
        }
        if(m-Y+1<=Z){
            for(int j=1;j<Y;j++)vec.emplace_back(1,j),ans[1][j]=1;
            print();
            return;
        }
    }
    if(n==2&&m-2+(Y==1||Y==m)<=Z){
        for(int i=1;i<=m;i++)vec.emplace_back(3-X,i),ans[3-X][i]=1;
        if(Y>1)vec.emplace_back(X,1),ans[X][1]=1;
        if(Y<m)vec.emplace_back(X,m),ans[X][m]=1;
        print();
        return;
    }
    puts("No");
}
int main(){
    // freopen("1.in","r",stdin);
    // freopen("1.out","w",stdout);
    int _;scanf("%d",&_);
    while(_--)solve();
    return 0;
}