#!/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