#include<bits/stdc++.h>
using namespace std;
const int QSZ=10000;
int q[11111],ql,qr;
int main()
{
	ios_base::sync_with_stdio(false);
	int T;
	cin>>T;
	while(T--)
	{
		int n,m;
		cin>>n>>m;
		vector<int> L(m+5),R(m+5);
		map<int,int> mp;
		for(int i=1;i<=m;i++)
		{
			cin>>L[i]>>R[i];
			mp[L[i]]=mp[R[i]+1]=0;
		}
		int idx=0;
		for(auto &[x,y]:mp)
		{
			y=++idx;
		}
		for(int i=1;i<=m;i++)
		{
			L[i]=mp[L[i]];
			R[i]=mp[R[i]+1];
		}
		vector<vector<pair<int,int>>> G(idx+5);
		for(int i=1;i<=m;i++)
		{
//				cerr<<L[i]<<' '<<R[i]<<endl;
			G[L[i]].emplace_back(R[i],1);
//				cerr<<R[i]<<' '<<L[i]<<endl;
			G[R[i]].emplace_back(L[i],-1);
		}
		for(int i=1;i<idx;i++)
		{
//				cerr<<i+1<<' '<<i<<endl;
			G[i+1].emplace_back(i,0);
		}
		//[li,ri)
		auto check=[&](int P,int Q)
		{
			vector<long long> dis(idx+5),cnt(idx+5),inq(idx+5);
			ql=0;qr=0;
			for(int i=1;i<=idx;i++)
			{
				q[qr]=i;
				qr=(qr+1)%QSZ;
				inq[i]=1;
			}
			while(ql!=qr)
			{
				int u=q[ql];
				ql=(ql+1)%QSZ;
				inq[u]=0;
				for(auto [v,w]:G[u])
				{
					if(w==1)w=P;else if(w==-1)w=-Q;
					if(dis[v]>dis[u]+w)
					{
						dis[v]=dis[u]+w;
						if(not inq[v])
						{
							cnt[v]++;
							if(cnt[v]>idx)return false;
							inq[v]=1;
							q[qr]=v;
							qr=(qr+1)%QSZ;
						}
					}
				}
			}
			return true;
		};
		/*
		for(int i=1;i<=5;i++)
			for(int j=1;j<=i;j++)
			{
				if(gcd(i,j)!=1)continue;
				cerr<<i<<' '<<j<<' '<<check(i,j)<<endl;
			}
		*/
		const int lim=4000*4000*1.1;
		struct frac
		{
			int p,q;
			frac operator+(const frac &f)const{return {p+f.p,q+f.q};}
			frac operator*(const int k)const{return {p*k,q*k};}
			bool operator<(const frac &f)const{return 1ll*p*f.q<1ll*q*f.p;}
		}ans{lim,1},l{1,1},r{1,0};
		if(check(1,1))
		{
			ans={1,1};
		}
		int dir=0,fir=1;
		while(1)
		{
//			cerr<<"gp "<<l.p<<"/"<<l.q<<'~'<<r.p<<"/"<<r.q<<' '<<dir<<endl;
			if((l+r).p>lim or (l+r).q>lim)break;
			if(dir==0)
			{
				int kl=(!fir),kr=min((lim-r.p)/max(l.p,1),(lim-r.q)/max(l.q,1));
			//	fir=0;
				while(kl<kr)
				{
					int km=(kl+kr+1)/2;
					frac tmp=l*km+r;
					if(check(tmp.p,tmp.q))
					{
						if(tmp<ans)ans=tmp;
						kl=km;
					}
					else
						kr=km-1;
				}
				r=l*kl+r;
				if(r<ans and check(r.p,r.q))ans=r;
				dir^=1;
			}
			else
			{
				int kl=(!fir),kr=min((lim-l.p)/max(r.p,1),(lim-l.q)/max(r.q,1));
				while(kl<kr)
				{
					int km=(kl+kr+1)/2;
					frac tmp=l+r*km;
					if(check(tmp.p,tmp.q))
					{
						if(tmp<ans)ans=tmp;
						kr=km-1;
					}
					else
						kl=km;
				}
				l=l+r*kl;
				if(l<ans and check(l.p,l.q))ans=l;
				dir^=1;
			}
		}
		const int MOD=998244353;
		auto poww=[&](long long x,int y)
		{
			long long ret=1;
			while(y)
			{
				if(y&1)ret*=x,ret%=MOD;
				x*=x,x%=MOD;
				y>>=1;
			}
			return ret;
		};
//		cerr<<"ans = "<<ans.p<<'/'<<ans.q<<endl;
		cout<<1ll*ans.p*poww(ans.q,MOD-2)%MOD<<"\n";
	}
	
	return 0;
}