@@ -210,6 +210,7 @@ def main():
210210
211211
212212 args = parser .parse_args ()
213+
213214 json_dict , criteria_labels = parse_qc_file (args .input_file )
214215
215216 # Write JSON
@@ -232,6 +233,7 @@ def main():
232233 else :
233234 print ("[!] Skipping histogram section — required length_count.txt files not found." )
234235 args .hist_out = None
236+
235237 # GC content files
236238 base1 = os .path .join (qc_stats_dir , "qa.QC.base_content.txt" )
237239 base2 = os .path .join (qc_stats_dir , "QC.base_content.txt" )
@@ -244,6 +246,7 @@ def main():
244246 print ("[!] Skipping GC content section — required base_content.txt files not found." )
245247 args .gc1_out = None
246248 args .gc2_out = None
249+
247250 # ATCG composition plots
248251 base_matrix1 = os .path .join (qc_stats_dir , "qa.QC.base.matrix" )
249252 base_matrix2 = os .path .join (qc_stats_dir , "QC.base.matrix" )
@@ -270,7 +273,7 @@ def main():
270273 combined_qual = faqcs_quality_histogram .combine_quality_histograms (qh_fig1 , qh_fig2 , qa_annotation , main_annotation , qh_min1 , qh_max1 , qh_min2 , qh_max2 )
271274 combined_qual .write_html (args .qual_out )
272275 print (f"[✓] Quality histogram written to { args .qual_out } " )
273- print ( "past quality histograms \n " )
276+
274277 # Quality boxplot and 3D plots
275278 qa_matrix = os .path .join (qc_stats_dir , "qa.QC.quality.matrix" )
276279 trim_matrix = os .path .join (qc_stats_dir , "QC.quality.matrix" )
@@ -290,7 +293,7 @@ def main():
290293 qbar_fig2 , qbar_anno2 = faqcs_quality_plots .quality_count_histogram (trim_matrix , qh_max1 , "Trimmed Reads Q score" , "" )
291294 faqcs_quality_plots .combine_quality_histograms (qbar_fig1 , qbar_fig2 , qbar_anno1 , qbar_anno2 ).write_html (args .qhist_out )
292295 print (f"[✓] Quality score histogram written to { args .qhist_out } " )
293- print ( "past quality boxplots and histograms" )
296+
294297 # Merge into final report
295298 sections = [("QC Summary Plots" , args .html_out )]
296299 if os .path .isfile (args .hist_out ):
@@ -313,5 +316,6 @@ def main():
313316 sections .append (("Quality Score Histogram" , args .qhist_out ))
314317
315318 merge_html_plots (sections , args .final_out )
319+
316320if __name__ == "__main__" :
317- main ()
321+ main ()
0 commit comments