# HG changeset patch
# User estrain
# Date 1539914644 14400
# Node ID 53bfb3b2c026d3050a3c483f6d8f1e5923226cd8
# Parent  e386e916efa1cd0f185fd28765ca907c29ab07ec
Uploaded

diff -r e386e916efa1 -r 53bfb3b2c026 sum_fastqc.pl
--- a/sum_fastqc.pl	Thu Oct 18 22:03:49 2018 -0400
+++ b/sum_fastqc.pl	Thu Oct 18 22:04:04 2018 -0400
@@ -16,9 +16,9 @@
 $qscore=~s/\s+//g;
 my(@qlist)=split(/\,/,$qscore);
 
-print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC\%\tMaxN\%\tMeanLen\tMeanQ";
+print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC\%\tMax_N\%\tAvg_Len\tMean_Q";
 foreach(@qlist) {
-  print "\tQ".$_;
+  print "\tQ".$_."\%";
 }
 print "\n";
 
@@ -87,6 +87,7 @@
    return $qmean;
 }
 
+# Calculate mean read Q score
 sub readmean {
    $nreads=shift(@_);
    @qarray=@{$_[0]};
@@ -101,6 +102,7 @@
    return $readq;
 }
 
+# Find position with hights fraction of Ns 
 sub maxn {
    @narray=@{$_[0]};
    my($max_nval)=0;
@@ -111,10 +113,11 @@
        $max_nval=$nval;
      }
    }
-   $max_nval = sprintf("%.2f", 100*$max_nval);
+   $max_nval = sprintf("%.4f", $max_nval);
    return $max_nval;
 }  
 
+# Calculate mean read length
 sub meanlen {
    $nreads=shift(@_);
    @larray=@{$_[0]};