changeset 7:53bfb3b2c026 draft

Uploaded
author estrain
date Thu, 18 Oct 2018 22:04:04 -0400
parents e386e916efa1
children 5a9a44e23dad
files sum_fastqc.pl
diffstat 1 files changed, 6 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- 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]};