Is there a function to tabulate exact match seq scoring using pwm weight table for Biopython?
0
0
Entering edit mode
3.3 years ago
ah07seven • 0
x = [Seq("AATTG"), Seq("AGGTC"), Seq("AGGAT"), Seq("AGGCG")]
countSetTable = motifs.create(x)

pwm = countSetTable.counts.normalize(pseudocounts= 0)

seq to matching total score: Seq("AATTG")

result:

total score = 2.5 ( A: 1.0, A: 0.25 , T: 0.25, T: 0.50, G: 0.50)

python biopython • 758 views
ADD COMMENT

Login before adding your answer.

Traffic: 2326 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6