This function measures the similarity of two PWM matrix in three measurements:
"normalised Euclidean distance", "Pearson correlation" and "Kullback Leibler divergence".
When pwmSubject and pwmQuery have different number of columns,
the smaller PWM will be shifted from the start position of larger PWM
and compare all the possible alignments.
Only the smallest distance, divergence or largest correlation will be reported.
References
Linhart, C., Halperin, Y., & Shamir, R. (2008). Transcription factor and microRNA motif discovery: The Amadeus platform and a compendium of metazoan target sets. Genome Research, 18(7), 1180-1189. doi:10.1101/gr.076117.108