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 )