Skip to content

HTTPS clone URL

Subversion checkout URL

You can clone with HTTPS or Subversion.

Download ZIP
Ruby binding for the tabix file indexing routines within the samtools package http://samtools.sourceforge.net/
Ruby
tree: ab30617c46

Fetching latest commit…

Cannot retrieve the latest commit at this time

Failed to load latest commit information.
ext/tabix
lib
test
.document
.gitignore
.rvmrc
Gemfile
LICENSE.txt
README.rdoc
Rakefile
VERSION
bio-tabix.gemspec

README.rdoc

bio-tabix

Ruby binding for the tabix file indexing routines within the samtools package samtools.sourceforge.net/

Tabix provides utilities for indexing and subsequently querying regions of interest from large tab delimited files.

Files are indexed on three columns: [group, pos1, pos2] and must be position sorted

This gem was modeled on the bio-samtools gem: github.com/helios/bioruby-samtools

Installation

gem install bio-tabix

Usage

Compress the file

The following command will run the bgzip utility. You must supply input and output filenames. If you don't compress your data first, it will be compressed when the file is opened appending '.bgzf' to the filename. Bio::Tabix::TFile.compress(in_file,compressed_file)

Index the file

Build an index by supplying the required columns. If you don't create an index the open method will use the default options. The opts parameter take a hash of index build options - :s => sequence/group column [1] - :b => beginning range column [2] - :e => ending range column. Can equal :b. [3] - :meta_char => comment character [#] - :line_skip => number of initial lines to ignore [0] tabix_file = Bio::Tabix::TFile.build_index(compressed_file, {:s => group_col, :b => pos1_col, :e => pos2_col})

Open the file

Create a new TFile instance and open it.

	 tabix_file = Bio::Tabix::TFile.open(compressed_file)

Create a proc or lambda

This will be called with the value of each fileline my_func = lambda do |line| # convert text to array and print column 6 puts line.split(“t”)[6] end

Process a region

Choose a group and range for your function and process it tabix_file.process_region('FirstGroup', 1, 500, my_func)

Dependencies:

-FFI (github.com/ffi/ffi)

Contributing to bio-tabix

  • Check out the latest master to make sure the feature hasn't been implemented or the bug hasn't been fixed yet.

  • Check out the issue tracker to make sure someone already hasn't requested it and/or contributed it.

  • Fork the project.

  • Start a feature/bugfix branch.

  • Commit and push until you are happy with your contribution.

  • Make sure to add tests for it. This is important so I don't break it in a future version unintentionally.

  • Please try not to mess with the Rakefile, version, or history. If you want to have your own version, or is otherwise necessary, that is fine, but please isolate to its own commit so I can cherry-pick around it.

Copyright

Copyright © 2012 Nicholas A Thrower. See LICENSE.txt for further details.

Something went wrong with that request. Please try again.