Skip to content
Snippets Groups Projects
Commit 91b0ea39 authored by Brandi Cantarel's avatar Brandi Cantarel
Browse files

full path in fpkm subset

parent b7202191
No related merge requests found
......@@ -8,14 +8,15 @@ use Getopt::Long;
my ($genelist,$fpkm_file) = ("") x 2;
GetOptions( 'g|genes=s' => \$genelist,
'f|fpkm=s' => \$fpkm_file);
'f|fpkm=s' => \$fpkm_file,
'i|refinfo=s'=> \$refinfo);
my $fpkm_out = $fpkm_file;
$fpkm_out =~ s/\.fpkm.+txt/\.fpkm.capture.txt/;
open OUT, ">$fpkm_out" or die $!;
#gene_info.human.txt file is required for gene alias information
open GENEINFO, "</project/shared/bicf_workflow_ref/human/gene_info.human.txt" or die $!;
open GENEINFO, "<$refinfo/gene_info.human.txt" or die $!;
my %geneinfo;
while(my $gline = <GENEINFO>){
chomp $gline;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment