Galaxy |

Changeset 7:53bfb3b2c026 (2018-10-18)

Previous changeset 6:e386e916efa1 (2018-10-18) Next changeset 8:5a9a44e23dad (2018-10-19)

Commit message:
Uploaded

modified:
sum_fastqc.pl

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]};