[Bioperl-l] AlignIO and SimpleAlign

Jason Stajich jason.stajich at duke.edu
Tue Feb 15 18:20:08 EST 2005

Nope - it's different. It removes every column which has a gap in any 
sequence - so it effectively removes all indels.

If you want the original sequence

  for my $seq ( $aln->each_seq ) {
   my $str = $seq->seq;
  $str =~ s/\-//g;


Jason Stajich
jason.stajich at duke.edu

On Feb 16, 2005, at 3:09 AM, michael watson ((IAH-C)) wrote:

> Hi
> Is it safe to say that if I read in alignment using AlignIO, then call
> remove_gaps on the SimpleAlign object, then access each sequence using
> each_seq, and print each sequence out as fasta using SeqIO, that I will
> have the original sequences that went to make up the alignment 
> (exactly,
> with gaps removed etc)?
> Thanks
> Mick
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PGP.sig
Type: application/pgp-signature
Size: 186 bytes
Desc: This is a digitally signed message part
Url : http://portal.open-bio.org/pipermail/bioperl-l/attachments/20050216/c4b71589/PGP-0001.bin

More information about the Bioperl-l mailing list