medirect - a multiprocessed utility for retrieving records and parsing feature tables from ncbi
.. contents:: Table of Contents
Chris Rosenthal <crosenth@gmail.com>
_As a bioinformatician I build a lot of bacterial dna reference databases.
Part of my job is to gather sequences data where it is available from outside
sources. One of those sources is the NCBI nucleotide database. I designed this
package to help me gather data from NCBI quickly by utilizing multiple processors
to make multiple data requests from the NCBI database servers. The utilities
mefetch and ftract are designed to work like efetch and xtract and
can be slotted in along with other ncbi utilities and follow the same edirect
documentation <https://www.ncbi.nlm.nih.gov/books/NBK25501/>
__,
guidelines and requirements <https://www.ncbi.nlm.nih.gov/books/NBK25497/#_chapter2_Usage_Guidelines_and_Requiremen_>
and
usage policies <https://www.ncbi.nlm.nih.gov/home/about/policies.shtml>
. The
utilities have primarily been tested on the nucleotide database but should work
on any type of data available through the NCBI servers.
The mefetch utility is designed to be fast and can easily overwhelm the NCBI servers.
For this reason I highlight two points from the usage policy:
The ftract utility pattern matches features based on the three column table
structure described here: feature tables <http://www.ncbi.nlm.nih.gov/projects/Sequin/table.html>
It is designed to parse data and coordinates from feature tables which are
magnitudes smaller and faster to parse than the xml tables and xtract parser utility
available as part of the standard edirect package. The entire edirect package
is available here: ftp downloads <https://ftp.ncbi.nlm.nih.gov/entrez/entrezdirect/>
biopython <https://pypi.python.org/pypi/biopython>
_ >= 1.68retrying <https://pypi.python.org/pypi/retrying>
_ >= 1.3.3medirect can be installed in two ways:
For regular users::
% pip3 install medirect
For developers::
% pip3 install git://github.com/crosenth/medirect.git
% git clone git://github.com/crosenth/medirect.git % cd medirect % python3 setup.py install
The mefetch executable works exactly like edirect
efetch <https://www.ncbi.nlm.nih.gov/books/NBK179288/efetch>
_ with a an
additional multiprocessing argument -proc
and a few more features.
By allowing additional processes to download records The -proc
argument
allows a linear download speed increase downloading large datasets.
Here is an example downloading 255,303 Rhizobium sequence accessions using one processor::
% esearch -db nucleotide -query 'Rhizobium' | time mefetch -email user@ema.il -mode text -format acc -proc 1 > accessions.txt 0.53s user 0.11s system 0% cpu 12:43.11 total
Which is equivalent to ncbi efetch::
% esearch -db nucleotide -query 'Rhizobium' | time efetch -mode text -format acc > accessions.txt 0.53s user 0.11s system 0% cpu 12:47.54 total
Adding another processor -proc 2
::
% esearch -db nucleotide -query 'Rhizobium' | time mefetch -email user@ema.il -proc 2 -mode text -format acc > accessions.txt 0.46s user 0.08s system 0% cpu 5:17.51 total
And another -proc 3
(default)::
% esearch -db nucleotide -query 'Rhizobium' | time mefetch -email user@ema.il -proc 3 -mode text -format acc > accessions.txt 0.35s user 0.10s system 0% cpu 2:57.01 total
And -proc 4
(see usage policy)::
% esearch -db nucleotide -query 'Rhizobium' | time mefetch -email user@ema.il -proc 4 -mode text -format acc > accessions.txt 0.35s user 0.08s system 0% cpu 1:40.54 total
Results can be returned in the exact order they intended by the NCBI
server using the -in-order
argument. Otherwise, the order will be
determined by how fast ncbi returns results per process.
The -retmax
argument (or chunksize) determines the number of results
returned per -proc
. By default, it is set to the 10,000 max records per
documentation <https://www.ncbi.nlm.nih.gov/books/NBK25499/#_chapter4_EFetch_>
__.
Setting the -retmax
to higher than 10,000 will automatically be set back
down to 10,000.
By default the -id
reads stdin xml output from esearch
. The -id
argument can also take input as a comma delimited list of ids or text
file of ids. When coupled with the -csv
argument the input can be a csv
file with additional argument columns. This is useful for bulk downloads with
different positional arguments.
ftract allows csv output of different features from ncbi
feature tables <http://www.ncbi.nlm.nih.gov/projects/Sequin/table.html>
_.
The required -feature
argument is comma separated
feature_key:qualifier_key:qualifier_value
::
% mefetch -id KN150849 -db nucleotide -email user@ema.il -format ft | ftract --feature rrna:product:16s id,seq_start,seq_stop,strand KN150849.1,594136,595654,2 KN150849.1,807985,809503,2 KN150849.1,2227751,2229271,1
And pipe this back into mefetch to download these three regions in genbank format::
% mefetch -id KN150849 -db nucleotide -email user@ema.il -format ft | ftract --feature rrna:product:16s | mefetch -db nucleotide -email crosenth@uw.edu -csv -format gb
And finally combining all these concepts, return all the Burkholderia gladioli 16s rrna products in fasta format using the default -proc 3
like this::
% esearch -query 'Burkholderia gladioli AND sequence_from_type[Filter]' -db 'nucleotide' | mefetch -email user@ema.il -format ft | ftract --feature rrna:product:16s | mefetch -db nucleotide -email user@ema.il -csv -format fasta 0.24s user 0.05s system 1% cpu 18.596 total
Please use the Issue Tracker(s) available on Github or Bitbucket to report any bugs
or feature requests. For all other inquiries email Chris Rosenthal <crosenth@gmail.com>
_.
Copyright (c) 2016 Chris Rosenthal
Released under the GPLv3 <http://www.gnu.org/copyleft/gpl.html>
_ License