[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;
  $seq->seq($str);
  }

-jason


--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/

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