diff --git a/genect_rnaseq/concat_fpkm.pl b/genect_rnaseq/concat_fpkm.pl index c3be1e3af0af0a8f8b9f50777e92a416207ee235..18ee9997e7aa41eb92d1224e1ea0ff2db37a35b0 100755 --- a/genect_rnaseq/concat_fpkm.pl +++ b/genect_rnaseq/concat_fpkm.pl @@ -15,7 +15,7 @@ while (my $line = <SYM>) { chomp($line); my ($chrom,$start,$end,$ensembl,$symbol,$type) = split(/\t/,$line); $ensembl = (split(/\./,$ensembl))[0]; - $names{$ensembl} = {symbol=>$symbol,type=>$type}; + $names{$ensembl} = {chr=>$chrom,symbol=>$symbol,type=>$type}; } my @files = @ARGV; @@ -32,8 +32,13 @@ foreach $file (@files) { while (my $line = <IN>) { chomp($line); my ($ensid,$gene,$chr,$strand,$start,$end,$cov,$fpkm,$tmp) = split(/\t/,$line); - my ($ens,$version) = split(/\./,$ensid); - $cts{$ens}{$sample} = $fpkm; + my $ens = $ensid; + $ens =~ s/\.\d+//; + if ($chr eq $names{$ens}{chr}) { + $cts{$ens}{$sample} = $fpkm; + }else { + warn "unable to map to genenames\n"; + } } close IN; }