sub blast_parse{
	my($maid,$maid_dir) = @_;
	my $url_hu = "http://hu_seq/";
	my $hu = get($url_hu.$maid);
	my $ltvec_small = $maid_dir.$maid."Ltvec_small.fa";
#=cut
# convert the raw sequence to fasta
	open ($hu);
	my $seq=<$in>;
	close($in);
	open ($hu, ">".$maid.".fa");
	print $hu_fa ">$maid\n$seq";
	close($hu_fa);
#=cut
	# syntax
	# bl2seq -p blastn -i nucleotide1 -j nucleotide2 -F F -D 1 
	my $command = "bl2seq -p blastn -i $ltvec_small -j $hu_fa -F F -D 1";
	print $command,"\n";
	open OUTPUT, '>',  "$maid_dir\\".$maid."_bl2seq.out" ;
	STDOUT->fdopen( \*OUTPUT, 'w');
	system($command); 
	bl2seq_parse();
}

# the next method

sub blast_hd_parse{
	my($maid,$maid_dir) = @_;
	my $url_hd = "http://hd_seq/";
	my $hd = get($url_hd.$maid);
	my $ltvec_small = $maid_dir.$maid."Ltvec_small.fa";
# convert the raw sequence to fasta	
	open ($hu);
	my $seq=<$in>;
	close($in);
	open ($hu, ">".$maid.".fa");
	print $hu_fa ">$maid\n$seq";
	close($hu_fa);
	
	my $command = "bl2seq -p blastn -i $ltvec_small -j $hd_fa -F F -D 1";
	print $command,"\n";
	open OUTPUT, '>',  "$maid_dir\\".$maid."_bl2seq.out" ;
	STDOUT->fdopen( \*OUTPUT, 'w');
	system($command); 
	bl2seq_parse();
}