Hi Perlmonks,

My intention to sort by the third column in the output is creating me trouble with the code I wrote. This code is a follow up code on the output generated with another code (Title: Reverse complement) bluray's writeups. Now, my aim is to create unique identifiers for each of the 11 characters (tag) in the 1st column of the input file. The 2nd column is the frequency of each tag. I used this code to create a new first column. Each line of this column starts with >HWTI_frequency_randomnumber. Though, I was able to get the results with this code, I have trouble in sorting with the second column of my input file. 

#!usr/bin/perl
use strict;
use warnings;

my $input_file="file1.seq";
open my $FILE1, "<", $input_file or die "Cannot open $input_file .$!";
my $output_file="file2.csv";


open my $FILE2, ">", $output_file or die "Cannot open $output_file .$!
+";
my $line=<$FILE1>;
chomp $line;
my @columnheadings=split(/\t/, $line);
 unshift(@columnheadings, ("Header"));
my $heading=join("\t", @columnheadings);
 print $FILE2 "$heading\n";
my %tag;
while (my $line=<$FILE1>)
{
    chomp $line;
    $line=~s/\t/,/g;
    my @columns=split(/,/, $line);
    my $tags=$columns[0];
    
     $tag{$tags}=$line;
}
    foreach my $tags (sort keys %tag){
    my $header;
    my @columns=split(/,/,$tag{$tags});
        $tags=$columns[0];
       my $freq=$columns[1];
       my $range=500000;
    my $random_number=int(rand($range));
        $header=">HWTI_".$freq."_".$random_number;
    my $printline=$tag{$tags};
    $printline=$header.",".$printline;
        print $FILE2 "$printline\n";    
}

[download]

In addition to the sorting issue, I am also thinking about doing a BLAST for each of these tags (#nucleotide length of 11). I will appreciate any suggestions in this matter.

In reply to Sorting issue by bluray

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":


  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.