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

How to get the score of two alignments? #43

Open
aadimator opened this issue Sep 21, 2020 · 0 comments
Open

How to get the score of two alignments? #43

aadimator opened this issue Sep 21, 2020 · 0 comments

Comments

@aadimator
Copy link

Is there a way to only get the score of two alignments? Let's say that I have two alignments; "MQD--RV--KRP" and "MKKL-K-K-H-P" and I want to get there score. Currently, I'm trying to do:

scoremodel = AffineGapScoreModel(BLOSUM62, gap_open=-3, gap_extend=-1)
res = pairalign(GlobalAlignment(), S1, S2, scoremodel, score_only=true)

But this is giving me an error because of - present in the S1 and S2 sequences. Is there a way to only compute the score of different alignments, without actually aligning them? I have a list of alignments already computed and I want to compute their scores only. Any help/guidance would be much appreciated.

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