None available.
sub options
{
return ('--model', 'affine:local', '--subopt', 'no', '--bestn', '1'); } |
sub submit_exonerate
{ my ($self, $query, $target, $root_dir, @options) = @_;
my $num_jobs
= XrefMapper::Methods::ExonerateBasic::calculate_num_jobs($query);
$self->jobcount($self->jobcount()+$num_jobs);
my $queryfile = basename($query);
my $targetfile = basename($target);
my $query_alphabet = $queryfile =~ /peptide/ ? 'peptide' : 'dna';
my $outfile = join( '_',
$self->get_class_name() ,
$query_alphabet,
q[%s.map] );
my $options_str = join(" ", @options);
my $exonerate_command = join
( "\\ \n ",
$exonerate_path,
qq[--target $target],
qq[--query $query],
qq[--querychunktotal $num_jobs],
q [--querychunkid %s], q [--showvulgar false],
q [--showalignment FALSE],
q [--ryo "%s"],
$options_str );
my $ryo_format = 'xref:%qi:%ti:%ei:%ql:%tl:%qab:%qae:%tab:%tae:%C:%s\n';
my $shell_command = qq[%s\\ \n| grep '^xref'\\ \n> %s/%s];
foreach( my $i=1; $i<=$num_jobs; $i++ ){
my $command = sprintf( $shell_command,
sprintf($exonerate_command, $i, $ryo_format),
$root_dir,
sprintf($outfile, $i) );
my $status = system( $command );
unless( $status == 0 ){ warn( "Bad exit status ($?) from:\n $command" )}
} } |