comparison vsnp_statistics.py @ 21:377c1a96aae9 draft

Uploaded
author greg
date Tue, 24 Aug 2021 12:32:50 +0000
parents 1becb6606626
children b34843f09f9f
comparison
equal deleted inserted replaced
20:d4ef7003c539 21:377c1a96aae9
108 else: 108 else:
109 columns = ['Reference', 'FASTQ', 'File Size', 'Mean Read Length', 'Mean Read Quality', 'Reads Passing Q30', 109 columns = ['Reference', 'FASTQ', 'File Size', 'Mean Read Length', 'Mean Read Quality', 'Reads Passing Q30',
110 'Total Reads', 'All Mapped Reads', 'Unmapped Reads', 'Unmapped Reads Percentage of Total', 110 'Total Reads', 'All Mapped Reads', 'Unmapped Reads', 'Unmapped Reads Percentage of Total',
111 'Reference with Coverage', 'Average Depth of Coverage', 'Good SNP Count'] 111 'Reference with Coverage', 'Average Depth of Coverage', 'Good SNP Count']
112 with open(output_file, "w") as outfh: 112 with open(output_file, "w") as outfh:
113 outfh.write("%s\n" % "\t".join(columns)) 113 # Make sure the header starts with a # so
114 # MultiQC can properly handle the output.
115 outfh.write("# %s\n" % "\t".join(columns))
114 line_items = [] 116 line_items = []
115 # Get the current stats and associated files. 117 # Get the current stats and associated files.
116 # Get and output the statistics. 118 # Get and output the statistics.
117 line_items.append(read1_stats.reference) 119 line_items.append(read1_stats.reference)
118 line_items.append(read1_stats.fastq_file) 120 line_items.append(read1_stats.fastq_file)