in reply to Re: Multiple Key Problem help
in thread Multiple Key Problem help
My updated code is:-Can't use string("72")as an array reference while "strict refs" in use + at test.pl line 21, <$fh> line 2
#!/usr/bin/perl use strict; use warnings; my $qfn1 = "wah.txt"; my $qfn2 = "taj.txt"; my %positions; { open(my $fh, '<', $qfn1) or die("Cannot open file \"$qfn1\": $!\n"); while (<$fh>) { chomp; my ($key, $pos) = split /\s+/; if (!$positions{$key}) { $positions{$key} = $pos; } else { my $ref = $positions{$key}; push @$ref,$pos; } } } my %sequences; { open(my $fh, '<', $qfn2) or die("Cannot open file \"$qfn2\": $!\n"); my $key; while (<$fh>) { if ( s/^>// ) { $key = ( split /\|/ )[1]; } else { chomp; $sequences{$key} .= $_; } } } for my $key ( sort keys %positions ) { my $ref = $positions{$key}; foreach my $value (@$ref) { if ( ! exists( $sequences{$key} )) { warn "KEY $key in File1 not found in File2\n"; next; } if ( length( $sequences{$key} ) < $positions{$key} ) { warn "KEY $key: File2 string length too short for File1 positi +on value\n"; next; } my $index = rindex( $sequences{$key}, "ATG", $positions{$key} ); if ( $index < 0 ) { warn sprintf( "KEY %s: No ATG in File2 string prior to positio +n %d\n", $key, $positions{$key} ); next; } $index += 3 while ( ($index + 3) < $positions{$key} ); print "$key $positions{$key} " . substr($sequences{$key}, $index, +3) . "\n"; } }
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^3: Multiple Key Problem help
by lostjimmy (Chaplain) on Dec 23, 2008 at 19:13 UTC | |
by Anonymous Monk on Dec 24, 2008 at 16:51 UTC | |
|