#!usr/bin/env python3
from collections import defaultdict
from collections import deque
from heapq import heappush, heappop
import sys
import math
import bisect
import random
def LI(): return list(map(int, sys.stdin.readline().split()))
def I(): return int(sys.stdin.readline())
def LS():return list(map(list, sys.stdin.readline().split()))
def S(): return list(sys.stdin.readline())[:-1]
def IR(n):
    l = [None for i in range(n)]
    for i in range(n):l[i] = I()
    return l
def LIR(n):
    l = [None for i in range(n)]
    for i in range(n):l[i] = LI()
    return l
def SR(n):
    l = [None for i in range(n)]
    for i in range(n):l[i] = S()
    return l
def LSR(n):
    l = [None for i in range(n)]
    for i in range(n):l[i] = LS()
    return l
sys.setrecursionlimit(1000000)
mod = 1000000007

#A
def A():
    x,n = LI()
    if x < n:
        print(0)
    else:
        print(10)
    return

#B
def B():
    n,x = LI()
    l = LI()
    d = 0
    ans = 1
    for i in range(n):
        d += l[i]
        if d <= x:
            ans += 1
    print(ans)
    return

#C
def C():
    w,h,x,y = LI()
    ans = w*h/2
    k = 0
    if x == w/2 and y == h/2:
        k = 1
    print(ans,k)
    return

#D
def D():
    n,k = LI()
    a = LI()
    b = [a[i] for i in range(n)]
    for i in range(n-1):
        b[i+1] += b[i]
    b.insert(0,0)
    ans = 0
    for l in range(n):
        r = bisect.bisect_left(b,b[l]+k)
        ans += n-r+1
    print(ans)
    return

#E
def E():
    n,m = LI()
    s = LI()
    t = LI()
    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 j in range(m+1):
        dp[0][j] = 1
    for i in range(1,n+1):
        for j in range(1,m+1):
            if s[i-1] == t[j-1]:
                dp[i][j] = dp[i-1][j]+dp[i][j-1]
            else:
                dp[i][j] = dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]
            dp[i][j] %= mod
    print(dp[n][m])
    return

#F
def F():
    n = I()
    x = [None for i in range(n)]
    y = [None for i in range(n)]
    f = [None for i in range(n)]
    for i in range(n):
        a,b,d = input().split()
        print(a,b,d)
        a = int(a)
        b = int(b)
        x[i] = a
        y[i] = b
        if d == "R":
            f[i] = (1,0)
        elif d == "L":
            f[i] = (-1,0)
        elif d == "U":
            f[i] = (0,1)
        else:
            f[i] = (0,-1)
    dx = max(x)-min(x)
    dy = max(x)-min(y)
    ans = dx*dy
    if ans == 0:
        print(0)
        quit()
    ma = 100000000
    p = [0,ma/3,2*ma/3,ma]
    for i in range(100):
        x1 = [x[i]+f[i][0]*p[1] for i in range(n)]
        x2 = [x[i]+f[i][0]*p[2] for i in range(n)]
        dx1 = max(x1)-min(x1)
        dx2 = max(x2)-min(x2)
        if dx1 < dx2:
            p[3] = p[2]
            p[1] = (p[3]-p[0])/3+p[0]
            p[2] = (p[3]-p[0])/3*2+p[0]
        else:
            p[0] = p[1]
            p[1] = (p[3]-p[0])/3+p[0]
            p[2] = (p[3]-p[0])/3*2+p[0]
        print(x1,p)
    a = sum(p)/4
    print(p)
    p = [-ma,2*ma/3-ma,2*ma/3*2-ma,ma]
    for i in range(100):
        x1 = [y[i]+f[i][1]*p[1] for i in range(n)]
        x2 = [y[i]+f[i][1]*p[2] for i in range(n)]
        dx1 = max(x1)-min(x1)
        dx2 = max(x2)-min(x2)
        if dx1 < dx2:
            p[3] = p[2]
            p[1] = (p[3]-p[0])/3+p[0]
            p[2] = (p[3]-p[0])/3*2+p[0]
        else:
            p[0] = p[1]
            p[1] = (p[3]-p[0])/3+p[0]
            p[2] = (p[3]-p[0])/3*2+p[0]
    b = sum(p)/4
    print(a,b)
    return

#Solve
if __name__ == "__main__":
    E()


Measured memory_bytes: 75996