sub process_and_parse {
    my $command = shift;
    my $output_file = shift or die "Missing argument(s)!";

    print $command,"\n";
    open OUTPUT, '>',  $output_file or die $!;
    STDOUT->fdopen( \*OUTPUT, 'w');
    system($command); 
    bl2seq_parse();
}