#!/usr/bin/perl -w
#######################################################################
### Huang He                                                        ###
### take fasta file from original fasta by blasted(unfinished) file ###
#######################################################################
($b)=@ARGV;
$c=$b+1;
for($a=1;$a<$c;$a++){
	open OUT,">test$a\.sh" or die "can't creat test$a\.sh";
	print OUT "blastall -p blastx -d /for_blast/ncbi/nr/nr -i /for_blast/temp/predicted_mRNA_$a.fa -o /home/huanghe/research/50000440/Sample_Im1749/treatment_for_mRNA/predicted_mRNA_$a.blastx -v 10 -b 10 -a 12";
}