biosyntax 0.1.0 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 3d25bb3d93262912c86faa5bc27143a5d4f1f513aa66f728543d07fec0bc5221
-  data.tar.gz: bc548ee7509ea154df6126bacf794233d8c2d523b53afe0dafe2ddb195e87e48
+  metadata.gz: '08a0ce98e7f92f7c08a03c56ebcd37fee25ca4bfd8717bc6b6a0b4873d6ad609'
+  data.tar.gz: 7fb63a24759d4f1a095b5b33a30437c857f73aad4e6f5a0cdbd45ed75ccb1590
 SHA512:
-  metadata.gz: 5ab4fb91d4c75df89b737fcd1354b16637c1c631f8e87950cd1c15532041d6f9dc1e382ad0b90f2efb6ae34cf9013f790808797b31d33fd4d23ad246417ecc5c
-  data.tar.gz: 4cbd7875d8ca86d40c9d21ddc00c662e8bc5319cb02862ae55c10fe5e77b005a794793ee7763bf515c9c7a8aa2ec692a7026aebe5d19ee503a27887f95bc4552
+  metadata.gz: 2862c64ea68a56625424b382e4069e02b775a1085652a29e6f60da9d52d05415e30a461beeb86789a815a7a22a2a8598d17744da92e229404da0dc9f130fc06c
+  data.tar.gz: 61963733732c0acf300ef3edb168de21246842001112eb49dfa71c4ef7ff08a780cede0d36c64488716a767b97f1fdd5d3fe38d08aa41c5fe9c8386cb31b9199

data/README.md

@@ -1,7 +1,10 @@
 # ruby-biosyntax
 
 [![CI](https://github.com/kojix2/ruby-biosyntax/actions/workflows/ci.yml/badge.svg)](https://github.com/kojix2/ruby-biosyntax/actions/workflows/ci.yml)
+[![Gem Version](https://badge.fury.io/rb/biosyntax.svg)](https://badge.fury.io/rb/biosyntax)
 [![Lines of Code](https://img.shields.io/endpoint?url=https%3A%2F%2Ftokei.kojix2.net%2Fbadge%2Fgithub%2Fkojix2%2Fruby-biosyntax%2Flines)](https://tokei.kojix2.net/github/kojix2/ruby-biosyntax)
+[![DOI](https://zenodo.org/badge/1269890086.svg)](https://doi.org/10.5281/zenodo.20698478)
+
 
 
 :dna: [bioSyntax](https://github.com/bioSyntax/bioSyntax) - Syntax highlighting for biological data formats - for Ruby.
@@ -98,20 +101,18 @@ BioSyntax.guess_format("a.vcf.gz") # :vcf
 The metadata is generated from `libbiosyntax` at load time. The Ruby side does
 not maintain a separate hand-written table of formats or kinds.
 
-## Examples
+## Command line
 
-This gem does not install a CLI. See `examples/` for small scripts:
+Installing the gem also installs `biocat`:
 
 ```sh
-ruby examples/bcat.rb sample.vcf
-ruby examples/bcat.rb -l fastq reads.fastq
-ruby examples/bcat.rb -l
-ruby examples/inspect_spans.rb sample.vcf
+biocat sample.vcf
+biocat --format fastq reads.fastq
+biocat -l
 ```
 
-`bcat.rb` guesses the format from the file name when possible. Use `-l` /
-`--language` to pass a format explicitly. Calling `-l` without an argument
-prints the supported format names.
+`.gz`/`.bgz` are decompressed automatically. BAM/CRAM/BCF require optional
+`ruby-htslib` (`gem install htslib`).
 
 ## Development tasks
 

data/exe/biocat

@@ -0,0 +1,131 @@
+#!/usr/bin/env ruby
+# frozen_string_literal: true
+
+require 'optparse'
+require 'zlib'
+require 'biosyntax'
+
+options = { format: nil, list_formats: false }
+
+BINARY_FORMATS = {
+  '.bam' => :sam,
+  '.cram' => :sam,
+  '.bcf' => :vcf
+}.freeze
+
+def binary_format_for(path, requested_format)
+  unless path == '-'
+    ext = File.extname(path).downcase
+    return BINARY_FORMATS[ext] if BINARY_FORMATS.key?(ext)
+  end
+
+  case requested_format&.to_s&.downcase
+  when 'bam', 'cram' then :sam
+  when 'bcf' then :vcf
+  end
+end
+
+def require_htslib!
+  require 'htslib'
+rescue LoadError => e
+  warn "[biocat] ruby-htslib is required to read BAM/CRAM/BCF files: #{e.message}"
+  warn '[biocat] install it with `gem install htslib`'
+  exit 1
+end
+
+def emit_colored(text, highlighter)
+  text.each_line do |line|
+    print(highlighter ? highlighter.colorize(line) : line)
+  end
+end
+
+def emit_record(record, highlighter)
+  line = record.to_s
+  line = "#{line}\n" unless line.end_with?("\n")
+  emit_colored(line, highlighter)
+end
+
+def emit_htslib_file(path, format, highlighter)
+  require_htslib!
+
+  case format
+  when :sam
+    HTS::Bam.open(path) do |bam|
+      emit_colored(bam.header.to_s, highlighter)
+      bam.each { |record| emit_record(record, highlighter) }
+    end
+  when :vcf
+    HTS::Bcf.open(path) do |bcf|
+      emit_colored(bcf.header.to_s, highlighter)
+      bcf.each { |record| emit_record(record, highlighter) }
+    end
+  end
+end
+
+parser = OptionParser.new do |opts|
+  opts.banner = 'usage: biocat [options] [FILE ...]'
+
+  opts.on('-f', '--format FORMAT', 'Highlight as FORMAT') do |format|
+    options[:format] = format
+  end
+
+  opts.on('-l', '--list-formats', 'Print supported format names') do
+    options[:list_formats] = true
+  end
+
+  opts.on('-h', '--help', 'Print this help') do
+    puts opts
+    exit
+  end
+end
+
+begin
+  parser.parse!
+
+  if options[:list_formats]
+    puts BioSyntax::FORMAT_NAMES.join("\n")
+    exit
+  end
+
+  paths = ARGV.empty? ? ['-'] : ARGV
+
+  paths.each do |path|
+    binary_format = binary_format_for(path, options[:format])
+    format = binary_format || options[:format]
+    format ||= BioSyntax.guess_format(path) unless path == '-'
+
+    highlighter = format ? BioSyntax[format] : nil
+    if binary_format
+      emit_htslib_file(path, binary_format, highlighter)
+      next
+    end
+
+    input =
+      if path == '-'
+        $stdin
+      elsif path.downcase.end_with?('.gz', '.bgz')
+        Zlib::GzipReader.open(path)
+      else
+        File.open(path, 'rb')
+      end
+
+    begin
+      emit_colored(input, highlighter)
+    ensure
+      input.close unless path == '-'
+    end
+  end
+rescue OptionParser::ParseError, BioSyntax::Error => e
+  warn "[biocat] #{e.message}"
+  warn parser
+  exit 2
+rescue SystemCallError => e
+  warn "[biocat] #{e.message}"
+  exit 1
+rescue Zlib::GzipFile::Error => e
+  warn "[biocat] #{e.message}"
+  exit 1
+rescue StandardError => e
+  warn "[biocat] #{e.message}"
+  exit 1
+end

data/lib/biosyntax/version.rb

@@ -1,5 +1,5 @@
 module BioSyntax
   # Ruby gem version.
   # @return [String]
-  VERSION = '0.1.0'
+  VERSION = '0.1.1'
 end

data/lib/biosyntax.rb

@@ -436,8 +436,8 @@ module BioSyntax
               else
                 name = value.to_s.downcase
                 FORMATS[name.to_sym] ||
-                  FORMATS[name.tr('_', '-').to_sym] ||
-                  FORMATS_BY_ID[Native.format_id_from_name(name)]
+                FORMATS[name.tr('_', '-').to_sym] ||
+                FORMATS_BY_ID[Native.format_id_from_name(name)]
               end
 
       return found if found

metadata

@@ -1,22 +1,24 @@
 --- !ruby/object:Gem::Specification
 name: biosyntax
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.1
 platform: ruby
 authors:
 - kojix2
-bindir: bin
+bindir: exe
 cert_chain: []
 date: 1980-01-02 00:00:00.000000000 Z
 dependencies: []
 email: 2xijok@gmail.com
-executables: []
+executables:
+- biocat
 extensions:
 - ext/biosyntax/extconf.rb
 extra_rdoc_files: []
 files:
 - LICENSE.md
 - README.md
+- exe/biocat
 - ext/biosyntax/biosyntax.c
 - ext/biosyntax/biosyntax.h
 - ext/biosyntax/biosyntax_ext.c