Mercurial > repos > sanbi-uwc > summarize_poliovirus_alignment
comparison summarize_alignment.py @ 2:5c8dfc4d9d68 draft
planemo upload for repository https://github.com/pvanheus/polio_report commit a99e10fec2fac5aae70974c977eb3b362a1a8429-dirty
author | sanbi-uwc |
---|---|
date | Thu, 21 Jul 2022 16:43:55 +0000 |
parents | 7e49c6b19f5e |
children | 7bcfc97b9284 |
comparison
equal
deleted
inserted
replaced
1:f2a7aaf93a1d | 2:5c8dfc4d9d68 |
---|---|
12 parser.add_argument("--summary_output_filename", help="Path to summary output file") | 12 parser.add_argument("--summary_output_filename", help="Path to summary output file") |
13 parser.add_argument("--variant_list_outputs", nargs="+") | 13 parser.add_argument("--variant_list_outputs", nargs="+") |
14 parser.add_argument("--datasets", nargs="+") | 14 parser.add_argument("--datasets", nargs="+") |
15 args = parser.parse_args() | 15 args = parser.parse_args() |
16 | 16 |
17 galaxy_metadata = {} | |
17 summary_output = open(args.summary_output_filename, "w") | 18 summary_output = open(args.summary_output_filename, "w") |
18 for i, json_filename in enumerate(args.datasets): | 19 for i, json_filename in enumerate(args.datasets): |
19 alignment_summary = json.load(open(json_filename)) | 20 alignment_summary = json.load(open(json_filename)) |
20 mismatch_list = ";".join( | 21 mismatch_list = ";".join( |
21 [ | 22 [ |
25 ) | 26 ) |
26 print( | 27 print( |
27 alignment_summary["sample_name"], | 28 alignment_summary["sample_name"], |
28 alignment_summary["best_reference"], | 29 alignment_summary["best_reference"], |
29 alignment_summary["mismatches"], | 30 alignment_summary["mismatches"], |
31 alignment_summary["perc_mismatches"], | |
30 round(alignment_summary["quality"], 2), | 32 round(alignment_summary["quality"], 2), |
31 mismatch_list, | 33 mismatch_list, |
32 sep="\t", | 34 sep="\t", |
33 file=summary_output, | 35 file=summary_output, |
34 ) | 36 ) |