#!/usr/bin/perl -w
use strict;

my $p; #protein
my $seq; #fasta sequence

#First the program asks for the name of the protein we want to work with, in order to name the files and to know which query to take

print "Hi! What protein are you working with? ";
$p = <STDIN>;
chomp ($p);

print “Running blast…\n”;

system ("tblastn -query $p.anole.fa -db /cursos/20428/BI/genomes/2016/Gekko_japonicus/genome.fa -outfmt 7 -out $p.blast");

print "Blast ended!:)\n”;

#The program shows us the result from blast in the screen so we can choose which hits to analyze

system ("more $p.blast");

print "\n";

print "PUGRAMA:\n";

#Finally the program opens the second program, that will perform the rest of the analysis

system ("perl pugrama.pl"); 


