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

update concat fpkm on ensid

parent 75a4aa68
Branches
Tags
No related merge requests found
......@@ -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{$symbol} = {ensembl=>$ensembl,type=>$type};
$names{$ensembl} = {symbol=>$symbol,type=>$type};
}
my @files = @ARGV;
......@@ -32,7 +32,7 @@ foreach $file (@files) {
while (my $line = <IN>) {
chomp($line);
my ($ensid,$gene,$chr,$strand,$start,$end,$cov,$fpkm,$tmp) = split(/\t/,$line);
$cts{$gene}{$sample} = $fpkm;
$cts{$ensid}{$sample} = $fpkm;
}
close IN;
}
......@@ -57,7 +57,7 @@ foreach $gene (keys %cts) {
#next unless ($sortct[0] >= 1);
next if ($names{$gene}{type} eq 'rRNA');
next if ($exc{$gene});
print OUT join("\t",$names{$gene}{ensembl},$gene,
print OUT join("\t",$gene,$names{$gene}{symbol},
$names{$gene}{type},@line),"\n";
}
close OUT;
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