lib/validator/transmem.rb in mspire-0.4.7 vs lib/validator/transmem.rb in mspire-0.4.9
- old
+ new
@@ -134,10 +134,10 @@
# what the avg overlap is with transmembrane sequences.
def pep_is_transmem?(pep)
prts = pep.prots
prts_w_keys = 0
sum_of_fractions = 0.0
- prts.each do |prot|
+ prts.each do |prot|
key = @transmem_index.reference_to_key(prot.reference)
ans = @transmem_index.avg_overlap(key, pep.aaseq, :fraction)
if ans
sum_of_fractions += ans
prts_w_keys += 1