# HG changeset patch # User dereeper # Date 1546960489 18000 # Node ID bba379e32116b2bff9f3bacb29ff019c1b4630ba # Parent 36e5445b7807cf310afbce0be75e1f7b34745342 Uploaded diff -r 36e5445b7807 -r bba379e32116 complete_names.pl --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/complete_names.pl Tue Jan 08 10:14:49 2019 -0500 @@ -0,0 +1,35 @@ +use strict; + +my $fasta = $ARGV[0]; +my $phylip = $ARGV[1]; +my $out = $ARGV[2]; + +my @ids; +open(F,$fasta) or die "Can not open file $fasta"; +while(){ + if (/>(.*)/){ + my $id = $1; + push(@ids,$id); + } +} +close(F); + +open(O,">$out") or die "Can not open and write into file $out"; +open(P,$phylip) or die "Can not open file $phylip"; +my $numline = 0; +while(

){ + my $line = $_; + $numline++; + if (/^([\w\-]+)\s+/ && $numline > 1){ + my $reported_id = $1; + print "ok"; + my $id = $ids[$numline-2]; + $line =~s/$reported_id/$id/g; + print O $line; + } + else{ + print O $line; + } +} +close(P); +close(O);