#!/bin/bash

export PATH=/cursos/BI/bin/ncbiblast/bin:$PATH

cp /cursos/BI/bin/ncbiblast/.ncbirc ~/

export PATH=/cursos/BI/bin/exonerate/i386/bin:$PATH

export PATH=/cursos/BI/bin:$PATH

export WISECONFIGDIR=/cursos/BI/soft/wise-2.2.0/wisecfg

# Creamos los siguientes directorios para almacenar los resultados de este programa.

mkdir tcoffee_genewise

for genome in Crithidia_fasciculata Trypanosoma_congolense Leishmania_tarentolae Leishmania_donovani_BPK282A1 Ichthyophthirius_multifiliis_strain_G5 Phytophthora_capsici Sphaeroforma_arctica Fragilariopsis_cylindrus Dictyostelium_fasciculatum Dictyostelium_discoideum_AX4 Albugo_laibachii_Nc14; do {

for seq in selprotT_human.fa selprotT_mouse.fa selR1_human.fa selR2_human.fa selR3_human.fa selR_cerevisiae.fa selR_droso.fa selR_elegans.fa; do {

# Ejecutamos el programa t_coffee para realizar un alineamiento global de las querys con la proteína predicha a partir de nuestra anotación.

t_coffee $seq seleccion_genewise/genewise_${seq}.CONTRA.${genome}.fa -outfile tcoffee_genewise/tcoffee_genewise_${seq}.CONTRA.${genome}.txt

} done

} done