modify annotatePTC.R to fit my case