Contents

#!/usr/bin/env ruby
require 'bio'
require 'optparse'
$: << File.expand_path(File.dirname(__FILE__) + '/../lib')
$: << File.expand_path('.')
path= File.expand_path(File.dirname(__FILE__) + '/../lib/bioruby-polyploid-tools.rb')
require path

options = {}
options[:identity] = 95
options[:covered] = 90
OptionParser.new do |opts|
  
  opts.banner = "Usage: filter_exonerate_by_identity.rb [options]"

  opts.on("-e", "--exo FILE", "Exonerate alignment produced by polymarker or with the following ryo: 'RESULT:\\t%S\\t%pi\\t%ql\\t%tl\\t%g\\t%V\\n'") do |o|
    options[:exo_file] = o.upcase
  end
  opts.on("-i", "--identity FLOAT", "Minimum percentage identity") do |o|
    options[:identity] = o.to_f
  end
  opts.on("-c", "--covered FLOAT", "Minimum percentage coverage") do |o|
    options[:covered] = o.to_f
  end
  
end.parse!


exo_file = options[:exo_file]
min_identity = options[:identity];
min_coverage = options[:covered]
File.foreach(exo_file) do |line|  
  aln = Bio::DB::Exonerate::Alignment.parse_custom(line)
  if aln.identity > min_identity and aln.query_coverage > min_coverage
	puts aln.line
  end
end

Version data entries

33 entries across 33 versions & 1 rubygems

Version	Path
bio-polyploid-tools-1.2.1	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-1.2.0	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-1.1.2	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-1.1.1	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-1.1.0	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-1.0.0	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.10.1	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.10.0	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.10	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.9	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.8	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.7	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.6	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.5	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.4	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.3	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.2	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.1	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.9.0	bin/filter_exonerate_by_identity.rb
bio-polyploid-tools-0.8.9	bin/filter_exonerate_by_identity.rb