calculate I'
calI(freq1, freq2, P)
freq1
position frequence for matrix 1 position j
freq2
position frequence for matrix 2 position j
P
background of profile1
numeric(1)