Re^4: fasta hash

I have changed the code, but now my problem seems to be that I can only access the last line of the output outside the loops. I wonder if there's a way to store the variables inside the loop to be accessible outside. The code looks like this now.

if(@ARGV < 3){
die "Not enough arguments\n";
}
$sequence="";
$fastaID;
open(FILE1,"$ARGV[0]") or die "No fasta file provided in command line:
+ $!\n";
while ($line=<FILE1>){
chomp($line);
if ($line=~/^\s*$/){
next;
}elsif ($line=~/^.*$/){
$fastaID=$line;
$fastahash{$fastaID}=1;
}
}
open(FILE2,"$ARGV[2]") or die "No fasta file provided in command line:
+ $!\n";
while($line2=<FILE2>){
chomp($line2);
if ($line2=~/^>/){
@data=split(" ",$line2);

$fasta=$data[1];

$sequence="";
}else{
$sequence.=$line2;
}
}
if (exists $fastahash{$fasta}){
print "$fastaID\t $sequence\n";
}


exit;
[download]

And the output, which is just the last key value in the fastahash is

2056360013       Musacgagchagshgashcgahcgacacsasasasacsacsasasacacaasc
+assacsaascascascascac
[download]

Comment on Re^4: fasta hash Select or Download Code

Replies are listed 'Best First'.
Re^5: fasta hash by moritz (Cardinal) on Aug 26, 2011 at 16:38 UTC
Sadly you ignored half of my advice; it wasn't given without reason. Please reconsider your position. You `print` can only be executed once because it's not inside a loop (should be inside the while loop). Proper indention would have made that obvious, and `use strict` would have told you that `$fasta` isn't in scope in the `if (exists $fastahash{$fasta})` (if you would have declared it properly inside the while loop). Perl 6 - second systems done right	[reply] [d/l] [select]
Re^6: fasta hash by morio56 (Initiate) on Aug 26, 2011 at 19:10 UTC
Thanks for your help, but non of this things work even if i follow your instructions to the letter. Putting the if loop within the while loop only prints the ids and no sequences. Putting it in the big while loop prints several lines repeatedly	[reply]
Re^7: fasta hash by moritz (Cardinal) on Aug 26, 2011 at 19:40 UTC
So, what's your current code? And please indent it properly, and `use strict; use warnings;`. Otherwise I'm not going to comment on it. Perl 6 - second systems done right	[reply] [d/l]
Re^7: fasta hash by Marshall (Canon) on Aug 28, 2011 at 15:32 UTC
Instead of processing the second file line by line, it is possible to use the '>' character as an end-of-line input separator. This simplifies processing of the file. #!/usr/bin/perl -w use strict; open (IDS, '<', "fastaIds") or #with your data 2056360012 2056360013 die "cannot open fastaIds $!\n"; #-------------# my %ids; while (<IDS>) #process first file with only ID's { chomp; $ids{$_}=1; } #-------------# $/='>'; #input record separator is now '>' while (<DATA>) { chomp; # now works on '>' not \n next if /^\s$/; # first record will be blank my ($id) = /^\s(\d+)/; print ">$_" if $ids{$id}; #print all lines for this id } #-------------# =prints > 2056360012 1047627436237 yyyacgagchagshgashcgahcgac acsasasasacsacsasasacaca ascassacsaascascascascac > 2056360013 1047627436238 xxxxcgagchagshgashcgahcgac acsasasasacsacsasasacaca ascassacsaascascascascac =cut __DATA__ > 2056360012 1047627436237 yyyacgagchagshgashcgahcgac acsasasasacsacsasasacaca ascassacsaascascascascac > 2056360013 1047627436238 xxxxcgagchagshgashcgahcgac acsasasasacsacsasasacaca ascassacsaascascascascac [download]	[reply] [d/l]