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";

    convert_raw_to_fasta($hu, $maid);

    # 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();
}

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_raw_to_fasta($hu, $maid);
    
    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();
}