[Biopython-dev] [Biopython] Bio.Sequencing.Ace
Jose Blanca
jblanca at btc.upv.es
Tue Jun 30 04:31:06 EDT 2009
> What I was thinking of was a contig class as an alignment subclass,
> holding a list of SeqRecord objects and offsets. The consensus might
> just be one element of this list - but could be handled specially. This
> sounds simpler than having to introduce a whole new object system,
> related to but different to SeqFeature objects. However, I don't yet
> have a sample implementation to demonstrate this.
I thought about that implementation and I created some code. The problem I
found with that approach is that the contig class code got too messy. Take
into account that besides the offset you also need the masks and that some
sequences could be reversed. That's why I decided to split the part that
calculates the offset and the mask into a separate class.
> One important thing I think we should do BEFORE adding any contig
> class to Biopython, is get it working with at least one other contig file
> format in addition to Ace. I don't want to end up with a class which
> is too specialised for how ace contigs work.
>
> Peter
Well, In fact my contig class is modeled after the caf file format. The ace
parsing was just an afterthought, my primary interest was the caf format.
--
Jose M. Blanca Postigo
Instituto Universitario de Conservacion y
Mejora de la Agrodiversidad Valenciana (COMAV)
Universidad Politecnica de Valencia (UPV)
Edificio CPI (Ciudad Politecnica de la Innovacion), 8E
46022 Valencia (SPAIN)
Tlf.:+34-96-3877000 (ext 88473)
More information about the Biopython-dev
mailing list