X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Ffindsim_seq;h=edfe980f5cb1e63c8bbef2851bdbd4a41baa302b;hb=4eed3b661f80010e9dd2936951c747ae4d7382a9;hp=fec7772269aed40b6fd0bccde1bf1b9a78faeaa3;hpb=652917ae712661f5679addc6e0f087c342a20167;p=biopieces.git diff --git a/bp_bin/findsim_seq b/bp_bin/findsim_seq index fec7772..edfe980 100755 --- a/bp_bin/findsim_seq +++ b/bp_bin/findsim_seq @@ -34,14 +34,15 @@ require 'maasha/fasta' require 'maasha/findsim' casts = [] -casts << {:long=>'database', :short=>'d', :type=>'file!', :mandatory=>true, :default=>nil, :allowed=>nil, :disallowed=>nil} -casts << {:long=>'kmer', :short=>'k', :type=>'uint', :mandatory=>false, :default=>8, :allowed=>"4,5,6,7,8,9", :disallowed=>nil} -casts << {:long=>'step', :short=>'s', :type=>'uint', :mandatory=>false, :default=>1, :allowed=>nil, :disallowed=>"0"} -casts << {:long=>'min_score', :short=>'m', :type=>'float', :mandatory=>false, :default=>0.5, :allowed=>nil, :disallowed=>nil} -casts << {:long=>'max_hits', :short=>'h', :type=>'uint', :mandatory=>false, :default=>7, :allowed=>nil, :disallowed=>"0"} -casts << {:long=>'query_ids', :short=>'Q', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} -casts << {:long=>'subject_ids', :short=>'S', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} -casts << {:long=>'realign', :short=>'r', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'database', :short=>'d', :type=>'file!', :mandatory=>true, :default=>nil, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'kmer', :short=>'k', :type=>'uint', :mandatory=>false, :default=>8, :allowed=>"4,5,6,7,8,9", :disallowed=>nil} +casts << {:long=>'step', :short=>'s', :type=>'uint', :mandatory=>false, :default=>1, :allowed=>nil, :disallowed=>"0"} +casts << {:long=>'min_score', :short=>'m', :type=>'float', :mandatory=>false, :default=>0.5, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'max_hits', :short=>'h', :type=>'uint', :mandatory=>false, :default=>20, :allowed=>nil, :disallowed=>"0"} +casts << {:long=>'max_diversity', :short=>'M', :type=>'float', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'query_ids', :short=>'Q', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'subject_ids', :short=>'S', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} +casts << {:long=>'realign', :short=>'r', :type=>'flag', :mandatory=>false, :default=>nil, :allowed=>nil, :disallowed=>nil} options = Biopieces.options_parse(ARGV, casts) tmpdir = Biopieces.mktmpdir