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); process_and_parse( "bl2seq -p blastn -i $ltvec_small -j $hu_fa -F F -D 1", "$maid_dir/" . $maid . "_bl2seq.out" ); } 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); process_and_parse( "bl2seq -p blastn -i $ltvec_small -j $hd_fa -F F -D 1", "$maid_dir/" . $maid . "_bl2seq.out" ); }