#include<iostream>
#include<map>
#include<vector>
#include<algorithm>
#include<cstring>
using namespace std;
const int mod=1e9+7;
#define int long long
int s[1<<22];
int id[3000];
vector<int>v[3000];
int e[3000][3000];
int solve(vector<int>a){
    memset(id,-1,sizeof(id));
    for(int i=0;i<a.size();i++) id[a[i]]=i;
    int b=a.size()/2;
    int num=(a.size()-b);
    memset(s,0,sizeof(int)*(1<<b));
    for(int u=0;u<(1<<num);u++){
        bool fl=0;
        for(int j=0;j<num;j++){
            if(!((u>>j)&1)) continue;
            for(int k=0;k<num;k++){
                if(!((u>>k)&1)) continue;
                if(e[a[b+j]][a[b+k]]==0){fl=1;break;}
            }
            if(fl==1) break;
        }
        if(fl) continue;
        int mask=(1<<b)-1;
        for(int j=0;j<num;j++){
            if(!((u>>j)&1)) continue;
            int z=0;
            for(int k=0;k<b;k++){
                if(e[a[k]][a[b+j]]) z|=(1<<k);
            }
            mask&=z;
        }
        s[mask]++;
        // cerr<<mask<<endl;
    }
    for(int mid=1;mid<(1<<b);mid<<=1){
        for(int j=0;j<(1<<b);j+=(mid<<1)){
            for(int k=0;k<mid;k++){
                s[j+k]+=s[j+k+mid];
            }
        }
    }
    long long ans=0;
    for(int u=0;u<(1<<b);u++){
        bool fl=0;
        for(int i=0;i<b;i++){
            if(!((u>>i)&1)) continue;
            for(int j=0;j<b;j++){
                if(!((u>>j)&1)) continue;
                if(e[a[i]][a[j]]==0){fl=1;break;}
            }
            if(fl) break;
        }
        if(fl) continue;
        ans+=s[u];ans%=mod;
    }

    return ans%mod;
    // for(int i=0;i<b;i++){
        
    // }
    // for(int i=b;i<a.size();i++){
        
    // }
}
int n,m;
struct e{
    int x,y;
}d[10000];int deg[3020];
int sol(vector<int>a){
    
    sort(a.begin(),a.end(),[&](int u,int v){
        return deg[u]>deg[v];
    });
    int z=a.back();
    vector<int>tmp;
    for(int i=1;i<=n;i++) if(e[i][z]&&i!=z) deg[i]--,tmp.push_back(i);
    a.pop_back();
    int cur=solve(tmp);
    return cur+solve(a);

}
int32_t main(){
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int x,y;cin>>x>>y;
        e[x][y]=e[y][x]=1;
        // d[i]=e{x,y};
        deg[x]++,deg[y]++;
    }
    for(int i=1;i<=n;i++) e[i][i]=1;
    vector<int>p;
    for(int i=1;i<=n;i++) p.push_back(i);
    cout<<(sol(p)-1+mod)%mod;
}