Skip to content

Instantly share code, notes, and snippets.

@gisleDK
Last active January 9, 2018 09:49
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save gisleDK/bb03a2efb20536eb506f779b0e302fdf to your computer and use it in GitHub Desktop.
Save gisleDK/bb03a2efb20536eb506f779b0e302fdf to your computer and use it in GitHub Desktop.
Merging several DADA2 denoised run feature tables and representative sequences for Qiime version 2017.12
#!/usr/bin/env ruby
# Qiime2: Merging several DADA2 denoised run feature tables and representative sequences for Qiime version 2017.12
# Feature tables must be name *_table.qza and representative sequences *_rep-seqs.qza
# V1.00 Written by Gisle Vestergaard
require 'pp'
require 'optparse'
tables = []
seqs = []
table = `ls | grep _table.qza`
table.each_line do |line|
line.chomp!
tables << line
end
seq = `ls | grep _rep-seqs.qza`
seq.each_line do |line|
line.chomp!
seqs << line
end
system 'qiime feature-table merge --i-tables ' + tables.each {|x| x}.join(" --i-tables ") + ' --o-merged-table merged_table.qza'
system 'qiime feature-table merge-seqs --i-data ' + seqs.each {|x| x}.join(" --i-data ") + ' --o-merged-data merged_rep-seqs.qza'
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment