filename=${1} # without extension
input="${filename}.paml"
output="${filename}.fasta"
grep ^SEQ ${input} | sort -d | sed 's/ \{8\}/:/' | tr -d " " | awk -F":" '$1==last {printf ".%s",$2; next} NR>1 {print "";} {last=$1; printf "%s",$0} END{print "";}' | tr -d "." | sed -e 's/^/>/' | tr ":" "\n" > ${output}