ptt.txt
Location Strand Length PID Gene Synonym Code COG + Product 1..1524 + 507 15607143 dnaA Rv0001 - COG0 +593L chromosomal replication initiator protein DnaA 2052..3260 + 402 15607144 dnaN Rv0002 - COG0592L + DNA polymerase III subunit beta 3280..4437 + 385 15607145 recF Rv0003 - COG1195L + DNA replication/repair protein RecF
####
mapped.txt
gi|50953765|ref|NC_002755.2|-587160 83 gi|448814763|ref|NC_00096 +2.3| 1329307
####
$file="mapped.txt";
open (FH1, "$file");
open (FA, ">>allassignment.txt");
while ($seq1=)
{
@seq1=split (/\t/, $seq1);
$pos1=@seq1[3];
open (FH, "ptt.txt");
while ($seq=)
{
@seq=split (/\t/, $seq);
$pos=@seq[0];
$pos8=@seq[8];
chomp $pos8;
@pos=split(/\.\./, $pos);
$pos_start=@pos[0];
$pos_end=@pos[1];
for ($i=$pos_start, $i<=$pos_end, $i++)
{if ($i=~m/^$pos1$/)
{$position1=@seq1[0]; $position8=@seq1[8]; $position9=@seq1[9];
print FA $i."\t".$pos1."\t".$pos."\t".$pos8."\t".$pos1."\t".$position1."\t".$position8."\t".$position9."\n" ;
last;}
}
}}}