#include<bits/stdc++.h>
#define MN 55
const int gx[] = {0, 0, 1, -1};
const int gy[] = {1, -1, 0, 0};
int n, m, N, M;
int a[MN][MN];
char s[MN << 1][MN << 1], ans[MN << 1][MN << 1];
std::vector<int> go[MN * MN];
int ha(int x, int y) {
    return (x - 1) * m + y;
}
bool check(int x, int y) {
    return 1 <= x && x <= n && 1 <= y && y <= m;
}
int mat[MN * MN];
bool vis[MN * MN];
bool find(int x) {
    for(int y : go[x]) {
        if(!vis[y]) {
            vis[y] = 1;
            if(mat[y] == 0 || find(mat[y])) {
                mat[y] = x;
                return 1;
            }
        }
    }
    return 0;
}
void solve() {
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= 2 * n - 1; i++) scanf("%s", s[i] + 1);
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            int x = (i << 1) - 1, y = (j << 1) - 1;
            int c = s[x][y] - '0';
            int o = 0;
            if(i == 1 || i == n) o++;
            if(j == 1 || j == m) o++;
            if(o == 2) {
                if(c != 2 && c != 3) {
                    puts("NO");
                    return ;
                }
                a[i][j] = c - 2;
            } else if(o == 1) {
                if(c != 4 && c != 5) {
                    puts("NO");
                    return ;
                }
                a[i][j] = c - 4;
            } else {
                if(c != 7 && c != 8) {
                    puts("NO");
                    return ;
                }
                a[i][j] = c - 7;
            }
        }
    }
    // std::cerr << "xxx\n";
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            if(!a[i][j]) {
                if((i ^ j) & 1) ++N;
                else ++M;
                // std::cerr << N << " " <<  M << "\n";
            }
        }
    }
    if(N != M) {
        puts("NO");
        return ;
    }
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            if(!a[i][j]) {
                if((i ^ j) & 1) ++N;
                else ++M;
                for(int u = 0; u < 4; u++) {
                    int ni = i + gx[u], nj = j + gy[u];
                    if(i == 1 && u == 2 && 1 < j && j < m) continue;
                    if(i == n && u == 3&& 1 < j && j < m) continue;
                    if(j == 1 && u == 0 && 1 < i && i < n) continue;
                    if(j == m && u == 1 && 1 < i && i < n) continue;

                    if(check(ni, nj) && a[ni][nj] == 0) go[ha(i, j)].push_back(ha(ni, nj));
                // std::cerr << i << " " << j << " " << a[i][j] << "\n";;
                }
            }
        }
    }
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            memset(vis, 0, sizeof vis);
            if(!a[i][j] && ((i ^ j) & 1)) {
                if(!find(ha(i, j))) {
                    puts("NO");
                    return ;
                }
            }
        }
    }
    for(int i = 1; i < 2 * n; i++) {
        for(int j = 1; j < 2 * m; j++) {
            ans[i][j] = s[i][j];
            if((i & 1 ^ 1) || (j & 1 ^ 1)) ans[i][j] = '#';
        }
    }
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            if(!a[i][j] && ((i ^ j) & 1 ^ 1)) {
                int u = mat[ha(i, j)], v = ha(i, j);
                if(u == v - m) {
                    ans[i * 2 - 2][j * 2 - 1] = '.';
                } else if(u == v + m) {
                    ans[i * 2][j * 2 - 1] = '.';
                } else if(u == v - 1) {
                    ans[i * 2 - 1][j * 2 - 2] = '.';
                } else 
                    ans[i * 2 - 1][j * 2] = '.';
            }
        }
    }
    puts("YES");
    for(int i = 1; i < 2 * n; i++) printf("%s\n", ans[i] + 1);
}
int main() {
    solve();
}