#! /usr/bin/perl my @prevrow=(); my @row=(); my @id=""; my $id2=""; my $chr=""; my $chr2=""; my $seleno=""; my $prevseleno=""; while () { chomp($_); @row=split(/\t/,$_); $seleno=$row[0]; if ($seleno eq $prevseleno) { push @id, $row[4]; $chr = substr ($row [1], 1); $chr2 = substr ($chr, 0, -1); } else { print "$prevseleno\t"; foreach $id (@id) { $id2 = substr ($id, 0, -1); print "~/gallopavodb/selenoprots/Uselenoprots/$id2.U.aa.fa ~/our_results/exonerate/aa/$id2.$chr2.rexonerate.aa ~/our_results/genewise/aa/$id2.$chr2.genewise.aa " ; } print "\n"; @id=""; $chr=""; $id[0]= $row [4]; } $prevseleno=$seleno; @prevrow=@row; }