Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Pfam uses #=GC seq_cons, Easel expects #=GC RF #18

Open
cryptogenomicon opened this issue Jan 18, 2017 · 0 comments
Open

Pfam uses #=GC seq_cons, Easel expects #=GC RF #18

cryptogenomicon opened this issue Jan 18, 2017 · 0 comments

Comments

@cryptogenomicon
Copy link
Member

Consensus column annotation in Pfam is on a #=GC seq_cons line. Easel expects reference annotation on a #=GC RF line. The Pfam seq_cons lines are handled by Easel as "unparsed annotation", which is usually fine. An example where it's not fine is in esl_msa_SequenceSubset(), which only propagates recognized GC annotations, not unparsed ones, for fear of propagating annotation that would not be valid for a subsetted alignment. Elena and Grey ran across this.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant