import sys
input=sys.stdin.readline
mod=1e9+7

if __name__ == '__main__':
    N,M=list(map(int,input().split()))
    S=list(map(int,input().split()))
    T=list(map(int,input().split()))
    dp=[[0 for i in range(M+1)] for j in range(N+1)]
    for i in range(N+1):
        dp[i][0]=1
    for i in range(M+1):
        dp[0][i]=1


    for i in range(N):
        for j in range(M):
            dp[i + 1][j + 1] = dp[i + 1][j]+dp[i ][j + 1]-dp[i][j]
            if S[i]==T[j]:
                dp[i+1][j+1]+=dp[i][j]

    print(int(dp[N][M]%mod))
