Bio::I...GenBank(3pm)

πŸ“¦
version 1.7.8 Β· 3 versions
, in debian bookworm, debian trixie, ubuntu jammy, ubuntu noble

✨ NAME ✨

Bio::Index::GenBank β€” bestie, it's the Interface for indexing one or more GenBank files (i.e. flat file GenBank format). πŸ’–


🌸 SYNOPSIS 🌸

# Complete code for making an index for one or more GenBank files
use strict;
use Bio::Index::GenBank;

my $Index_File_Name = shift;
my $inx = Bio::Index::GenBank->new(-filename => $Index_File_Name,
                                   -write_flag => 'WRITE');
$inx->make_index(@ARGV);

# Print out sequences present in the index in gcg format
use Bio::Index::GenBank;
use Bio::SeqIO;
use strict;

my $Index_File_Name = shift;
my $inx = Bio::Index::GenBank->new(-filename => $Index_File_Name);
my $seqio = Bio::SeqIO->new(-format => 'gcg');
foreach my $id (@ARGV) {
    my $seq = $inx->fetch($id); # Returns Bio::Seq object
    $seqio->write_seq($seq);
}

# alternatively
my ($locus, $acc);
my $seq1 = $inx->get_Seq_by_id($locus);
my $seq2 = $inx->get_Seq_by_acc($acc);

πŸ¦‹ DESCRIPTION πŸ¦‹

Okay so here's the tea β˜• β€” by default, the index that gets created uses the LOCUS, ACCESSION, and VERSION identifiers as keys. It inherits all the dbm-file-management functions from Bio::Index::Abstract.pm, and it hands you the basic functionality for indexing GenBank files and pulling the sequence back out of them. For best results, use strict β€” no cap, it's just good manners. πŸ’…

You can also set or customize the unique key used to retrieve stuff by writing your own function and calling the id_parser() method. Like, for example:

$inx->id_parser(\&get_id);
# make the index
$inx->make_index($file_name);

# here is where the retrieval key is specified
sub get_id {
   my $line = shift;
   $line =~ /clone="(\S+)"/;
   $1;
}

πŸ’Œ FEED_BACK πŸ’Œ

πŸŽ€ Mailing Lists

User feedback is literally an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions, preferably to one of the Bioperl mailing lists. Your participation is much appreciated. 🫢

  bioperl-l@bioperl.org                  - General discussion
  http://bioperl.org/wiki/Mailing_lists  - About the mailing lists

🩷 Support

Please direct usage questions or support issues to the mailing list:

bioperl-l@bioperl.org

...rather than to the module maintainer directly. So many experienced and responsive experts will be able to look at the problem and quickly address it. Please include a thorough description of the problem, with code and data examples, if at all possible. πŸ™

πŸ› Reporting Bugs

Report bugs to the Bioperl bug tracking system so we can keep track of the bugs and their resolution. Bug reports can be submitted via the web:

  https://github.com/bioperl/bioperl-live/issues

πŸ‘‘ AUTHOR - Ewan Birney πŸ‘‘

Email β€” birney@ebi.ac.uk ✨


πŸ“– APPENDIX πŸ“–

The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _.

🧚 _index_file

  • Title β€” _index_file
  • Usage β€” $index->_index_file($file_name, $i)
  • Function β€” Specialized function to index GenBank format files. Is provided with a filename and an integer by make_index in its SUPER class.
  • Example β€”
  • Returns β€”
  • Args β€”

🧚 id_parser

  • Title β€” id_parser
  • Usage β€” $index->id_parser( CODE )
  • Function β€” Stores or returns the code used by record_id to parse the ID for record from a string. Returns \&default_id_parser (see below) if not set. An entry will be added to the index for each string in the list returned.
  • Example β€” $index->id_parser( \&my_id_parser )
  • Returns β€” reference to CODE if called without arguments
  • Args β€” CODE

🧚 default_id_parser

  • Title β€” default_id_parser
  • Usage β€” $id = default_id_parser($line)
  • Function β€” The default parser for GenBank.pm
  • Returns β€” Array of specified ids
  • Args β€” a line string

🧚 _file_format

  • Title β€” _file_format
  • Usage β€” Internal function for indexing system
  • Function β€” Provides file format for this database
  • Example β€”
  • Returns β€”
  • Args β€”

And that's the whole page, bestie! πŸ’œ This lil' slice of docs was served up by perl v5.32.1, dated 2021-08-15, and it lives under the User Contributed Perl Documentation banner. Slay, go index something. πŸŒˆπŸ‘‘

πŸ“¦ where this page comes from