# HG changeset patch
# User estrain
# Date 1539808151 14400
# Node ID d47775122e782694f1cbc28418bca1a32a83ddea
# Parent  8256c1d0d63b3f82dbede7335e72ab13580e4a49
Uploaded

diff -r 8256c1d0d63b -r d47775122e78 sum_fastqc.pl
--- a/sum_fastqc.pl	Wed Oct 17 11:10:03 2018 -0400
+++ b/sum_fastqc.pl	Wed Oct 17 16:29:11 2018 -0400
@@ -16,7 +16,7 @@
 $qscore=~s/\s+//g;
 my(@qlist)=split(/\,/,$qscore);
 
-print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC";
+print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC\tMeanQ";
 foreach(@qlist) {
   print "\tQ".$_;
 }
@@ -53,7 +53,8 @@
   print $pass[1]."\t";
   print $nreads[1]."\t";
   print $npoor[1]."\t";
-  print $gc[1];
+  print $gc[1]."\t";
+  print readmean($nreads[1],\@qlines);
   foreach $qs (@qlist) {
     print "\t";
     print qcal($nreads[1],$qs,\@qlines);
@@ -77,3 +78,17 @@
    $qmean = sprintf("%.2f", 100 * $sum / $nreads);
    return $qmean;
 }
+
+sub readmean {
+   $nreads=shift(@_);
+   @qarray=@{$_[0]};
+   my($sum) = 0;
+ 
+   foreach $item (@qarray) {
+      my($qval,$q)=split(/\t/,$item);
+      $sum += $q*$qval;
+   }
+
+   $readq = sprintf("%.2f", $sum / $nreads);
+   return $readq;
+}