powered by
cat for phase information
cat_phase( input.seq, input.ph, all.ph, ct, seq.num, twopt.phase.number, hmm.phase.number )