[Biopython-dev] More SeqRecord methods

Peter biopython at maubp.freeserve.co.uk
Thu Jun 3 17:24:43 UTC 2010


On Thu, Jun 3, 2010 at 5:43 PM, Peter wrote:
> On Mon, May 31, 2010 at 3:53 PM, Peter wrote:
>
>> ..., how about adding another Seq method to the
>> SeqRecord, the newish ungap method?
>> http://bugzilla.open-bio.org/show_bug.cgi?id=3060
>
> This one I would like some feedback on first. I'm sure the
> implementation could be made much more efficient too.

Maybe I should mention that I also envisage a similar method for the
alignment object, to give a new alignment with any all-gap-columns
removed (perhaps with an optional argument to specify a threshold
for the number of gaps required - defaulting to only removing columns
which are all gaps).

Again, the simplest way to implement this is to re-use the new
alignment slicing and addition features - much as how I did it
for the proposed SeqRecord ungap method.

Peter



More information about the Biopython-dev mailing list