Mercurial > repos > estrain > sum_fastqc
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]};