Mercurial > repos > estrain > amrformat
annotate amrformat.pl @ 7:abe5635dd846 tip
Uploaded
| author | estrain |
|---|---|
| date | Thu, 30 Jan 2020 08:42:23 -0500 |
| parents | 514bd7256422 |
| children |
| 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@7 | 16 $iso_hash{$els[0]}{$vals[1]}=1; |
| estrain@0 | 17 } |
| estrain@0 | 18 |
| estrain@0 | 19 |
| estrain@0 | 20 print OUT "Source.Name\tGene symbol\n"; |
| estrain@0 | 21 foreach(sort keys %iso_hash){ |
| estrain@7 | 22 print OUT $_; |
| estrain@7 | 23 foreach my $gene (sort keys %{$iso_hash{$_}}){ |
| estrain@7 | 24 print OUT " ",$gene; |
| estrain@7 | 25 } |
| estrain@7 | 26 print OUT "\n"; |
| estrain@0 | 27 } |
| estrain@0 | 28 |