Guest User

Untitled

a guest
Jan 16th, 2018
165
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.74 KB | None | 0 0
  1. #!/usr/bin/env ruby
  2. # Qiime2: Merging several DADA2 denoised run feature tables and representative sequences for Qiime version 2017.12
  3. # Feature tables must be name *_table.qza and representative sequences *_rep-seqs.qza
  4. # V1.00 Written by Gisle Vestergaard
  5.  
  6. require 'pp'
  7. require 'optparse'
  8.  
  9. tables = []
  10. seqs = []
  11.  
  12. table = `ls | grep _table.qza`
  13. table.each_line do |line|
  14. line.chomp!
  15. tables << line
  16. end
  17.  
  18. seq = `ls | grep _rep-seqs.qza`
  19. seq.each_line do |line|
  20. line.chomp!
  21. seqs << line
  22. end
  23.  
  24. system 'qiime feature-table merge --i-tables ' + tables.each {|x| x}.join(" --i-tables ") + ' --o-merged-table merged_table.qza'
  25. system 'qiime feature-table merge-seqs --i-data ' + seqs.each {|x| x}.join(" --i-data ") + ' --o-merged-data merged_rep-seqs.qza'
Add Comment
Please, Sign In to add comment