[ensembl-dev] bugs in blastz-chain-net pipeline

Sharon Wei weix at cshl.edu
Thu Feb 2 15:39:35 GMT 2012


Hi Ensembl Compara team,

I recently updated our local blastz pipeline to ensembl-65 code, and 
have found two places where there seem to be a bug.

1. In 
ensembl-analysis/modules/Bio/EnsEMBL/Analysis/Runnable/AlignmentNets.pm
line 394:      $fp->level_id($level_id);
The $fp is a GenomicAlignBlock object, and this object has no function 
level_id(), only GenomicAlign has this function.
So I changed it to
    map {$_->level_id($level_id)} @{$fp->genomic_align_array()};


2. In ensembl-hive/modules/Bio/EnsEMBL/Hive/DBSQL/AnalysisJobAdaptor.pm
In function CreateNewJob(), the following doesn't make sense and I had 
to comment it out to make the pipeline work.
throw("Please specify prev_job object instead of input_job_id if 
available") if ($prev_job_id);

sub CreateNewJob {
   my ($class, @args) = @_;

   return undef unless(scalar @args);

   my ($input_id, $analysis, $prev_job, $prev_job_id, $blocked, 
$semaphore_count, $semaphored_job_id) =
      rearrange([qw(INPUT_ID ANALYSIS PREV_JOB INPUT_JOB_ID BLOCK 
SEMAPHORE_COUNT SEMAPHORED_JOB_ID)], @args);

   throw("must define input_id") unless($input_id);
   throw("must define analysis") unless($analysis);
   throw("analysis must be [Bio::EnsEMBL::Analysis] not a [$analysis]")
     unless($analysis->isa('Bio::EnsEMBL::Analysis'));
   throw("analysis must have adaptor connected to database")
     unless($analysis->adaptor and $analysis->adaptor->db);
   throw("Please specify prev_job object instead of input_job_id if 
available") if ($prev_job_id);   # 'obsolete' message

   $prev_job_id = $prev_job && $prev_job->dbID();

         # if the user did not specifically ask for a new fan, consider 
propagation:


Sharon





More information about the Dev mailing list