#include <bits/stdc++.h>
using namespace std;
const int INF = 1e9 + 7;
long long dp[1005][1005];
string S[55];
vector<int> G[1005];
int sz;
bool vis[1005];
void dfs(int v) {
  vis[v] = true;
  sz++;
  for (int k = 0; k < (G[v].size()); k++) {
    int u = G[v][k];
    if (!vis[u]) dfs(u);
  }
}
int main() {
  ios::sync_with_stdio(false);
  int n, m;
  cin >> n >> m;
  dp[0][0] = 1;
  for (int i = 1; i <= (n); i++)
    for (int j = 1; j <= (i); j++)
      dp[i][j] = (dp[i - 1][j] * j + dp[i - 1][j - 1]) % INF;
  for (int k = 0; k < (m); k++) cin >> S[k];
  for (int i = 0; i < (n); i++)
    for (int j = 0; j < (n); j++) {
      bool ok = true;
      for (int k = 0; k < (m); k++)
        if (S[k][i] != S[k][j]) ok = false;
      if (ok) G[i].push_back(j);
    }
  long long ans = 1LL;
  for (int k = 0; k < (n); k++)
    if (!vis[k]) {
      sz = 0;
      dfs(k);
      long long ca = 0LL;
      for (int k = 1; k <= (sz); k++) ca += dp[sz][k];
      ca %= INF;
      ans = ans * ca % INF;
    }
  cout << ans << '\n';
  return 0;
}
