Mercurial > repos > iuc > hyphy_fade

changeset 28:7ae1a772f08e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hyphy/ commit 9fa2234b56facaf70fce12e5c60638d801997594"
author iuc
date Fri, 18 Jun 2021 15:25:05 +0000
parents d71d20c0ac09
children fccaa7ba05fd
files scripts/hyphy_summary.py
diffstat 1 files changed, 5 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/hyphy_summary.py	Wed Jun 09 17:26:30 2021 +0000
+++ b/scripts/hyphy_summary.py	Fri Jun 18 15:25:05 2021 +0000
@@ -336,6 +336,9 @@
                 report = self.annotation_json[self.include_in_annotation[i]]
                 report['fade'] = {}
                 for residue, info in self.fade['MLE']['content'].items():
+                    # Skip anything past the end of this section.
+                    if i >= len(info['0']):
+                        continue
                     if len(residue) == 1:
                         report['fade'][residue] = {'rate': info['0'][i][1], 'BF': info['0'][i][-1]}
 
@@ -373,11 +376,13 @@
 
                 def ignore_record(x):
                     pass
+
                 for s in self.ref_genes:
                     _align_par(SeqRecord(Seq(s[1]), id=s[0]), [SeqRecord(Seq(ref_seq), id='ref')],
                                self.score_matrix_, False, False, 0.8, ignore_record, output_record)
                     if (self.aligned_str is not None):
                         break
+
                 self.ref_map = self.aligned_str.seq.strip('-')
                 c = 0
                 i = 0