[Bioperl-l] ace to msf format?

[Wes Barris]

Can anyone give me a hint as to how I could use bioperl to read in
an ACE assembly and write out an MSF formatted alignment?  This shows
what I have figured out so far:

#!/usr/local/bin/perl -w
#
use strict;
use Bio::Assembly::IO;
#
my $usage = "Usage: $0 <infile.ace>\n";
my $infile = shift or die $usage;

my $io = new Bio::Assembly::IO(-file=>$infile, -format=>'ace');
my $assembly = $io->next_assembly;

my $aln = $assembly->all_contigs();

-- 
Wes Barris
E-Mail: Wes.Barris at csiro.au