#!/usr/bin/perl

#use SeqStats;
use Data::Dumper;
use Bio::SeqIO ;
use Bio::Seq ;
use Bio::Tools::SeqStats;
use Bio::Root::Root;  

# Calculate mean Kyte-Doolittle hydropathicity (aka "gravy" score)
  my $prot = Bio::PrimarySeq->new(-seq=>'MSFVLVAPDMLATAAADVVQIGSAVSAGS',-format => 'Fasta',-alphabet=>'protein', -id=>'test');
   $gravy = Bio::Tools::SeqStats->hydropathicity($seqobj);
  print "might be hydropathic" if $gravy > 1;