my $runnable = Bio::EnsEMBL::Analysis::Runnable::TRF->
  new(
      -query => $slice,
      -program => 'trf',
     );
  $runnable->run;
  my @repeats = @{$runnable->output};

    my ($self, $arg) = @_;

    $self->{'_delta'} = $arg if defined $arg;

    return $self->{'_delta'};

    my ($self, $arg) = @_;

    $self->{'_match'} = $arg if defined $arg;
    
    return $self->{'_match'};

    my ($self, $arg) = @_;

    $self->{'_maxperiod'} = $arg if defined $arg;

    return $self->{'_maxperiod'};

    my ($self, $arg) = @_;

    $self->{'_minscore'} = $arg if defined $arg;

    return $self->{'_minscore'};

    my ($self, $arg) = @_;

    $self->{'_mismatch'} = $arg if defined $arg;
    return $self->{'_mismatch'};

  my ($class,@args) = @_;
  my $self = $class->SUPER::new(@args);
  my ($match, $mismatch, $delta, $pm, $pi, $minscore,
      $maxperiod) = rearrange(['MATCH', 'MISMATCH', 'DELTA', 'PM', 'PI',
                               'MINSCORE', 'MAX_PERIOD'], @args);

  ######################
  #SETTING THE DEFAULTS#
  ######################
  $match = 2 unless(defined $match);
  $mismatch = 5 unless(defined $mismatch);
  $delta = 7 unless(defined $delta);
  $pm = 80 unless(defined $pm);
  $pi = 10 unless(defined $pi);
  $minscore = 40 unless(defined $minscore);
  $maxperiod = 500 unless(defined $maxperiod);
  ######################
  $self->program('trf') if(!$self->program);
  $self->match     ($match)     if defined $match;
  $self->mismatch  ($mismatch)  if defined $mismatch;
  $self->delta     ($delta)     if defined $delta;
  $self->pm        ($pm)        if defined $pm;
  $self->pi        ($pi)        if defined $pi;
  $self->minscore  ($minscore)  if defined $minscore;
  $self->maxperiod ($maxperiod) if defined $maxperiod;
  return $self;
}


# Please use: trf File Match Mismatch Delta PM PI Minscore MaxPeriod 
# [options]
# Where: (all scores are positive)
#   File = sequences input file
#   Match  = matching weight
#   Mismatch  = mismatching penalty
#   Delta = indel penalty
#   PM = match probability (whole number)
#   PI = indel probability (whole number)
#   Minscore = minimum alignment score to report
#   MaxPeriod = maximum period size to report
#   [options] = one or more of the following :
#                -m    masked sequence file
#                -f    flanking sequence
#                -h    avoid html output 

  my ($self, $results) = @_;
  if(!$results){
    $results = $self->resultsfile;
  }
  my $ff = $self->feature_factory;
  open(OUT, "<".$results) or throw("FAILED to open ".$results.
                                   "TRF:parse_results");
  my  @output;
 LINE:while(<OUT>){
    my $seqname;
    if (/Sequence: (\w+)/) {
      $seqname = $1;
    }
    next LINE unless (/^\d/); # ignore introductory lines
    my ($start, $end, $period_size,$copy_number, $consensus_size, 
        $percent_matches, $percent_indels, $score, $A, $C, $G, $T,
        $entropy, $mer) = split;
    if (($score < 50 && $percent_matches >= 80 && $copy_number > 2 && $period_size < 10) || ($copy_number >= 2 && $percent_matches >= 70 && $score >= 50)){
      my $rc = $ff->create_repeat_consensus('trf', 'trf', '', $mer);
      my $rf = $ff->create_repeat_feature($start, $end, 0, $score, 1,
                                         ($end - $start +1), $rc);
      push(@output, $rf);
    }
  }
  print "No tandem repeats found" if(@output == 0);
  $self->output(\@output);
  close(OUT) or throw("FAILED to close ".$results.
                      "TRF:parse_results");
}

1;

    my ($self, $arg) = @_;

    $self->{'_pi'} = $arg if defined $arg;

    return $self->{'_pi'};

    my ($self, $arg) = @_;

    $self->{'_pm'} = $arg if defined $arg;

    return $self->{'_pm'};

  my ($self, $program) = @_;
  if(!$program){
    $program = $self->program;
  }
  my @arg_list = ($self->match, $self->mismatch, $self->delta, $self->pm,
                  $self->pi, $self->minscore, $self->maxperiod);
	$self->options(join(' ', @arg_list));
  throw($program." is not executable TRF::run_analysis ") 
    unless($program && -x $program);
  # Use this command if you are using TRF version 4.0.0
  my $command = $program." ".$self->queryfile." ".$self->options." -d -h";

  # Use this command if your TRF program is previous to the 4.0.0 version
  # my $command = $program." ".$self->queryfile." ".$self->options." -d";

  print "Running analysis ".$command."\n";
  my $exit = system($command);
  foreach my $file(glob $self->queryfile."*"){
    $self->files_to_delete($file);
  }
  $self->resultsfile($self->queryfile.".".join('.', @arg_list).".dat");
  if(!-e $self->resultsfile){
    throw("FAILED to run TRF on ".$self->queryfile." ".
          $self->resultsfile." has not been produced ".
          "TRF:run_analysis");
  }

  Arg [1]   : Bio::EnsEMBL::Analysis::Runnable::CPG
  Arg [2]   : int,
  Function  : container for specified variable. This pod refers to the
  7 methods below match, mismatch, delta, pm, pi, minscore and maxperiod. 
  These are simple containers which dont do more than hold and return an 
  given value nothing is defined
  Returntype: int
  Exceptions: 
  Example   :