lib/lederhosen/tasks/rep_reads.rb in lederhosen-0.2.9 vs lib/lederhosen/tasks/rep_reads.rb in lederhosen-0.2.10
- old
+ new
@@ -15,9 +15,10 @@
def rep_reads
input = options[:clusters]
output = options[:output]
joined_reads = options[:joined]
+ ohai "getting represntative reads for #{input} w/ reads #{joined_reads} and saving to #{output}"
# Load cluster table!
clstr_info = Helpers.load_uc_file input
clstr_counts = clstr_info[:clstr_counts] # clstr_counts[:clstr][sample.to_i] = reads
seed_to_clstrnr = clstr_info[:seed_to_clstrnr]