[Biopython-dev] [Biopython - Bug #3387] (New) Generic per column annotation from stockholm alignment are not stored in alignment object

redmine at redmine.open-bio.org redmine at redmine.open-bio.org
Wed Oct 17 13:27:18 UTC 2012


Issue #3387 has been reported by saverio vicario.

----------------------------------------
Bug #3387: Generic per column annotation from stockholm alignment are not stored in alignment object
https://redmine.open-bio.org/issues/3387

Author: saverio vicario
Status: New
Priority: Normal
Assignee: Biopython Dev Mailing List
Category: Main Distribution
Target version: 
URL: 


Stockholm format includes 4 types of annotations
#=GF <feature> <Generic per-File annotation, free text>
#=GC <feature> <Generic per-Column annotation, exactly 1 char per column>
#=GS <seqname> <feature> <Generic per-Sequence annotation, free text>
#=GR <seqname> <feature> <Generic per-Sequence AND per-Column markup, exactly 1 char per column>
GC and GF annotation are not pickup by AlignIO and not supported in Bio.Align.MultipleSeqAlignment because no annotation is available at alignment level. In fact Bio.Align.MultipleSeqAlignment.annotations or Bio.Align.MultipleSeqAlignment.letter_annotations do not exist, only Bio.Align.MultipleSeqAlignment._annotations that is generated from the single records annotations and letter_annotations.

GC annotation in stockholm contain the quality score of the sites (columns of the alignment) that is a quite important parameters to decide if to trim the sites or not. 




-- 
You have received this notification because you have either subscribed to it, or are involved in it.
To change your notification preferences, please click here and login: http://redmine.open-bio.org




More information about the Biopython-dev mailing list