annotate sum_fastqc.pl @ 4:d47775122e78 draft

Uploaded
author estrain
date Wed, 17 Oct 2018 16:29:11 -0400
parents 8256c1d0d63b
children 7df018757d26
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
1 #!/usr/bin/perl
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
2
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
3 ####################################################
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
4 ##
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
5 ## sum_fastqc.pl
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
6 ##
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
7 ## Errol Strain (estrain@gmail.com)
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
8 ##
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
9 ## Description: Takes raw FASTQC output and produces
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
10 ## simple table summary
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
11 ##
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
12 ####################################################
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
13
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
14 my($inname)=shift(@ARGV);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
15 my($qscore)=shift(@ARGV);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
16 $qscore=~s/\s+//g;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
17 my(@qlist)=split(/\,/,$qscore);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
18
4
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
19 print "Input\tFile\tFastQC\tPass-Fail\tReads\tPoor_Reads\tGC\tMeanQ";
3
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
20 foreach(@qlist) {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
21 print "\tQ".$_;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
22 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
23 print "\n";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
24
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
25 foreach (@ARGV) {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
26 print_stats($_);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
27 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
28
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
29 sub print_stats {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
30 $infile = shift;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
31 # First 10 lines of raw FASTQC contain basic overview
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
32 @sumlines=`head -n 10 $infile`;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
33 chomp(@sumlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
34
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
35 # Sequence level Q scores are buried in the middle of the file
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
36 @qlines=`awk '/#Quality\tCount/,/>>END_MODULE/' $infile | head -n -1 | tail -n +2`;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
37 chomp(@qlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
38
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
39 @fastqc = split(/[\n\t]/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
40 @pass = split(/\t/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
41 shift(@sumlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
42 @fn = split(/\t/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
43 shift(@sumlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
44 shift(@sumlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
45 @nreads = split(/\t/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
46 @npoor = split(/\t/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
47 shift(@sumlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
48 @gc = split(/\t/,shift(@sumlines));
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
49
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
50 print $inname."\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
51 print $fn[1]."\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
52 print $fastqc[1]."\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
53 print $pass[1]."\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
54 print $nreads[1]."\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
55 print $npoor[1]."\t";
4
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
56 print $gc[1]."\t";
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
57 print readmean($nreads[1],\@qlines);
3
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
58 foreach $qs (@qlist) {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
59 print "\t";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
60 print qcal($nreads[1],$qs,\@qlines);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
61 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
62 print "\n";
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
63 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
64
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
65 # Sum reads w/ Q scores > cutoff and divide by number of reads
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
66 sub qcal {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
67 $nreads=shift(@_);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
68 $cutoff=shift(@_);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
69 @qarray=@{$_[0]};
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
70 $sum = 0;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
71
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
72 foreach $item (@qarray) {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
73 my($qval,$q)=split(/\t/,$item);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
74 if($qval>=$cutoff) {
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
75 $sum += $q;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
76 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
77 }
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
78 $qmean = sprintf("%.2f", 100 * $sum / $nreads);
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
79 return $qmean;
8256c1d0d63b Uploaded
estrain
parents:
diff changeset
80 }
4
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
81
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
82 sub readmean {
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
83 $nreads=shift(@_);
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
84 @qarray=@{$_[0]};
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
85 my($sum) = 0;
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
86
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
87 foreach $item (@qarray) {
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
88 my($qval,$q)=split(/\t/,$item);
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
89 $sum += $q*$qval;
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
90 }
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
91
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
92 $readq = sprintf("%.2f", $sum / $nreads);
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
93 return $readq;
d47775122e78 Uploaded
estrain
parents: 3
diff changeset
94 }