#include<bits/stdc++.h>
#define ll   long long
#define pb   emplace_back
#define mp   make_pair
#define orz  1000000007
#define fi   first
#define se   second
using namespace std;
struct E{int x,y,z;};
struct F{
    int n,s,t,i[21111],l[21111],q[21111],ql,qr;
    vector<E> v[21111];
    void add(int x,int y,int z){
        E X,Y;
        X.x=y,X.y=z,X.z=v[y].size();
        Y.x=x,Y.y=0,Y.z=v[x].size();
        v[x].pb(X),v[y].pb(Y);
    }
    void bfs(){
        memset(l,-1,sizeof(l));
        ql=0,qr=1,q[1]=s,l[s]=0;
        while(ql<qr){
            int x=q[++ql];
            for(int i=0;i<v[x].size();++i)if(v[x][i].y&&l[v[x][i].x]==-1)q[++qr]=v[x][i].x,l[v[x][i].x]=l[x]+1;
        }
    }
    int dfs(int x,int y){
        if(x==t) return y;
        int ret=0;
        for(;i[x]<v[x].size();++i[x]){
            E &e=v[x][i[x]];
            if(e.y&&l[e.x]==l[x]+1){
                int d=dfs(e.x,min(e.y,y));
                if(!d) continue;
                e.y-=d,y-=d,ret+=d,v[e.x][e.z].y+=d;
                if(!y) return ret;
            }
        }
        l[x]=-1;
        return ret;
    }
    int dinic(){
        int ans=0;
        while(1){
            bfs();
            if(l[t]==-1) return ans;
            memset(i,0,sizeof(i));
            ans+=dfs(s,orz);
        }
    }
}_;        
int n,m,num[205][205],cnt,_n;
vector<pair<int,int>> v[40005];
char c[205][205];
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n*2-1;++i){
		for(int j=1;j<=m*2-1;++j)scanf(" %c",c[i]+j);
	}
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j)if((i+j)&1){
			int x=i*2-1,y=j*2-1;
			if(c[x][y]=='2'||c[x][y]=='4'||c[x][y]=='7')num[x][y]=++cnt;
		}
	}
	_n=cnt;
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j)if((i+j+1)&1){
			int x=i*2-1,y=j*2-1;
			if(c[x][y]=='2'||c[x][y]=='4'||c[x][y]=='7')num[x][y]=++cnt;
		}
	}
	if(cnt!=_n*2){
		printf("NO\n");
		return 0;
	}
	_.s=cnt+1;
	_.n=_.t=_.s+1;
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j)if((i+j)&1){
			int x=i*2-1,y=j*2-1;
			if(!num[x][y]) continue;
			if(x>1&&!(c[x][y]=='4'&&i==n)&&num[x-2][y])_.add(num[x][y],num[x-2][y],1),v[num[x][y]].pb(mp(x-1,y));
			if(y>1&&!(c[x][y]=='4'&&j==m)&&num[x][y-2])_.add(num[x][y],num[x][y-2],1),v[num[x][y]].pb(mp(x,y-1));
			if(i<n&&!(c[x][y]=='4'&&x==1)&&num[x+2][y])_.add(num[x][y],num[x+2][y],1),v[num[x][y]].pb(mp(x+1,y));
			if(j<m&&!(c[x][y]=='4'&&y==1)&&num[x][y+2])_.add(num[x][y],num[x][y+2],1),v[num[x][y]].pb(mp(x,y+1));
		}
	}
    for(int i=1;i<=_n;++i)_.add(_.s,i,1),_.add(_n+i,_.t,1);
    int ans=_.dinic();
    if(ans!=_n){
    	printf("NO\n");
    	return 0;
    }
    printf("YES\n");
    for(int i=1;i<=n*2-1;++i){
    	for(int j=1;j<=m*2-1;++j)if(c[i][j]=='.')c[i][j]='#';
    }
    for(int i=1;i<=_n;++i){
        for(int j=0;j<_.v[i].size();++j){
            if(!_.v[i][j].y&&_.v[i][j].x!=_.s){
                c[v[i][j].fi][v[i][j].se]='.';
                break;
            }
        }
    }
    for(int i=1;i<=n*2-1;++i){
    	for(int j=1;j<=m*2-1;++j)printf("%c",c[i][j]);
    	printf("\n");
    }
	return 0;
}