[Bioperl-l] Remove gap only columns from alignment

[Jason Stajich]

Passing in a second parameter, will only remove the all-gapped column
Like this:
$aln->remove_gaps('-',1);

or
$aln->remove_gaps(undef,1);


  remove_gaps

  Title     : remove_gaps
  Usage     : $aln2 = $aln->remove_gaps('-'[,$all_gaps_columns])
  Function  : Creates an aligment with gaps removed
  Returns   : a Bio::SimpleAlign object
  Args      : a gap character(optional) if no specified,
                 taken from $self->gap_char,
                 optional $all_gaps_columns flag
                 indicates that only all-gaps columns should be deleted



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

On Feb 16, 2005, at 3:21 AM, Goel, Manisha wrote:

> Hi All,
> I am currently using Bio::SimpleAlign to manipulate my alignments.
> I need to remove columns which are represented by gap characters in all
> seqs (similar to remove gap only columns in clustal), but the
> remove_gaps command seems to remove all columns even if a single
> sequence in the alignment has a gap at that position.
> Is there any other workaround to get this done ?
>
> Thanks,
> -Manisha
>
> _______________________________________________
> 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/2ac8da1d/PGP.bin