# HG changeset patch # User estrain # Date 1539774711 14400 # Node ID b769c810924e44544b2faa9102eec7d94a3e2ea8 # Parent 9a9fdfef8b90bacb9221353f43934b04b626c4f6 Added option for multiple Q scores diff -r 9a9fdfef8b90 -r b769c810924e sum_fastqc.pl --- a/sum_fastqc.pl Wed Oct 17 07:11:24 2018 -0400 +++ b/sum_fastqc.pl Wed Oct 17 07:11:51 2018 -0400 @@ -13,8 +13,14 @@ my($inname)=shift(@ARGV); my($qscore)=shift(@ARGV); +$qscore=~s/\s+//g; +my(@qlist)=split(/\,/,$qscore); -print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC\tQ".$qscore."\n"; +print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC"; +foreach(@qlist) { + print "\tQ".$_; +} +print "\n"; foreach (@ARGV) { print_stats($_); @@ -47,8 +53,11 @@ print $pass[1]."\t"; print $nreads[1]."\t"; print $npoor[1]."\t"; - print $gc[1]."\t"; - print qcal($nreads[1],$qscore,\@qlines); + print $gc[1]; + foreach $qs (@qlist) { + print "\t"; + print qcal($nreads[1],$qs,\@qlines); + } print "\n"; }