diff complete_names.pl @ 2:bba379e32116 draft

Uploaded
author dereeper
date Tue, 08 Jan 2019 10:14:49 -0500
parents
children 637ed4e3aa1d
line wrap: on
line diff
--- /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(<F>){
+	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(<P>){
+	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);