Next Trivial Utilities 57

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();
        }


Next Menu Copyright © 2012 M. J. Dominus