I think you want something more like this. I'm guessing a little on some of the code.
Otherwise I factored out all the common code, which as far as I can tell is everything but the url.
sub blast_hu_parse {
my($maid,$maid_dir) = @_;
my $url_hu = "http://hu_seq/";
blast_parse($maid,$maid_dir,$url_hu);
}
sub blast_hd_parse {
my($maid,$maid_dir) = @_;
my $url_hd = "http://hd_seq/";
blast_parse($maid,$maid_dir,$url_hd);
}
sub blast_parse{
my($maid,$maid_dir,$url) = @_;
my $h_file = get($url.$maid);
my $ltvec_small = $maid_dir.$maid."Ltvec_small.fa";
#what is $hu where is it declared?
#open ($hu);
#what is $in where is it declared
my $seq=<$in>;
close($in);
open (my $hu, ">".$maid.".fa");
print $hu ">$maid\n$seq";
close($hu);
# syntax
# bl2seq -p blastn -i nucleotide1 -j nucleotide2 -F F -D 1
my $command = "bl2seq -p blastn -i $ltvec_small -j $h_file -F F -D
+
+1";
print $command,"\n";
open OUTPUT, '>', "$maid_dir\\".$maid."_bl2seq.out" ;
STDOUT->fdopen( \*OUTPUT, 'w');
system($command);
bl2seq_parse();
}
update added some comment on the file handling stuff, which appears to have some errors, after looking over the original again
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.