python3 version with pvalues for peaks in workflow.txt report