Mercurial > repos > romaingred > pirna_pipeline

changeset 46:a03408665904 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author romaingred
date Fri, 01 Dec 2017 08:14:22 -0500
parents ffc66227835f
children 851c8bbb3214
files bin/align.pm
diffstat 1 files changed, 5 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
--- a/bin/align.pm	Fri Dec 01 08:06:35 2017 -0500
+++ b/bin/align.pm	Fri Dec 01 08:14:22 2017 -0500
@@ -46,8 +46,8 @@
 	my $funi = $col_prefix.'unique_mappers.fastq';
 	my $frej = $col_prefix.'unmapped.fastq';
 	
-  my $repartition = $prefix.'distribution.txt';
-	my $png_rep = $prefix.'distribution.png';
+  my $repartition = $out_prefix.'distribution.txt';
+	my $png_rep = $out_prefix.'distribution.png';
 	my ( %duplicates, %genome_hits) ;
 
 	#alignement to the first reference
@@ -77,9 +77,9 @@
 		histogram($size_mappedHashR, $png_rep, $total);
 
 
-		my $dup = $prefix.'dup_mapnum.txt';
-		my $dup_u = $prefix .'dup_unique.txt';
-		my $dup_r = $prefix .'dup_nonmapp.txt';
+		my $dup = $out_prefix.'dup_mapnum.txt';
+		my $dup_u = $out_prefix .'dup_unique.txt';
+		my $dup_r = $out_prefix .'dup_nonmapp.txt';
 		open(my $tab,">".$dup) || die "cannot open output txt file\n";
 		open(my $tab_r,">".$dup_r) || die "cannot open output txt file\n";
 		open(my $tab_u,">".$dup_u) || die "cannot open output txt file\n";