Contents

require 'set'

module Lederhosen

  class CLI

    desc 'separate_unclassified',
         'separate unclassified reads (with or without strict pairing)'

    method_option :uc_file, :type => :string, :required => true
    method_option :reads,   :type => :string, :required => true
    method_option :output,  :type => :string, :required => true
    method_option :strict,  :type => :string, :default => false

    def separate_unclassified
      uc_file = options[:uc_file]
      reads   = options[:reads]
      output  = options[:output]
      strict  = options[:strict]

      unclassifieds = Set.new
      handle = File.open(uc_file)
      uc = UCParser.new(handle)

      if not strict
        uc.each do |result|
          unclassifieds << result.query if result.miss?
        end

      elsif strict

        uc.each_slice(2) do |left, right|
          if left.miss? || right.miss? # at least one is a miss
            unclassifieds << left.query
            unclassifieds << right.query
          # both are hits, check taxonomies
          else
            ta = parse_taxonomy(right.target)
            tb = parse_taxonomy(left.target)
            # inconsistent assignment or at least one is a miss
            if (ta[strict] != tb[strict])
              unclassifieds << left.query
              unclassifieds << right.query
            end
          end
        end

      end

      ohai "found #{unclassifieds.size} unclassified #{'(strict pairing)' if strict} reads."

      handle.close

      # open fasta file, output unclassified reads
      out = File.open(output, 'w')
      Dna.new(File.open(reads)).each do |record|
        if unclassifieds.include? record.name
          out.puts record
        end
      end
      out.close

    end
  end
end

Version data entries

5 entries across 5 versions & 1 rubygems

Version	Path
lederhosen-2.0.4	lib/lederhosen/tasks/separate_unclassified.rb
lederhosen-2.0.3	lib/lederhosen/tasks/separate_unclassified.rb
lederhosen-2.0.2	lib/lederhosen/tasks/separate_unclassified.rb
lederhosen-2.0.1	lib/lederhosen/tasks/separate_unclassified.rb
lederhosen-2.0.0	lib/lederhosen/tasks/separate_unclassified.rb