#!/usr/local/bin/bash # setup work-*.psq (ptn sequence data in FASTA format) for i in ND1 ND2 COX1 COX2 ATP8 ATP6 COX3 ND3 ND4L ND4 ND5 CYTB ; do echo $i j=work-$i fasta2ptn -t 2 $j.nsq > $j.psq done