working code, but runs slowly on cit-HepPh.seq