changeset 38:15fbf3bdb553 draft

Uploaded
author romaingred
date Tue, 28 Nov 2017 10:11:00 -0500 (2017-11-28)
parents af4ee65eaffb
children e36c4f4fe5fd
files bin/subgroups.pm
diffstat 1 files changed, 7 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/bin/subgroups.pm	Tue Nov 28 10:05:52 2017 -0500
+++ b/bin/subgroups.pm	Tue Nov 28 10:11:00 2017 -0500
@@ -13,7 +13,7 @@
 
 sub subgroups
 {
-	my ($fin, $dir, $mis, $mis_TE, $proc, $tRNAs, $rRNAs, $snRNAs, $miRNAs, $exons, $TE, $min_si, $max_si, $min_pi, $max_pi, $report ) = @_;
+	my ($fin, $dir, $mis, $mis_TE, $proc, $tRNAs, $rRNAs, $snRNAs, $miRNAs, $transcripts, $TE, $min_si, $max_si, $min_pi, $max_pi, $report ) = @_;
 	my (@files, $sum, $pie, $repar, %ismapped, %isjunk, %repartition, @junk_ref, @all_ref );
 
 	srand();
@@ -52,18 +52,18 @@
 	my $bo = $tmp[1];
 	unlink $sam, $sam.'_aln.err', $sam.'_samse.err';
 
-	my $sam_exons = $dir.'exons.sam'; 
-	my $reject_exons = $dir.'rejected_exons.fastq';
-	@tmp = get_hash_alignment($exons, $mis, 0, 1, 'NA', $reject_exons, $bonafide, $proc, 'exons', $sam_exons, $report, $dir.'exons.fai');
-	$repartition{'exons'} = $tmp[0];
+	my $sam_transcripts = $dir.'transcripts.sam'; 
+	my $reject_transcripts = $dir.'rejected_transcripts.fastq';
+	@tmp = get_hash_alignment($transcripts, $mis, 0, 1, 'NA', $reject_transcripts, $bonafide, $proc, 'transcripts', $sam_transcripts, $report, $dir.'transcripts.fai');
+	$repartition{'transcripts'} = $tmp[0];
 
 	
 	my $sam_TEs = $dir.'TEs.sam';
 	my $reject_TEs = $dir.'rejected.fastq';
-	@tmp = get_hash_alignment($TE, $mis_TE, 0, 1, 'NA', $reject_TEs, $reject_exons, $proc, 'TEs', $sam_TEs, $report, $dir.'TEs.fai' );
+	@tmp = get_hash_alignment($TE, $mis_TE, 0, 1, 'NA', $reject_TEs, $reject_transcripts, $proc, 'TEs', $sam_TEs, $report, $dir.'TEs.fai' );
 	$repartition{'TEs'} = $tmp[0] ; $repartition{'others'} = $tmp[1];
 	unlink $sam, $sam.'_aln.err', $sam.'_samse.err';
-	unlink $reject_exons;
+	unlink $reject_transcripts;
 	unlink $reject_rRNAs;
 	unlink $reject_miRNAs;
 	unlink $reject_tRNAs;