#!/usr/bin/perl -w
#time perl ../../run-blast.pl ../../di1-di2-tri-arch.fasta 0.01

open(IN,"/db2/data/kcsabadata/dUTPase/AllGenomeFileNames-wo-plasmids.txt");
while (<IN>) {
        chomp;
        $file=$_;
        if ( $file =~ /(N._\d+\.fna)/ ) {
            my $actual = $1;
            system("/db2/data/kcsabadata/ncbi-blast-2.2.28+/bin/tblastn -query $ARGV[0] -evalue $ARGV[1] -db /db2/data/kcsabadata/dUTPase/blast-allgenomes/db/$file > $actual");
        }
}
close IN;

system("ls *fna > list");
