XrefMapper::Methods LocalExonerateGappedBest1
Included librariesPackage variablesGeneral documentationMethods
Toolbar
WebCvsRaw content
Package variables
No package variables defined.
Included modules
File::Basename
XrefMapper::Methods::ExonerateBasic
Inherit
XrefMapper::Methods::ExonerateBasic
Synopsis
No synopsis!
Description
No description!
Methods
options
No description
Code
query_identity_threshold
No description
Code
submit_depend_job
No description
Code
submit_exonerate
No description
Code
target_identity_threshold
No description
Code
Methods description
None available.
Methods code
optionsdescriptionprevnextTop
sub options {
  return ('--model', 'affine:local', '--subopt', 'no', '--bestn', '1');
}
query_identity_thresholddescriptionprevnextTop
sub query_identity_threshold {
  #return 90;
return 50;
}
submit_depend_jobdescriptionprevnextTop
sub submit_depend_job {
  die;
}

1;
}
submit_exoneratedescriptionprevnextTop
sub submit_exonerate {
  my ($self, $query, $target, $root_dir, @options) = @_;

  # Exonerate can run individual chunks of a large job.
# Determine the number of chunks that will be used, and add to total
my $num_jobs = XrefMapper::Methods::ExonerateBasic::calculate_num_jobs($query); $self->jobcount($self->jobcount()+$num_jobs); # Determine the output filename
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] # Substitute the chunk ID
); # Build a template representing the exonerate command
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], # Substitute the chunk ID
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]; # Run exonerate for each job chunk
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" )} }
}
target_identity_thresholddescriptionprevnextTop
sub target_identity_threshold {
  #return 90;
return 50;
}
General documentation
No general documentation available.