mark: seq_read

        sub seq_read {
          my ($seqname) = @_;

          for my $n (0 .. $#SEQUENCE) {
            for ($SEQUENCE[$n]) {
              if (/^\Q$seqname\E:\s+(.*)/) {
                return Set::IntSpan->new($1);
              }
            }
          }
          return Set::IntSpan->new();
        }

• There's hidden magic here

• @SEQUENCE is tied to the database file

• Scanning the array reads the file

Next

Menu

Copyright © 2005 M. J. Dominus