#!/bin/bash

interruptor="0"

for query in ./query/*; do {

for genome in  A.taiwanensis B.hominis C.merolae C.muris C.owczarzaki C.parvum E.siliculosus H.arabidopsidis N.gruberi P.pallidum P.tricornutum P.ultimum S.parasitica T.trahens; do {
	if [ "$interruptor" = "0" ]; then
		mkdir ./blast_extended/$genome
	fi
          blastdb=`grep $genome /cursos/BI/genomes/protists/genomes_list_and_info.tab | cut -f 4`
          blastall -p tblastn -d $blastdb -i ./$query.txt -o ./blast_extended/$genome/${genome}_$query.txt 
	
	  echo ${query}_$genome fet

        } done
interruptor="1"
} done