Mercurial > repos > greg > vsnp_statistics
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) |