Mercurial > repos > sanbi-uwc > assess_poliovirus_alignment

diff assess_alignment.py @ 8:852e76e7d22a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/pvanheus/polio_report commit a99e10fec2fac5aae70974c977eb3b362a1a8429-dirty
author sanbi-uwc
date Wed, 14 Sep 2022 11:25:43 +0000
parents 9fd6dde72d2e
children acaaf49e2747
line wrap: on
line diff
--- a/assess_alignment.py	Wed Sep 14 11:11:02 2022 +0000
+++ b/assess_alignment.py	Wed Sep 14 11:25:43 2022 +0000
@@ -184,7 +184,7 @@
             best_match_mismatch_list = mismatch_list
             best_match_quality = quality
             best_match_reference = dataset_name
-            best_consensus = open(args.consensi[file_index]).read()
+            best_consensus = open(args.consensi[file_index]).read().replace('>Consensus', f'>{args.sample_name}')
             percent_mismatches = round(min_mismatches / lengths[best_match_reference] * 100, 2)
 
     info = {