??? Your first paragraph makes no sense. You can't tell what snp.txt looks like if you print @output. The contents of @output is processed, changed by a split and a grep. Use an editor or 'less' to look at a file.

You might try out the following (this is the same program, just added a Dumper-line before the foreach:

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

use Data::Dumper;

use strict;
my$line;
my@fields;
my@output;
open (FILE1, 'snp.txt') or die "can't open the file: $!";
open (FILE2,'chr22.txt') or die "can't open the file: $!";
open (FD, '>test.txt') or die "can't open the file: $!";
my $position = tell(FILE2);
my %rs;
while ($line=<FILE2>) {
  my ($key)= $line=~/^(rs\d{5,})\b/;
  if (defined $key) {
  $rs{$key}= $position;
  }
  $position=tell(FILE2);
}

while (defined ($line= <FILE1>)) {
    my@fields= split (/\s+/ ,$line);
    push @output, grep /^rs\d{5,}\b/ ,@fields;
}

print Dumper(\@output,\%rs);

        foreach (@output) {
              if (exists $rs{$_}) {
                seek(FILE2,$rs{$_},0);
                my $line= <FILE2>;
    print FD $line;

    }
}
close FILE1;
close FILE2;
close FD;

[download]

With the data I posted, I get the following output: 

$VAR1 = [
          'rs34569384',
          'rs123456',
          'rs234567',
          'rs753444'
        ];
$VAR2 = {
          'rs123456' => 43,
          'rs234567' => 15
        };

[download]

As you can see @output (==$VAR1) contains a lot of rs-numbers. %rs (==$VAR2) contains some of the same rs-numbers and corresponding file positions. You should see the same if you use my data.

Now try it with your data. What is different? Are there rs-numbers that are in both @output and %rs? Do the file position numbers look correct?

In reply to Re^9: General program and related problems by jethro
in thread General program and related problems by micky744monk

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.