Extract regions from a sequence emboss extractseq -sequence $input1 -outseq $out_file1 -regions $regions -separate $separate -osformat2 $out_format1 -auto You can view the original documentation here_. .. _here: http://emboss.sourceforge.net/apps/release/5.0/emboss/apps/extractseq.html