#!/usr/bin/perl -w use strict; use Bio::SeqIO; use analysis; my $filename=<>; my %seqhash=(); my $sequence=''; my $seqname=''; my $seq_no=0; my $stream = Bio::SeqIO->new(-file => $filename, -format => 'GenBank'); while (my $seqObj = $stream->next_seq() ) { my $seq = $seqObj->seq(); my $display_id = $seqObj->display_id; $seqhash{$sequence}=$seq; $seqhash{$seqname}=$display_id; ++$seq_no; foreach my $key (keys %seqhash) { print "\n\nSeq No.".$seq_no.":=".$seqhash{$seqname}; my $nucleotides=analysis->nucleotide_count($seqhash{seqname},$seq); print "Amino Acid = "; my $aminoacid=analysis->aminoanalysis($seqhash{seqname},$seq); } }