pwmtool_dir='scripts/pwm_tools' data_dir="data/pwm/jaspar_2018_clustering" for file_tf in $(ls $data_dir/*tf) do file_prob=$(basename $file_tf | sed s/.tf/_prob.mat/) python3.6 $pwmtool_dir/extract_transfac_pwm.py -i $file_tf --norm > $data_dir/$file_prob done