GSP
Quick Navigator

Search Site

Unix VPS
A - Starter
B - Basic
C - Preferred
D - Commercial
MPS - Dedicated
Previous VPSs
* Sign Up! *

Support
Contact Us
Online Help
Handbooks
Domain Status
Man Pages

FAQ
Virtual Servers
Pricing
Billing
Technical

Network
Facilities
Connectivity
Topology Map

Miscellaneous
Server Agreement
Year 2038
Credits
 

USA Flag

 

 

Man Pages


Manual Reference Pages  -  BP_HIVQ.PL (1)

.ds Aq ’

NAME

bp_hivq.PL - an interactive command-line interface to Bio::DB::HIV and Bio::DB::Query::HIVQuery

CONTENTS

SYNOPSIS



 $ perl bp_hivq.PL
 hivq> query C[subtype] SI[phenotype]
 hivq> prerun
 80 sequences returned
 Query: C[subtype] SI[phenotype]
 hivq> outfile csi.fas
 hivq> run
 Download complete.
 hivq> outfile dsi.fas
 hivq> run D[subtype] SI[phenotype]
 Download complete.
 hivq> count
 25 sequences returned
 Query: D[subtype] SI[phenotype]
 hivq> exit
 $



DESCRIPTION

The BioPerl modules Bio::DB::HIV and Bio::DB::Query::HIVQuery together allow batch queries against the Los Alamos National Laboratories’ HIV Sequence Database using a simple query language. bp_hivq.PL provides both an example of the use of these modules, and a standalone interactive command-line interface to the LANL HIV DB. Simple commands allow the user to retrieve HIV sequences and annotations using the query language implemented in Bio::DB::Query::HIVQuery. Visit the man pages for those modules for more details.

USAGE

Run the script using perl bp_hivq.PL or, in Unix, ./bp_hivq.PL. You will see the



 hivq>



prompt. Type commands with queries to retrieve sequence and annotation data. See the SYNOPSIS for a sample session. Available commands are described below.

    TIPS

The LANL database is pretty complex and extensive. Use the find facility to explore the available database tables and fields. To identify aliases for a particular field, use find alias [fieldname]. For example, to find a short alias to the weirdly named field seq_sample.ssam_second_receptor, do



 hivq> find alias seq_sample.ssam_second_receptor



which returns



 coreceptor             second_receptor



Now, instead of the following query



 hivq> run C[subtype] CCR5[seq_sample.ssam_second_receptor]



you know you can do



 hivq> run C[subtype] CCR5[coreceptor]



Use the outfile command to set the file that receives the retrieved sequences. You can change the current output file simply by issuing a new outfile command during the session. The output file defaults to standard output.

Use the query command to validate a query without hitting the database. Use the prerun or count commands to get a count of sequence hits for a query without retrieving the data. Use run or do to perform a complete query, retrieving sequence data into the currently set output files.

To process bp_hivq.PL commands in batch, create a text file (bp_hivq.cmd, for example) containing desired commands one per line. Then execute the following from the shell:



 $ cat bp_hivq.cmd | perl bp_hivq.PL



COMMANDS

Here is a complete list of commands. Options in single brackets ([req_option]) are required; options in double brackets ([[opt_option]]) are optional.



 confirm            : Toggle interactive confirmation before
                      executing queries
 exit               : Quit script
 find               : Explore database schema
  find tables                 Display all database tables
  find fields                 Display all database fields (columns)
  find fields [table]         Display all fields in [table]
  find alias [field]          Display valid aliases for [field]
 help [[command]]   : Show command help
                      if [[command]] not specified, list all
                      available commands
 id                 : Display current session id
 outfile [filename] : Set file for collecting retrieved data
 ping               : Check if LANL DB is available
 prerun [[query]]   : Execute query but retreive hit count only
                      if [[query]] not specified, use current query
 query [query]      : Validate and set current query
 run [[query]]      : Execute query and retrieve data
                      if [[query]] not specified, use current query
 state              : Display current state of the script

 bye                : Alias for exit
 config             : Alias for state
 count              : Alias for prerun
 do                 : Alias for run
 out                : Alias for outfile
 quit               : Alias for exit



OPTIONS



 -v : verbose; turns on the internal debug() function



FEEDBACK

    Mailing Lists

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



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



    Reporting Bugs

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



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



AUTHOR - Mark A. Jensen

Mark A. Jensen <maj@fortinbras.us>
Search for    or go to Top of page |  Section 1 |  Main Index


perl v5.20.3 BP_HIVQ (1) 2016-04-05

Powered by GSP Visit the GSP FreeBSD Man Page Interface.
Output converted with manServer 1.07.