#!/bin/bash

         export PATH=/cursos/BI/bin/ncbiblast/bin:$PATH             # pel NCBI Blast
         cp /cursos/BI/bin/ncbiblast/.ncbirc ~/                     # pel NCBI Blast

 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 {
          blastdb=`grep $genome /cursos/BI/genomes/protists/genomes_list_and_info.tab | cut -f 4`
          blastall -p tblastn -i query.fa -d $blastdb -o tblastnB_SelO_E.siliculosus_AND_$genome.fa -F F -m 9 -e 0.00001
        } done