// 用givest 定源点汇点
// addedge一次加了正反两条边
// init 慎用
// S必须是0
// 输出方案注意是head开头
#include <bits/stdc++.h>
using namespace std;
const int N=100010,M=1000000; 
int n,m,cnt;
string a[110];
bool in(int x,int y){
    if(x>=0&&x<n&&y>=0&&y<m)return true;
    return false;
}
void no(){
    cout<<"NO\n";exit(0);
}
void yes(){
    cout<<"YES\n";
    for(int i=0;i<2*n-1;i++){
        cout<<a[i]<<"\n";
    }
    exit(0);
}
int fx[4]={1,0,-1,0};
int fy[4]={0,1,0,-1};
bool Less(char c){
    if(c=='2'||c=='4'||c=='7')return true;
    return false;
}
class Maxflow{
private:
	int nedge=1,p[2*M],nex[2*M],head[N],c[2*M],cur[2*M];
    int dist[2*N];
    int S,T;
	void Addedge(int a,int b,int v){
        // cout<<a<<" "<<b<<" "<<v<<"\n";
		p[++nedge]=b;nex[nedge]=head[a];head[a]=nedge;
		c[nedge]=v;
	}
    bool bfs(){
		queue<int>q;
		for(int i=S;i<=T;i++)dist[i]=-1;
		dist[S]=0;q.push(S);
		while(!q.empty()){
			int now=q.front();q.pop();
			for(int k=head[now];k;k=nex[k])if(dist[p[k]]==-1&&c[k]>0){
				dist[p[k]]=dist[now]+1;
				q.push(p[k]);
			}
		}
		return dist[T]>-1;
	}
	int dfs(int x,int low){
		if(x==T)return low;
		if(low==0)return 0;
		int used=0;
		for(int &k=cur[x];k;k=nex[k])if(dist[p[k]]==dist[x]+1&&c[k]>0){
			int a=dfs(p[k],min(c[k],low-used));
			c[k]-=a;c[k^1]+=a;used+=a;
			if(low==used)break;
		}
		if(used==0)dist[x]=-1;
		return used;
	}
public:
    void givest(int s,int t){
        S=s;T=t;
    }
	void init(){
		memset(head,0,sizeof head);
		nedge=1;
	}
    void addedge(int a,int b,int v){
        Addedge(a,b,v);
        Addedge(b,a,0);
	}
	int dinic(){
		int flow=0;
		while(bfs()){
			for(int i=S;i<=T;i++)cur[i]=head[i];
			flow+=dfs(S,1e9);
		}
        // cout<<"flow="<<flow<<"\n";
        if(flow*2!=cnt)no();
        for(int i=1;i<=n*m;i++){
            int A=(i-1)/m,B=(i-1)%m;
            if((A+B)%2==0)continue;
            for(int k=head[i];k;k=nex[k])if(c[k]==0&&p[k]!=S&&p[k]!=T){
                int xx=(i-1)/m,yy=(i-1)%m;
                int x=(p[k]-1)/m,y=(p[k]-1)%m;
                // cout<<i<<" "<<p[k]<<"\n";
                a[xx*2+x-xx][yy*2+y-yy]='.';
            }
        }
		return flow;
	}
};
Maxflow mf;
int main()
{
    cin>>n>>m;
    for(int i=0;i<2*n-1;i++){
        cin>>a[i];
        for(int j=0;j<2*m-1;j++){
            if(a[i][j]=='.')a[i][j]='#';
        }
    }
    mf.givest(0,n*m+1);
    int s=0,t=n*m+1;
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++)if(Less(a[i*2][j*2])){
            cnt++;
            if((i+j)%2){
                mf.addedge(s,i*m+j+1,1);
                for(int z=0;z<4;z++){
                    int x=i+fx[z],y=j+fy[z];
                    if(in(x,y)&&Less(a[x*2][y*2])){
                        if(a[x*2][y*2]+a[i*2][j*2]=='7'+'4')continue;
                        if(a[x*2][y*2]+a[i*2][j*2]=='4'+'4'){
                            if(!(x==0&&i==0||x==n-1&&i==n-1||y==0&&j==0||y==m-1&&j==m-1))continue;
                        }
                        mf.addedge(i*m+j+1,x*m+y+1,1);
                    }
                }
            }
            else mf.addedge(i*m+j+1,t,1);
        }
    }
    if(cnt%2)no();
    mf.dinic();
    yes();
}