#include <bits/stdc++.h>
using namespace std;
const long long MOD = 1000000007;
string andstring(const string& a, const string& b) {
  string res(a);
  for (int i = 0; i < b.size(); i++)
    res[i] = ((a[i] - '0') & (b[i] - '0')) + '0';
  return res;
}
string invstring(const string& a) {
  string res(a);
  for (int i = 0; i < a.size(); i++) res[i] = '0' + '1' - res[i];
  return res;
}
const int MAXM = 1001;
long long dp[MAXM];
long long cnk[MAXM][MAXM];
int main() {
  int m, n;
  cin >> m >> n;
  vector<string> s(n);
  for (int i = 0; i < n; i++) cin >> s[i];
  for (int i = 0; i < n; i++) s.push_back(invstring(s[i]));
  sort(s.begin(), s.end());
  s.erase(unique(s.begin(), s.end()), s.end());
  n = s.size();
  vector<string> minstring(m);
  for (int i = 0; i < m; i++) {
    minstring[i] = string(m, '1');
    for (int j = 0; j < n; j++)
      if (s[j][i] == '1') minstring[i] = andstring(minstring[i], s[j]);
  }
  sort(minstring.begin(), minstring.end());
  minstring.erase(unique(minstring.begin(), minstring.end()), minstring.end());
  for (int i = 0; i <= m; i++) {
    cnk[i][0] = 1;
    for (int j = 1; j < i; j++)
      cnk[i][j] = (cnk[i - 1][j - 1] + cnk[i - 1][j]) % MOD;
    cnk[i][i] = 1;
  }
  dp[0] = 1;
  for (int i = 1; i <= m; i++) {
    dp[i] = 0;
    for (int j = 0; j < i; j++)
      dp[i] = (dp[i] + dp[i - j - 1] * cnk[i - 1][j] % MOD) % MOD;
  }
  long long res = 1;
  for (int i = 0; i < minstring.size(); i++) {
    int sum = 0;
    for (int j = 0; j < m; j++) sum += minstring[i][j] == '1' ? 1 : 0;
    res = res * dp[sum] % MOD;
  }
  printf("%lld\n", res);
  return 0;
}
