#include<bits/stdc++.h>
#define N 3005
using namespace std;
int n,m;
string s[105];
int mch[55][55];
int id(int x,int y){
    return (x-1)*m+y;
}
struct edge{
    int to,nxt,w;
}e[N<<3];
int head[N],cnt=1;
void adde(int a,int b,int c){
    // cerr<<a<<" "<<b<<" "<<c<<endl;
    e[++cnt]=(edge){b,head[a],c},head[a]=cnt;
    e[++cnt]=(edge){a,head[b],0},head[b]=cnt;
}
int ss,t;
int cur[N],dep[N],gap[N];
int mxfl=0;
queue<int> q;
void bfs(){
    for(int i=0;i<=t+5;++i) dep[i]=-1,gap[i]=0;
    dep[t]=0,gap[0]=1;
    q.push(t);
    while(!q.empty()){
        int u=q.front(); q.pop();
        for(int i=head[u];i;i=e[i].nxt){
            int v=e[i].to; if(dep[v]!=-1) continue;
            dep[v]=dep[u]+1; ++gap[dep[v]];
            q.push(v);
        }
    }
}
int dfs(int x,int fl){
    if(x==t){
        mxfl+=fl;
        return fl;
    }
    int usd=0;
    for(int i=cur[x];i;cur[x]=i=e[i].nxt){
        int v=e[i].to;
        if(e[i].w==0||dep[x]!=dep[v]+1) continue;
        int tmp=dfs(v,min(fl-usd,e[i].w));
        e[i].w-=tmp,e[i^1].w+=tmp;
        usd+=tmp;
        if(usd==fl) return fl;
    }
    --gap[dep[x]++]==0?dep[ss]=t+1:++gap[dep[x]];
    return usd;
}
void ISAP(){
    bfs();
    while(dep[ss]<=t){
        for(int i=1;i<=t;++i) cur[i]=head[i];
        dfs(ss,100000);
    }
}
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin>>n>>m;
    for(int i=1;i<=n*2-1;++i){
        cin>>s[i];
        for(int j=1;j<=m*2-1;++j)
            if(s[i][j-1]=='.') s[i][j-1]='#';
    }
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            int val=8;
            int cnt=0;
            if(i==1||i==n) ++cnt;
            if(j==1||j==m) ++cnt;
            if(cnt==1) val=5;
            if(cnt==2) val=3;
            if(val!=s[i*2-1][j*2-2]-'0')
                mch[i][j]=1;
        }
    }
    ss=n*m+1;
    t=n*m+2;
    int tot=0;
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j) if(mch[i][j]){
            ++tot;
            if((i+j)&1)
                adde(ss,id(i,j),1);
            else
                adde(id(i,j),t,1);
            if((j==1||j==m||(i>2&&i<n))&&mch[i-1][j]){
                if((i+j)&1)
                    adde(id(i,j),id(i-1,j),1);
                else
                    adde(id(i-1,j),id(i,j),1);
            }
            if((i==1||i==n||(j>2&&j<m))&&mch[i][j-1]){
                if((i+j)&1)
                    adde(id(i,j),id(i,j-1),1);
                else
                    adde(id(i,j-1),id(i,j),1);
            }
        }
    }
    ISAP();
    // cerr<<mxfl<<endl;
    if(mxfl*2!=tot){
        cout<<"NO\n";
        return 0;
    }
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j) if(mch[i][j]&&((i+j)&1)){
            int id1=id(i,j);
            int id2=id(i-1,j);
            int id3=id(i,j-1);
            int id4=id(i+1,j);
            int id5=id(i,j+1);
            for(int k=head[id1];k;k=e[k].nxt){
                if(e[k].to==id2&&e[k].w==0){
                    s[i*2-2][j*2-2]='.';
                    break;
                }
                if(e[k].to==id3&&e[k].w==0){
                    s[i*2-1][j*2-3]='.';
                    break;
                }
                if(e[k].to==id4&&e[k].w==0){
                    s[i*2][j*2-2]='.';
                    break;
                }
                if(e[k].to==id5&&e[k].w==0){
                    s[i*2-1][j*2-1]='.';
                    break;
                }
            }
        }
    cout<<"YES\n";
    for(int i=1;i<=n*2-1;++i)
        cout<<s[i]<<'\n';
    return 0;
}