PiperOrigin-RevId: 368415190
This commit is contained in:
Augustin Zidek
2021-04-14 14:07:38 +01:00
committed by Diego de Las Casas
parent 9c7a288afc
commit 6baaa6950c
+7 -3
View File
@@ -342,12 +342,16 @@ def extract_hmm_profile(hhm_file, sequence, asterisks_replace=0.0):
# The first and the last values in line_values are metadata, skip them. # The first and the last values in line_values are metadata, skip them.
for j, t in enumerate(line_values[2:-1]): for j, t in enumerate(line_values[2:-1]):
aa_profile[count_aa, j] = ( aa_profile[count_aa, j] = (
2**(-float(t) / 1000.) if t != '*' else asterisks_replace) 2**(-float(t) / 1000.0) if t != '*' else asterisks_replace)
count_aa += 1 count_aa += 1
elif len(line_values) == 10: elif len(line_values) == 10:
for j, t in enumerate(line_values): for j, t in enumerate(line_values):
gap_profile[count_gap, j] = ( if j <= 6:
2**(-float(t) / 1000.) if t != '*' else asterisks_replace) gap_profile[count_gap, j] = (
2**(-float(t) / 1000.0) if t != '*' else asterisks_replace)
else:
# Neff_M, Neff_I, and Neff_D are given in units of 0.001.
gap_profile[count_gap, j] = float(t) / 1000.0
count_gap += 1 count_gap += 1
elif not line_values: elif not line_values:
pass pass