Mercurial > repos > estrain > amrformat
annotate amrformat.pl @ 2:ae0d90df065f
Uploaded
author | estrain |
---|---|
date | Wed, 29 Jan 2020 16:33:31 -0500 |
parents | 514bd7256422 |
children | abe5635dd846 |
rev | line source |
---|---|
estrain@0 | 1 #!/usr/bin/perl |
estrain@0 | 2 |
estrain@0 | 3 open IN,$ARGV[0]; |
estrain@0 | 4 |
estrain@0 | 5 open OUT,">$ARGV[1]"; |
estrain@0 | 6 |
estrain@0 | 7 my %iso_hash; |
estrain@0 | 8 |
estrain@0 | 9 ## Skip header |
estrain@0 | 10 <IN>; |
estrain@0 | 11 |
estrain@0 | 12 while(<IN>) { |
estrain@0 | 13 chomp(); |
estrain@0 | 14 @vals=split(/\t/,$_); |
estrain@0 | 15 @els=split(/\_/,$vals[0]); |
estrain@0 | 16 if(exists $iso_hash{$els[0]}){ |
estrain@0 | 17 $iso_hash{$els[0]}="$iso_hash{$els[0]} $vals[1]"; |
estrain@0 | 18 } else { |
estrain@0 | 19 $iso_hash{$els[0]}=$vals[1]; |
estrain@0 | 20 } |
estrain@0 | 21 } |
estrain@0 | 22 |
estrain@0 | 23 |
estrain@0 | 24 print OUT "Source.Name\tGene symbol\n"; |
estrain@0 | 25 foreach(sort keys %iso_hash){ |
estrain@0 | 26 print OUT $_,"\t",$iso_hash{$_},"\n"; |
estrain@0 | 27 } |
estrain@0 | 28 |