#!/usr/bin/perl

# Copyright (C) 1993-1995 J. Adachi,  All rights reserved.  Dec 05 1995

# Usage: info2mol [sequence_file]

chop($_ = <>);  # first line
($numotu,$numsite) = ( $_ =~ /(\d+)\s+(\d+)/ );
($datainfo) = ( $_ =~ /\d+\s+\d+\s+(\S.*)$/ );
#warn "\n$numotu OTUs $numsite sites.\n";

for ($i = 0; $i < $numotu; $i++) {
	while(<>) { 
		next if /^\s*$/;  # white line skip
		#next if /^\s*#/;  # comment line skip
		chop;
		($name) = ( $_ =~ /(\w\S*)/ );
		die "abnormal identifier(name): $_\n" unless $name;
		die "Identifier \"$name\" is double defined\n" if grep(/^$name$/,@names);
		($info) = ( $_ =~ /\w\S*\s+(\S.*)$/ );
		push(@names,$name);
		$infos{$name} = $info;
		last;
	}
}

while($leng < $numsite) { 
	for ($i = 0; $i < $numotu; $i++) {
		while(<>) { 
			next if /^\s*$/;  # white line skip
			#next if /^\s*#/;  # comment line skip
			chop;
			s/\s//g;
			$seq = $_;
			$seqs{$names[$i]} .= $seq;
			$leng = length($seqs{$names[$i]});
			#print "$names[$i]rbseq);

				print "# ",$de,"\n";
				print $name," ",$or1," ",$or2," # ",$ac," ", $so, "\n";
				print $len, "\n";
				$offset = 0;
				while ($line = substr($subseq, $offset, 60)) {
					print $line, "\n";
					$offset += 60;
				}
				$subseq = "";
				$cds = "";
				$de = "";
				print "\n";
			}
}
                                                                                                                                                                                                                                                                 ./usr/bin/inf2mol                                                                                   0000755 0000000 0000000 00000004022 10303532434 012561  0                                                                                                    ustar   root                            root                                                                                                                                                                                                                   #!/usr/bin/perl

