#!/usr/bin/perl -w
use Bio::Perl;

use strict;

  # the databases you can get sequences from 
  # are 'swiss', 'genbank', 'genpept', 'embl', and 'refseq'

  my $str = get_sequence('genbank',"NC_001321")->subseq(5,10);
  # my $seq_object = get_sequence('genbank',"ROA1_HUMAN");
  # print $seq_object->subseq(5,10);
  print $str."\n";

  # write_sequence(">nc_001321.fa",'fasta',$seq_object);

  my $gb = new Bio::DB::GenBank();

  # this returns a Seq object :
  my $seq1 = $gb->get_Seq_by_id('MUSIGHBA1');

  # this also returns a Seq object :
  my $seq2 = $gb->get_Seq_by_acc('AF303112');

  # this returns a SeqIO object, which can be used to get a Seq object :
  my $seqio = $gb->get_Stream_by_id(["J00522","AF303112","2981014"]);
  my $seq3 = $seqio->next_seq;

  print "seq3:$seq3";
  write_sequence(">foo",'fasta',$seq3);