#include <bits/stdc++.h>
#define mod 1000000007
#define F first
#define S second
#define ll long long
#define N 100010
using namespace std;
int ksm(int x,int y){
	int ret=1;
	while(y){
		if(y&1) ret=(1ll*ret*x)%mod;
		x=(1ll*x*x)%mod;
		y>>=1;
	}
	return ret;
}
int n,m,inv[N],fact[N],f[N][2],ans;
vector<pair<int,int> > vt[N];
bool vis[N];
void dfs(int x,int lst){
	if(vt[x].size()==1&&lst!=-1){
		f[x][0]=1,f[x][1]=0;
		return;
	}
	int val=inv[vt[x].size()-1];
	f[x][0]=f[x][1]=0;
	for(int i=0;i<vt[x].size();i++){
		int y=vt[x][i].F,c=vt[x][i].S;
		if(y==lst) continue;
		dfs(y,x);
		int nw[2];
		nw[0]=nw[1]=0;
		nw[vis[c]]=(1ll*nw[vis[c]]+1ll*f[y][0]*inv[vt[y].size()-1])%mod;
		nw[1]=(1ll*nw[1]+1ll*f[y][1]*inv[vt[y].size()-1])%mod;
		int cur=(1ll*f[x][0]*nw[1]+1ll*f[x][1]*nw[0]+1ll*f[x][1]*nw[1])%mod;
		ans=(1ll*ans+1ll*cur*val)%mod;
		f[x][0]=(f[x][0]+nw[0])%mod,f[x][1]=(f[x][1]+nw[1])%mod;
	}
	return;
}
int main(){
	freopen("traverse.in","r",stdin);
	freopen("traverse.out","w",stdout);
	int _,T;
	inv[1]=fact[0]=fact[1]=1;
	for(int i=2;i<N;i++){
		inv[i]=(1ll*inv[mod%i]*(mod-mod/i))%mod;
		fact[i]=(1ll*fact[i-1]*i)%mod;
	}
	inv[0]=1;
	scanf("%d%d",&_,&T);
	while(T--){
		scanf("%d%d",&n,&m);
		for(int i=1;i<n;i++) vis[i]=false;
		for(int i=0;i<n;i++) vt[i].clear();
		for(int i=1;i<n;i++){
			int x,y;
			scanf("%d%d",&x,&y);
			x--,y--;
			vt[x].push_back(make_pair(y,i));
			vt[y].push_back(make_pair(x,i));
		}
		while(m--){
			int x;
			scanf("%d",&x);
			vis[x]=true;
		}
		if(n==2){
			puts("1");
			continue;
		}
		ans=0;
		for(int i=0;i<n;i++) if(vt[i].size()>1){
			dfs(i,-1);
			break;
		}
		for(int i=0;i<n;i++) ans=(1ll*ans*fact[vt[i].size()-1])%mod;
		printf("%d\n",ans);
	}
	return 0;
}