Hi all,

I am dealing with a FASTA file with multiple IDs (the same ID can be repeated) and the corresponding sequence. My idea was to create an Hash with as keys the IDs, in order to have a list without repetitions, and as values an array containing the sequence/sequences of the same ID. I need that each sequence is an element of the array. My script works as long as it stores the first and the second sequence in the array, from the third the data are no longer correct...

Hope to have been clear enough! I tried everything but nothing works..and the best would be not using any modules.

Thanks 


my $tot = 0;
my $current_ID = 1;
while ( my $line = <IN> ){
    
    chomp( $line );
    
    if ( $line =~ /^>/ ) { 
            $title= $line;
            $tot++;
    }else{
        
            $seq.= $line;
        shift @{$hash{$title}};
        push (@{$hash{$title}}, $seq);
            
       }
                    
    if ($tot > $current_ID) {        
           $current_ID ++;
                  push (@{$hash{$title}}, $seq);
               $seq = "";
     }                       

}

[download]

In reply to Saving different values for the same key by using Hash of Arrays by beginner27

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.