#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int n,m,shu1,shu2,point[1000005],dp[1000005];
int ans,num[1000005];
int moc(int x)
{
	if(x>=mod) 
        return x-mod;
	if(x<0) 
        return x+mod;
	return x;
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&shu1,&shu2);
		shu1--,shu2--;
		point[shu1]|=1<<shu2;
	}
	int maxa=(1<<n)-1;
	dp[0]=1;
	for(int i=0;i<=maxa;i++)
	{
		int now=0;
		for(int j=0;j<n;j++)
		    if(i>>j&1) 
                now++;
		num[i]=1<<now;
	}
	for(int s=1;s<=maxa;s++)
	if((s&1)==(s>>1&1))
	{
		for(int ss=s;ss;ss=(ss-1)&s)
        {
            if((ss&1)==(ss>>1&1))
		    {
			    int t=s^ss,now=1;
			    for(int i=0;i<n;i++)
			    if(ss>>i&1)
			        now=1ll*now*num[point[i]&t]%mod;
			    else if(t>>i&1)now=1ll*now*(num[point[i]&ss]-1)%mod;
			        dp[s]=moc(dp[s]+1ll*now*dp[t]%mod);
		    }
        }	    
	}
	ans=1;
	for(int i=1;i<=m;i++)
	    ans=moc(ans<<1);
	ans=moc(ans-dp[maxa]);
	printf("%d",ans);
	return 0;
}