Align nucleic coding regions given the aligned proteins emboss tranalign -asequence $input1 -bsequence $input2 -outseq $out_file1 -table $table -osformat3 $out_format1 -auto You can view the original documentation here_. .. _here: http://emboss.sourceforge.net/apps/release/5.0/emboss/apps/tranalign.html