#include<bits/stdc++.h>
using namespace std;
const int N=1e3,S=1<<22;
const int mod=1e9+7;
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int n,m;
int d[N+5],e[N+5][N+5];
int p[N+5];
int f[S+5],g[N+5];
vector<int> ver;
int e1[50],e2[50];
bool cmp(int x,int y){ return d[x]<d[y]; }
int solve(){
    memset(e1,0,sizeof e1);
    memset(e2,0,sizeof e2);
    int n1=ver.size()/2,n2=(ver.size()+1)/2;
    int ans=0;
    for(int i=0;i<n1+n2;i++){
        e1[i]=0;
        for(int j=0;j<n1;j++)
            if(e[ver[i]][ver[j]])
                e1[i]|=1<<j;
    }
    for(int i=n1;i<n1+n2;i++){
        e2[i]=0;
        for(int j=n1;j<n1+n2;j++)
            if(e[ver[i]][ver[j]])
                e2[i]|=1<<j-n1;
    }
    f[0]=1;
    for(int i=1;i<(1<<n1);i++){
        int k=0; while(!(i>>k&1)) k++;
        f[i]=f[i^(1<<k)]+f[e1[k]&i];
    }
    g[0]=(1<<n1);
    for(int i=0;i<(1<<n2);i++){
        if(!g[i]) continue;
        ans=add(ans,f[g[i]-1]);
        for(int k=0;k<n2&&!(i>>k&1);k++)
            if((e2[k+n1]&i)==i)
                g[i|(1<<k)]=(g[i]-1&e1[k+n1])+1;
            else
                g[i|(1<<k)]=0;
    }
    return ans;
}
int main(){
    scanf("%d %d",&n,&m);
    for(int i=1,u,v;i<=m;i++){
        scanf("%d %d",&u,&v);
        e[u][v]=e[v][u]=1;
        d[u]++; d[v]++;
    }
    for(int i=1;i<=n;i++) p[i]=i;
    sort(p+1,p+1+n,cmp);
    int ans=0;
    for(int i=1;i<=n;i++){
        ver.clear();
        for(int j=i+1;j<=n;j++)
            if(e[p[i]][p[j]]) ver.push_back(p[j]);
        ans=add(ans,solve());
        // printf("%d\n",ans);
    }
    printf("%d\n",ans);
    return 0;
}