.\" Automatically generated by Pod::Man v1.37, Pod::Parser v1.32 .\" .\" Standard preamble: .\" ======================================================================== .de Sh \" Subsection heading .br .if t .Sp .ne 5 .PP \fB\\$1\fR .PP .. .de Sp \" Vertical space (when we can't use .PP) .if t .sp .5v .if n .sp .. .de Vb \" Begin verbatim text .ft CW .nf .ne \\$1 .. .de Ve \" End verbatim text .ft R .fi .. .\" Set up some character translations and predefined strings. \*(-- will .\" give an unbreakable dash, \*(PI will give pi, \*(L" will give a left .\" double quote, and \*(R" will give a right double quote. | will give a .\" real vertical bar. \*(C+ will give a nicer C++. Capital omega is used to .\" do unbreakable dashes and therefore won't be available. \*(C` and \*(C' .\" expand to `' in nroff, nothing in troff, for use with C<>. .tr \(*W-|\(bv\*(Tr .ds C+ C\v'-.1v'\h'-1p'\s-2+\h'-1p'+\s0\v'.1v'\h'-1p' .ie n \{\ . ds -- \(*W- . ds PI pi . if (\n(.H=4u)&(1m=24u) .ds -- \(*W\h'-12u'\(*W\h'-12u'-\" diablo 10 pitch . if (\n(.H=4u)&(1m=20u) .ds -- \(*W\h'-12u'\(*W\h'-8u'-\" diablo 12 pitch . ds L" "" . ds R" "" . ds C` "" . ds C' "" 'br\} .el\{\ . ds -- \|\(em\| . ds PI \(*p . ds L" `` . ds R" '' 'br\} .\" .\" If the F register is turned on, we'll generate index entries on stderr for .\" titles (.TH), headers (.SH), subsections (.Sh), items (.Ip), and index .\" entries marked with X<> in POD. Of course, you'll have to process the .\" output yourself in some meaningful fashion. .if \nF \{\ . de IX . tm Index:\\$1\t\\n%\t"\\$2" .. . nr % 0 . rr F .\} .\" .\" For nroff, turn off justification. Always turn off hyphenation; it makes .\" way too many mistakes in technical documents. .hy 0 .if n .na .\" .\" Accent mark definitions (@(#)ms.acc 1.5 88/02/08 SMI; from UCB 4.2). .\" Fear. Run. Save yourself. No user-serviceable parts. . \" fudge factors for nroff and troff .if n \{\ . ds #H 0 . ds #V .8m . ds #F .3m . ds #[ \f1 . ds #] \fP .\} .if t \{\ . ds #H ((1u-(\\\\n(.fu%2u))*.13m) . ds #V .6m . ds #F 0 . ds #[ \& . ds #] \& .\} . \" simple accents for nroff and troff .if n \{\ . ds ' \& . ds ` \& . ds ^ \& . ds , \& . ds ~ ~ . ds / .\} .if t \{\ . ds ' \\k:\h'-(\\n(.wu*8/10-\*(#H)'\'\h"|\\n:u" . ds ` \\k:\h'-(\\n(.wu*8/10-\*(#H)'\`\h'|\\n:u' . ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'^\h'|\\n:u' . ds , \\k:\h'-(\\n(.wu*8/10)',\h'|\\n:u' . ds ~ \\k:\h'-(\\n(.wu-\*(#H-.1m)'~\h'|\\n:u' . ds / \\k:\h'-(\\n(.wu*8/10-\*(#H)'\z\(sl\h'|\\n:u' .\} . \" troff and (daisy-wheel) nroff accents .ds : \\k:\h'-(\\n(.wu*8/10-\*(#H+.1m+\*(#F)'\v'-\*(#V'\z.\h'.2m+\*(#F'.\h'|\\n:u'\v'\*(#V' .ds 8 \h'\*(#H'\(*b\h'-\*(#H' .ds o \\k:\h'-(\\n(.wu+\w'\(de'u-\*(#H)/2u'\v'-.3n'\*(#[\z\(de\v'.3n'\h'|\\n:u'\*(#] .ds d- \h'\*(#H'\(pd\h'-\w'~'u'\v'-.25m'\f2\(hy\fP\v'.25m'\h'-\*(#H' .ds D- D\\k:\h'-\w'D'u'\v'-.11m'\z\(hy\v'.11m'\h'|\\n:u' .ds th \*(#[\v'.3m'\s+1I\s-1\v'-.3m'\h'-(\w'I'u*2/3)'\s-1o\s+1\*(#] .ds Th \*(#[\s+2I\s-2\h'-\w'I'u*3/5'\v'-.3m'o\v'.3m'\*(#] .ds ae a\h'-(\w'a'u*4/10)'e .ds Ae A\h'-(\w'A'u*4/10)'E . \" corrections for vroff .if v .ds ~ \\k:\h'-(\\n(.wu*9/10-\*(#H)'\s-2\u~\d\s+2\h'|\\n:u' .if v .ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'\v'-.4m'^\v'.4m'\h'|\\n:u' . \" for low resolution devices (crt and lpr) .if \n(.H>23 .if \n(.V>19 \ \{\ . ds : e . ds 8 ss . ds o a . ds d- d\h'-1'\(ga . ds D- D\h'-1'\(hy . ds th \o'bp' . ds Th \o'LP' . ds ae ae . ds Ae AE .\} .rm #[ #] #H #V #F C .\" ======================================================================== .\" .IX Title "Bio::ASN1::Sequence 3" .TH Bio::ASN1::Sequence 3 "2005-05-04" "perl v5.8.8" "User Contributed Perl Documentation" .SH "NAME" Bio::ASN1::Sequence \- Regular expression\-based Perl Parser for ASN.1\-formatted NCBI Sequences. .SH "SYNOPSIS" .IX Header "SYNOPSIS" .Vb 1 \& use Bio::ASN1::Sequence; .Ve .PP .Vb 5 \& my $parser = Bio::ASN1::Sequence->new('file' => "downloaded.asn1"); \& while(my $result = $parser->next_seq) \& { \& # extract data from $result, or Dumpvalue->new->dumpValue($result); \& } .Ve .PP .Vb 4 \& # a new way to get the $result data hash for a particular sequence id: \& use Bio::ASN1::Sequence::Indexer; \& my $inx = Bio::ASN1::Sequence::Indexer->new(-filename => 'seq.idx'); \& my $seq = $inx->fetch_hash('AF093062'); .Ve .PP .Vb 2 \& # for creation of .idx index files please refer to \& # Bio::ASN1::Sequence::Indexer perldoc .Ve .SH "PREREQUISITE" .IX Header "PREREQUISITE" None. .SH "INSTALLATION" .IX Header "INSTALLATION" Bio::ASN1::Sequence is part of the Bio::ASN1::EntrezGene package. Bio::ASN1::EntrezGene package can be installed & tested as follows: .PP .Vb 4 \& perl Makefile.PL \& make \& make test \& make install .Ve .SH "DESCRIPTION" .IX Header "DESCRIPTION" Bio::ASN1::Sequence is a regular expression-based Perl Parser for \s-1ASN\s0.1\-formatted \&\s-1NCBI\s0 sequences. It parses an \s-1ASN\s0.1\-formatted sequence record and returns a data structure that contains all data items from the sequence record. .PP The parser will report error & line number if input data does not conform to the \&\s-1NCBI\s0 Sequence annotation file format. .PP The sequence parser is basically a modified version of the high-performance Bio::ASN1::EntrezGene parser. However, I created a standalone module for sequence since it is more efficient to keep Sequence-specific code out of EntrezGene.pm. .PP In fact it is possible to provide reading of all \s-1NCBI\s0's \s-1ASN\s0.1\-formatted files through simple variations of the Entrez Gene parser (I need more investigation to be sure, but at least the sequence parser works well). .PP Since demand for parsing \s-1NCBI\s0 \s-1ASN\s0.1\-formatted sequences is much lower than EntrezGene, this module is more like a beta version that works on the examples I checked, but I did not check all available records or data definitions. The error-reporting function of this module has to be useful sometimes. :) .SH "SEE ALSO" .IX Header "SEE ALSO" The parse_sequence_example.pl script included in this package (please see the Bio\-ASN1\-EntrezGene\-x.xx/examples directory) shows the usage. .PP Please check out perldoc for Bio::ASN1::EntrezGene for more info. .SH "AUTHOR" .IX Header "AUTHOR" Dr. Mingyi Liu .SH "COPYRIGHT" .IX Header "COPYRIGHT" The Bio::ASN1::EntrezGene module and its related modules and scripts are copyright (c) 2005 Mingyi Liu, \s-1GPC\s0 Biotech \s-1AG\s0 and Altana Research Institute. All rights reserved. I created these modules when working on a collaboration project between these two companies. Therefore a special thanks for the two companies to allow the release of the code into public domain. .PP You may use and distribute them under the terms of the Perl itself or \&\s-1GPL\s0 (). .SH "CITATION" .IX Header "CITATION" Liu, M and Grigoriev, A (2005) \*(L"Fast Parsers for Entrez Gene\*(R" Bioinformatics. In press .SH "OPERATION SYSTEMS SUPPORTED" .IX Header "OPERATION SYSTEMS SUPPORTED" Any \s-1OS\s0 that Perl runs on. .SH "METHODS" .IX Header "METHODS" .Sh "new" .IX Subsection "new" .Vb 12 \& Parameters: maxerrstr => 20 (optional) - maximum number of characters after \& offending element, used by error reporting, default is 20 \& file or -file => $filename (optional) - name of the file to be \& parsed. call next_seq to parse! \& fh or -fh => $filehandle (optional) - handle of the file to be \& parsed. \& Example: my $parser = Bio::ASN1::Sequence->new(); \& Function: Instantiate a parser object \& Returns: Object reference \& Notes: Setting file or fh will reset line numbers etc. that are used \& for error reporting purposes, and seeking on file handle would \& mess up linenumbers! .Ve .Sh "maxerrstr" .IX Subsection "maxerrstr" .Vb 6 \& Parameters: $maxerrstr (optional) - maximum number of characters after \& offending element, used by error reporting, default is 20 \& Example: $parser->maxerrstr(20); \& Function: get/set maxerrstr. \& Returns: maxerrstr. \& Notes: .Ve .Sh "parse" .IX Subsection "parse" .Vb 18 \& Parameters: $string that contains Sequence record, \& $trimopt (optional) that specifies how the data structure \& returned should be trimmed. 2 is recommended and \& default \& $noreset (optional) that species that line number should not \& be reset \& DEPRECATED as external function!!! Do not call this function \& directly! Call next_seq() instead \& Example: my $value = $parser->parse($text); # DEPRECATED as \& # external function!!! Do not call this function \& # directly! Call next_seq() instead \& Function: Takes in a string representing Sequence record, parses \& the record and returns a data structure. \& Returns: A data structure containing all data items from the sequence \& record. \& Notes: DEPRECATED as external function!!! Do not call this function \& directly! Call next_seq() instead \& $string should not contain 'Seq-entry ::= set' at beginning! .Ve .Sh "input_file" .IX Subsection "input_file" .Vb 9 \& Parameters: $filename for file that contains Sequence record(s) \& Example: $parser->input_file($filename); \& Function: Takes in name of a file containing Sequence records. \& opens the file and stores file handle \& Returns: none. \& Notes: Attemps to open file larger than 2 GB even on Perl that \& does not support 2 GB file (accomplished by calling \& "cat" and piping output. On OS that does not have "cat" \& error message will be displayed) .Ve .Sh "next_seq" .IX Subsection "next_seq" .Vb 17 \& Parameters: $trimopt (optional) that specifies how the data structure \& returned should be trimmed. option 2 is recommended and \& default \& Example: my $value = $parser->next_seq(); \& Function: Use the file handle generated by input_file, parses the next \& the record and returns a data structure. \& Returns: A data structure containing all data items from the sequence \& record. \& Notes: Must pass in a filename through new() or input_file() first! \& For details on how to use the $trimopt data trimming option \& please see comment for the trimdata method. An option \& of 2 is recommended and default \& The acceptable values for $trimopt include: \& 1 - trim as much as possibile \& 2 (or 0, undef) - trim to an easy-to-use structure \& 3 - no trimming (in version 1.06, prior to version \& 1.06, 0 or undef means no trimming) .Ve .Sh "trimdata" .IX Subsection "trimdata" .Vb 17 \& Parameters: $hashref or $arrayref \& $trimflag (optional, see Notes) \& Example: trimdata($datahash); # using the default flag \& Function: recursively process all attributes of a hash/array \& hybrid and get rid of any arrayref that points to \& one-element arrays (trims data structure) depending on \& the optional flag. \& Returns: none - trimming happenes in-place \& Notes: This function is useful to compact a data structure produced by \& Bio::ASN1::Sequence::parse. \& The acceptable values for $trimopt include: \& 1 - trim as much as possibile \& 2 (or 0, undef) - trim to an easy-to-use structure \& 3 - no trimming (in version 1.06, prior to version \& 1.06, 0 or undef means no trimming) \& This function is duplicate to EntrezGene.pm's and code should \& be compressed in the future (using util module & subclass). .Ve .Sh "fh" .IX Subsection "fh" .Vb 7 \& Parameters: $filehandle (optional) \& Example: trimdata($datahash); # using the default flag \& Function: getter/setter for file handle \& Returns: file handle for current file being parsed. \& Notes: Use with care! \& Line number report would not be corresponding to file's line \& number if seek operation is performed on the file handle! .Ve .Sh "rawdata" .IX Subsection "rawdata" .Vb 8 \& Parameters: none \& Example: my $data = $parser->rawdata(); \& Function: Get the sequence data file that was just parsed \& Returns: a string containing the ASN1-formatted sequence record \& Notes: Must first parse a record then call this function! \& Could be useful in interpreting line number value in error \& report (if user did a seek on file handle right before parsing \& call) .Ve