#!/usr/bin/perl -w

open(IN, "< genomes.list") or die "cannot open genomes.list: $!";
$i=0;
@genome=();
while(defined($line=<IN>)){
    chomp($line);
    $genome[$i]=$line;
    $i++;

}
close(IN);

$n=@genome; #number of genomes to work on
$db_folder="/Users/jpgogarten/genomes/";


for($i=0;$i<$n;$i++){

    system("pepstats $db_folder$genome[$i].faa -outfile $genome[$i].pepstats");

    
}