annotate bin/sRNAPipe.pl @ 50:372ce358ef09 draft

Uploaded
author pierre.pouchin
date Thu, 06 Sep 2018 11:25:35 -0400
parents 4bc00caa60b4
children 9645d995fb3c
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
40
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
1 #!/usr/bin/perl
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
2 use strict;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
3 use warnings;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
4 use Getopt::Long;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
5 use Parallel::ForkManager;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
6 use File::Basename;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
7 use File::Copy::Recursive qw( dircopy );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
8 use POSIX;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
9 use FindBin;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
10 use lib $FindBin::Bin;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
11 use resize qw ( size_distribution );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
12 use subgroups qw (subgroups );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
13 use ppp qw ( ping_pong_partners );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
14 use Rcall qw (pie_chart bg_to_png );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
15 use align qw ( to_build get_unique sam_count sam_count_mis sam_sorted_bam rpms_rpkm rpms_rpkm_te BWA_call get_fastq_seq extract_sam sam_to_bam_bg );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
16 use html qw ( main_page details_pages menu_page ppp_page );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
17 use File::Copy;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
18
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
19 my ( @fastq, @fastq_n, $dir, $min, $max, $mis, $misTE, $help, $Pcheck, $mapnumf, $html_out);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
20 my ( $ref, $tRNAs, $rRNAs, $snRNAs, $miRNAs, $transcripts, $TE );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
21 my ( $si_min, $si_max, $pi_min, $pi_max );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
22 my ( $build_index, $build_tRNAs, $build_rRNAs, $build_snRNAs, $build_miRNAs, $build_transcripts, $build_TE );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
23 my $max_procs = 8;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
24
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
25 ( $build_index, $build_tRNAs, $build_rRNAs, $build_snRNAs, $build_miRNAs, $build_transcripts, $build_TE ) = (0,0,0,0,0,0,0);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
26 ( $min, $max, $mis, $misTE, $si_min, $si_max, $pi_min, $pi_max, $dir ) = ( 18, 29, 0, 3, 21, 21, 23, 29 );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
27 $Pcheck ='true';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
28
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
29 GetOptions (
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
30 "fastq=s" => \@fastq,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
31 "fastq_n=s" => \@fastq_n,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
32 "dir=s" => \$dir,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
33 "min:i" => \$min,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
34 "max:i" => \$max,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
35 "si_min:i" => \$si_min,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
36 "si_max:i" => \$si_max,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
37 "pi_min:i" => \$pi_min,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
38 "pi_max:i" => \$pi_max,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
39 "mis:i" => \$mis,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
40 "misTE:i" => \$misTE,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
41 "html:s" => \$html_out,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
42 "PPPon:s" => \$Pcheck,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
43 "help" => \$help,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
44 "ref:s" => \$ref,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
45 "tRNAs:s" => \$tRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
46 "rRNAs:s" => \$rRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
47 "snRNAs:s" => \$snRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
48 "miRNAs:s" => \$miRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
49 "transcripts:s" => \$transcripts,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
50 "TE:s" => \$TE,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
51 "build_index" => \$build_index,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
52 "build_tRNAs" => \$build_tRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
53 "build_snRNAs" => \$build_snRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
54 "build_miRNAs" => \$build_miRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
55 "build_transcripts" => \$build_transcripts,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
56 "build_rRNAs" => \$build_rRNAs,
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
57 "build_TE" => \$build_TE
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
58 );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
59
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
60 my $fq_collection = 'fastq_dir/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
61 mkdir $dir; mkdir $fq_collection;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
62 $dir = $dir.'/' unless $dir =~ /\/$/;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
63 mkdir $dir.'/css';mkdir $dir.'/js';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
64 dircopy( $FindBin::Bin.'/css', $dir.'/css' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
65 dircopy( $FindBin::Bin.'/js', $dir.'/js' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
66
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
67 my $file = $dir.'report.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
68 open my $report, '>', $file or die "Cannot open $file $!\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
69
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
70 my @toBuild = ( [$build_index, \$ref], [$build_tRNAs, \$tRNAs], [$build_rRNAs, \$rRNAs], [$build_snRNAs, \$snRNAs], [$build_miRNAs, \$miRNAs], [$build_transcripts, \$transcripts], [$build_TE, \$TE] );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
71 to_build ( \@toBuild, $report, $dir );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
72
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
73 my $proc_child = ceil($max_procs / scalar(@fastq));
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
74 my $proc_grand_child = ceil($proc_child/4);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
75 my $pm = Parallel::ForkManager->new($max_procs);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
76 my $pm2 = Parallel::ForkManager->new($proc_grand_child);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
77
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
78 $pm->run_on_finish( sub {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
79 my ($pid, $exit_code, $ident) = @_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
80 print $report "Fastq fork $ident just finished ".
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
81 "with PID $pid and exit code: $exit_code\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
82 die "Something went wrong!\n" if $exit_code != 0;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
83 });
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
84 $pm->run_on_start( sub {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
85 my ($pid,$ident)=@_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
86 print $report "Fastq fork : $ident started, pid: $pid\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
87 });
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
88 $pm2->run_on_finish( sub {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
89 my ($pid, $exit_code, $ident) = @_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
90 print $report "** Subgroup fork $ident just finished ".
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
91 "with PID $pid and exit code: $exit_code\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
92 die "Something went wrong!\n" if $exit_code != 0;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
93 });
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
94 $pm2->run_on_start( sub {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
95 my ($pid,$ident)=@_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
96 print $report "** Subgroup fork $ident started, pid: $pid\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
97 });
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
98
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
99
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
100 foreach my $child ( 0 .. $#fastq )
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
101 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
102 my @suffix = ('.fastq', '.fastq.gz,', '.fq', '.fq.gz', 'ref', '.dat', '.fa','.fas','.fasta', '.txt');
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
103 my ( $name, $path, $suffix ) = fileparse( $fastq[$child], @suffix );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
104 my ( $ref_name, $ref_path, $ref_suffix ) = fileparse( $ref, @suffix );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
105 my ( $TE_name, $TE_path, $TE_suffix ) = fileparse( $TE, @suffix );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
106 my ( $ex_name, $ex_path, $ex_suffix ) = fileparse( $transcripts, @suffix );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
107
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
108 $pm->start($fastq[$child]) and next;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
109
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
110 my $dir_fq = $dir.$fastq_n[$child].'/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
111 mkdir $dir_fq;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
112
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
113 my $gen_dir = $dir_fq.'genome/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
114 mkdir $gen_dir;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
115
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
116 my $size_dir = $dir_fq.'size/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
117 mkdir $size_dir;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
118
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
119 my $fastq_resized = $dir_fq.$name.'_'.$min.'-'.$max.'.fastq';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
120 size_distribution ( $fastq[$child], $fastq_resized, $size_dir, $min, $max );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
121
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
122 my $sam_genome = $gen_dir.$fastq_n[$child].'_'.$min.'-'.$max.'.sam';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
123 my $sam_genome_unique = $gen_dir.$fastq_n[$child].'_'.$min.'-'.$max.'_unique.sam';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
124 my $fastq_prefix = $gen_dir.$fastq_n[$child].'_'.$min.'-'.$max;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
125
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
126 BWA_call ( $ref, $fastq_resized, $sam_genome, $mis, $proc_child, $report );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
127 my ( $fai_ref_hashP, $ma, $ma_uni ) = get_unique ( $sam_genome, $sam_genome_unique, $gen_dir, $fq_collection.$fastq_n[$child], 1, $report );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
128
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
129 die "No Reads mapped on the genome reference!\n" if $ma == 0;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
130 my $scale = 1000000 / $ma;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
131 sam_to_bam_bg ( $sam_genome_unique, $scale, $proc_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
132 sam_to_bam_bg ( $sam_genome, $scale, $proc_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
133
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
134 my $Gviz_dir = $gen_dir.'Gviz/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
135 my $fai_file = $gen_dir.'fai';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
136 mkdir $Gviz_dir;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
137 my $Gviz_dir_rand = $Gviz_dir.'rand/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
138 mkdir $Gviz_dir_rand;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
139 my $Gviz_dir_uni = $Gviz_dir.'unique/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
140 mkdir $Gviz_dir_uni;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
141
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
142 open my $gfai, '>', $fai_file;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
143 foreach my $k ( sort keys %{$fai_ref_hashP} )
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
144 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
145 print $gfai "$k\t$fai_ref_hashP->{$k}\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
146 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
147 close $gfai;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
148 bg_to_png ( $fai_file, $fastq_prefix.'_unique_plus.bedgraph', $fastq_prefix.'_unique_minus.bedgraph', $Gviz_dir_uni, 'Mb' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
149 bg_to_png ( $fai_file, $fastq_prefix.'_plus.bedgraph', $fastq_prefix.'_minus.bedgraph', $Gviz_dir_rand, 'Mb' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
150
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
151 my $group_dir = $dir_fq.'subgroups/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
152 my $fastq_uni = $fq_collection.$fastq_n[$child].'_unique_mappers.fastq';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
153 my $fastq_all = $fq_collection.$fastq_n[$child].'_all_mappers.fastq';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
154 my ($bo, $mi, $pi) = subgroups ( $fastq_all, $group_dir, $mis, $misTE, $proc_child, $tRNAs, $rRNAs, $snRNAs, $miRNAs, $transcripts, $TE, $si_min, $si_max, $pi_min, $pi_max, $report);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
155
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
156 pie_chart($group_dir);
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
157
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
158 open (my $dupnum, $gen_dir.'dup_mapnum.txt') || die "cannot open dup_mapnum.txt $!";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
159 my %dupnum_genome;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
160 my $header = <$dupnum>;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
161 while (<$dupnum>)
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
162 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
163 chomp $_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
164 my @dupline = split /\t/, $_;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
165 $dupnum_genome{$dupline[0]} = [$dupline[1], $dupline[2]];
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
166 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
167 close $dupnum;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
168
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
169 my $mi_sam = $group_dir.'miRNAs.sam';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
170 mkdir $group_dir.'miRNAs/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
171 my $mi_count_file = $group_dir.'miRNAs/miRNAs_reads_counts.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
172 my ( $mi_count, $mi_ref_size ) = sam_count ( $mi_sam );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
173
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
174 rpms_rpkm( $mi_count, $mi_ref_size, $ma, $mi_count_file, $pi, $mi, $bo );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
175
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
176 my ( $sam_transcripts, $sam_TEs ) = ( $group_dir.'transcripts.sam', $group_dir.'TEs.sam' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
177 my @types = ($group_dir.'bonafide_reads.fastq', $group_dir.'miRNAs.fastq', $group_dir.'siRNAs.fastq', $group_dir.'piRNAs.fastq' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
178 my @types_names = ('bonafide_reads', 'miRNAs', 'siRNAs', 'piRNAs');
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
179 foreach my $grand_child ( 0 .. $#types )
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
180 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
181 my $type_dir = $group_dir.$types_names[$grand_child].'/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
182 my $type_prefix = $types_names[$grand_child].'-';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
183 mkdir $type_dir;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
184 $pm2->start($types[$grand_child]) and next;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
185 my ( $type_sam_genome, $type_sam_TEs, $type_sam_transcripts ) = ( $type_dir.$type_prefix.'genome.sam', $type_dir.$type_prefix.'TEs.sam', $type_dir.$type_prefix.'transcripts.sam' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
186 my ( $type_sam_uni_genome, $type_sam_uni_TEs, $type_sam_uni_transcripts ) = ( $type_dir.$type_prefix.'genome_unique.sam', $type_dir.$type_prefix.'TEs_unique.sam', $type_dir.$type_prefix.'transcripts_unique.sam' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
187 my ( $type_uni_genome_fastq, $type_uni_TEs_fastq, $type_uni_transcripts_fastq ) = ( $fq_collection.$fastq_n[$child].'-'.$type_prefix.'genome_uni.fastq', $fq_collection.$fastq_n[$child].'-'.$type_prefix.'TEs_uni.fastq', $fq_collection.$fastq_n[$child].'-'.$type_prefix.'transcripts_uni.fastq');
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
188 my ( $type_genome_fastq, $type_TEs_fastq, $type_transcripts_fastq ) = ( $fq_collection.$fastq_n[$child].'-'.$type_prefix.'genome.fastq', $fq_collection.$fastq_n[$child].'-'.$type_prefix.'TEs.fastq', $fq_collection.$fastq_n[$child].'-'.$type_prefix.'transcripts.fastq');
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
189 my $type_sequence_hashP = get_fastq_seq ( $types[$grand_child] );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
190
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
191 if ( $grand_child == 1 )
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
192 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
193 BWA_call ( $TE, $types[$grand_child], $type_sam_TEs, $misTE, $proc_child, $report );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
194 BWA_call ( $transcripts, $types[$grand_child], $type_sam_transcripts, $mis, $proc_child, $report );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
195 BWA_call ( $ref, $types[$grand_child], $type_sam_genome, $mis, $proc_child, $report );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
196 extract_sam ( undef, $type_sam_TEs, $type_sam_TEs, $type_sam_uni_TEs, $type_uni_TEs_fastq, $type_uni_TEs_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
197 extract_sam ( undef, $type_sam_transcripts, $type_sam_transcripts, $type_sam_uni_transcripts, $type_transcripts_fastq, $type_uni_transcripts_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
198 extract_sam ( undef, $type_sam_genome, $type_sam_genome, $type_sam_uni_genome, $type_genome_fastq, $type_uni_genome_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
199 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
200 else
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
201 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
202 extract_sam ( $type_sequence_hashP, $sam_TEs, $type_sam_TEs, $type_sam_uni_TEs, $type_TEs_fastq, $type_uni_TEs_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
203 extract_sam ( $type_sequence_hashP, $sam_transcripts, $type_sam_transcripts, $type_sam_uni_transcripts, $type_transcripts_fastq, $type_uni_transcripts_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
204 extract_sam ( $type_sequence_hashP, $sam_genome, $type_sam_genome, $type_sam_uni_genome, $type_genome_fastq, $type_uni_genome_fastq );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
205 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
206
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
207 my $ex_count_file = $type_dir.$type_prefix.'transcripts_reads_counts.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
208 my ( $ex_count, $ex_ref_size ) = sam_count_mis ( $type_sam_transcripts );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
209 rpms_rpkm_te( $ex_count, $ex_ref_size, $ma, $ex_count_file, $pi, $mi, $bo );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
210
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
211 my ( $TEs_count, $TEs_ref_size, $TEs_count_NoM, $TEs_count_M ) = sam_count_mis ( $type_sam_TEs );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
212 my $TEs_count_file = $type_dir.$type_prefix.'TEs_reads_counts.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
213 my $TEs_count_file_M = $type_dir.$type_prefix.'TEs_reads_counts_mismatches.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
214 my $TEs_count_file_noM = $type_dir.$type_prefix.'TEs_reads_counts_nomismatches.txt';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
215 rpms_rpkm_te( $TEs_count, $TEs_ref_size, $ma, $TEs_count_file, $pi, $mi, $bo );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
216 rpms_rpkm_te( $TEs_count_NoM, $TEs_ref_size, $ma, $TEs_count_file_noM, $pi, $mi, $bo );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
217 rpms_rpkm_te( $TEs_count_M, $TEs_ref_size, $ma, $TEs_count_file_M, $pi, $mi, $bo );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
218
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
219 sam_to_bam_bg ( $type_sam_TEs, $scale, $grand_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
220 sam_sorted_bam ( $type_sam_transcripts, $grand_child ); sam_sorted_bam ( $type_sam_uni_transcripts, $grand_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
221 sam_sorted_bam ( $type_sam_uni_TEs, $grand_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
222
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
223 my $Gviz_TEs = $type_dir.'Gviz_TEs/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
224 mkdir $Gviz_TEs;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
225 bg_to_png ( $group_dir.'TEs.fai', $type_dir.$type_prefix.'TEs_plus.bedgraph', $type_dir.$type_prefix.'TEs_minus.bedgraph', $Gviz_TEs, 'Kb' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
226
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
227 my $Gviz_genome= $type_dir.'Gviz_genome/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
228 my $Gviz_genome_rand = $Gviz_genome.'rand/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
229 my $Gviz_genome_uni = $Gviz_genome.'unique/';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
230 mkdir $Gviz_genome; mkdir $Gviz_genome_uni; mkdir $Gviz_genome_rand;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
231
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
232 sam_to_bam_bg ( $type_sam_genome, $scale, $grand_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
233 sam_to_bam_bg ( $type_sam_uni_genome, $scale, $grand_child );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
234
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
235 bg_to_png ( $fai_file, $type_dir.$type_prefix.'genome_unique_plus.bedgraph', $type_dir.$type_prefix.'genome_unique_minus.bedgraph', $Gviz_genome_uni, 'Mb' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
236 bg_to_png ( $fai_file, $type_dir.$type_prefix.'genome_plus.bedgraph', $type_dir.$type_prefix.'genome_minus.bedgraph', $Gviz_genome_rand, 'Mb' );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
237
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
238 #HTML Details
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
239 my $prefix_details_pages = $dir.$fastq_n[$child].'-'.$types_names[$grand_child];
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
240 details_pages ( $type_dir, $prefix_details_pages, \@fastq_n, $fastq_n[$child], $misTE, $dir, $Pcheck );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
241
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
242 $pm2->finish();
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
243 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
244 $pm2->wait_all_children;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
245
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
246 if ( $Pcheck eq 'true' )
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
247 {
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
248 my $ppp = $group_dir.'PPPartners/'; mkdir $ppp;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
249 print $report "ping_pong_partners $group_dir/piRNAs/TEs.sam $ppp\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
250 ping_pong_partners ( $group_dir.'TEs.fai', $group_dir.'piRNAs/piRNAs-TEs_sorted.bam', $ppp, $pi_min );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
251 my $ppp_page = $dir.$fastq_n[$child].'-piRNAs-PPP.html';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
252 ppp_page ( $group_dir, $ppp_page, \@fastq_n, $fastq_n[$child], $ppp, $dir );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
253 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
254
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
255 #HTML Main Webpage
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
256 my $index_page = $dir.$fastq_n[$child].'.html';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
257 main_page ( $gen_dir, $index_page, \@fastq_n, $fastq_n[$child], $ma, $ma_uni, $dir );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
258 copy ($index_page, $html_out) if $child == 0;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
259 #HTML Menu
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
260 my $menu_page = $dir.$fastq_n[$child].'-sub.html';
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
261 menu_page ( $group_dir, $menu_page, \@fastq_n, $fastq_n[$child], $min, $max, $si_min, $si_max, $pi_min, $pi_max, $dir );
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
262 unlink glob "$group_dir*.sam"; unlink glob "$group_dir*.fastq";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
263 $pm->finish(); # pass an exit code to finish
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
264 }
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
265 $pm->wait_all_children;
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
266 unlink glob $dir."dataset_*symlink.fa*";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
267 print $report "Job done!\n";
4bc00caa60b4 Uploaded
pierre.pouchin
parents: 36
diff changeset
268 close $report;