#!/usr/bin/perl-w
use strict;

  use Bio::Tools::Analysis::Protein::Mitoprot;

  use Bio::PrimarySeq;
  use Bio::SeqIO;
  
my  $inputstream = Bio::SeqIO->new(-file => "myseq2.fa",-format => 'Fasta');



while (my $seqobj = $inputstream->next_seq){   
    # print the sequence   
    print ">";
    print $seqobj->id,"\n";
    
    my $seq= $seqobj-> seq;
    my $id= $seqobj->id;
   

  	my $seq = new Bio::PrimarySeq(-seq=>$seq, -primary_id=>$id); 
 	
 	
     my $mitoprot = Bio::Tools::Analysis::Protein::Mitoprot->new(-seq => $seq); # sequence must be  >!5aa long and start with an M.

	$mitoprot->run();

  die "Could not get a result" unless $mitoprot->status =~ /^COMPLETED/;

  print $mitoprot->result;     # print raw prediction to STDOUT



 
  my $rawdata = my $analysis_object->result;
  
  print "rawdata: $rawdata\n"; 
  
  ## ....................................................
    
   
} # close while