miga-base 1.2.17.0 → 1.2.17.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (265) hide show
  1. checksums.yaml +4 -4
  2. data/lib/miga/version.rb +2 -2
  3. data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Archaea_SCG.hmm +41964 -0
  4. data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Bacteria_SCG.hmm +32439 -0
  5. data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Complete_SCG_DB.hmm +62056 -0
  6. data/utils/FastAAI/FastAAI +3659 -0
  7. data/utils/FastAAI/FastAAI-legacy/FastAAI +1336 -0
  8. data/utils/FastAAI/FastAAI-legacy/kAAI_v1.0_virus.py +1296 -0
  9. data/utils/FastAAI/README.md +84 -0
  10. data/utils/enveomics/Docs/recplot2.md +244 -0
  11. data/utils/enveomics/Examples/aai-matrix.bash +66 -0
  12. data/utils/enveomics/Examples/ani-matrix.bash +66 -0
  13. data/utils/enveomics/Examples/essential-phylogeny.bash +105 -0
  14. data/utils/enveomics/Examples/unus-genome-phylogeny.bash +100 -0
  15. data/utils/enveomics/LICENSE.txt +73 -0
  16. data/utils/enveomics/Makefile +52 -0
  17. data/utils/enveomics/Manifest/Tasks/aasubs.json +103 -0
  18. data/utils/enveomics/Manifest/Tasks/blasttab.json +790 -0
  19. data/utils/enveomics/Manifest/Tasks/distances.json +161 -0
  20. data/utils/enveomics/Manifest/Tasks/fasta.json +802 -0
  21. data/utils/enveomics/Manifest/Tasks/fastq.json +291 -0
  22. data/utils/enveomics/Manifest/Tasks/graphics.json +126 -0
  23. data/utils/enveomics/Manifest/Tasks/mapping.json +165 -0
  24. data/utils/enveomics/Manifest/Tasks/ogs.json +382 -0
  25. data/utils/enveomics/Manifest/Tasks/other.json +906 -0
  26. data/utils/enveomics/Manifest/Tasks/remote.json +356 -0
  27. data/utils/enveomics/Manifest/Tasks/sequence-identity.json +650 -0
  28. data/utils/enveomics/Manifest/Tasks/tables.json +308 -0
  29. data/utils/enveomics/Manifest/Tasks/trees.json +68 -0
  30. data/utils/enveomics/Manifest/Tasks/variants.json +111 -0
  31. data/utils/enveomics/Manifest/categories.json +165 -0
  32. data/utils/enveomics/Manifest/examples.json +162 -0
  33. data/utils/enveomics/Manifest/tasks.json +4 -0
  34. data/utils/enveomics/README.md +42 -0
  35. data/utils/enveomics/Scripts/AAsubs.log2ratio.rb +171 -0
  36. data/utils/enveomics/Scripts/Aln.cat.rb +221 -0
  37. data/utils/enveomics/Scripts/Aln.convert.pl +35 -0
  38. data/utils/enveomics/Scripts/AlphaDiversity.pl +152 -0
  39. data/utils/enveomics/Scripts/BedGraph.tad.rb +138 -0
  40. data/utils/enveomics/Scripts/BedGraph.window.rb +71 -0
  41. data/utils/enveomics/Scripts/BlastPairwise.AAsubs.pl +102 -0
  42. data/utils/enveomics/Scripts/BlastTab.addlen.rb +63 -0
  43. data/utils/enveomics/Scripts/BlastTab.advance.bash +48 -0
  44. data/utils/enveomics/Scripts/BlastTab.best_hit_sorted.pl +55 -0
  45. data/utils/enveomics/Scripts/BlastTab.catsbj.pl +104 -0
  46. data/utils/enveomics/Scripts/BlastTab.cogCat.rb +76 -0
  47. data/utils/enveomics/Scripts/BlastTab.filter.pl +47 -0
  48. data/utils/enveomics/Scripts/BlastTab.kegg_pep2path_rest.pl +194 -0
  49. data/utils/enveomics/Scripts/BlastTab.metaxaPrep.pl +104 -0
  50. data/utils/enveomics/Scripts/BlastTab.pairedHits.rb +157 -0
  51. data/utils/enveomics/Scripts/BlastTab.recplot2.R +48 -0
  52. data/utils/enveomics/Scripts/BlastTab.seqdepth.pl +86 -0
  53. data/utils/enveomics/Scripts/BlastTab.seqdepth_ZIP.pl +119 -0
  54. data/utils/enveomics/Scripts/BlastTab.seqdepth_nomedian.pl +86 -0
  55. data/utils/enveomics/Scripts/BlastTab.subsample.pl +47 -0
  56. data/utils/enveomics/Scripts/BlastTab.sumPerHit.pl +114 -0
  57. data/utils/enveomics/Scripts/BlastTab.taxid2taxrank.pl +90 -0
  58. data/utils/enveomics/Scripts/BlastTab.topHits_sorted.rb +123 -0
  59. data/utils/enveomics/Scripts/Chao1.pl +97 -0
  60. data/utils/enveomics/Scripts/CharTable.classify.rb +234 -0
  61. data/utils/enveomics/Scripts/EBIseq2tax.rb +83 -0
  62. data/utils/enveomics/Scripts/FastA.N50.pl +60 -0
  63. data/utils/enveomics/Scripts/FastA.extract.rb +152 -0
  64. data/utils/enveomics/Scripts/FastA.filter.pl +52 -0
  65. data/utils/enveomics/Scripts/FastA.filterLen.pl +28 -0
  66. data/utils/enveomics/Scripts/FastA.filterN.pl +60 -0
  67. data/utils/enveomics/Scripts/FastA.fragment.rb +100 -0
  68. data/utils/enveomics/Scripts/FastA.gc.pl +42 -0
  69. data/utils/enveomics/Scripts/FastA.interpose.pl +93 -0
  70. data/utils/enveomics/Scripts/FastA.length.pl +38 -0
  71. data/utils/enveomics/Scripts/FastA.mask.rb +89 -0
  72. data/utils/enveomics/Scripts/FastA.per_file.pl +36 -0
  73. data/utils/enveomics/Scripts/FastA.qlen.pl +57 -0
  74. data/utils/enveomics/Scripts/FastA.rename.pl +65 -0
  75. data/utils/enveomics/Scripts/FastA.revcom.pl +23 -0
  76. data/utils/enveomics/Scripts/FastA.sample.rb +98 -0
  77. data/utils/enveomics/Scripts/FastA.slider.pl +85 -0
  78. data/utils/enveomics/Scripts/FastA.split.pl +55 -0
  79. data/utils/enveomics/Scripts/FastA.split.rb +79 -0
  80. data/utils/enveomics/Scripts/FastA.subsample.pl +131 -0
  81. data/utils/enveomics/Scripts/FastA.tag.rb +65 -0
  82. data/utils/enveomics/Scripts/FastA.toFastQ.rb +69 -0
  83. data/utils/enveomics/Scripts/FastA.wrap.rb +48 -0
  84. data/utils/enveomics/Scripts/FastQ.filter.pl +54 -0
  85. data/utils/enveomics/Scripts/FastQ.interpose.pl +90 -0
  86. data/utils/enveomics/Scripts/FastQ.maskQual.rb +89 -0
  87. data/utils/enveomics/Scripts/FastQ.offset.pl +90 -0
  88. data/utils/enveomics/Scripts/FastQ.split.pl +53 -0
  89. data/utils/enveomics/Scripts/FastQ.tag.rb +70 -0
  90. data/utils/enveomics/Scripts/FastQ.test-error.rb +81 -0
  91. data/utils/enveomics/Scripts/FastQ.toFastA.awk +24 -0
  92. data/utils/enveomics/Scripts/GFF.catsbj.pl +127 -0
  93. data/utils/enveomics/Scripts/GenBank.add_fields.rb +84 -0
  94. data/utils/enveomics/Scripts/HMM.essential.rb +351 -0
  95. data/utils/enveomics/Scripts/HMM.haai.rb +168 -0
  96. data/utils/enveomics/Scripts/HMMsearch.extractIds.rb +83 -0
  97. data/utils/enveomics/Scripts/JPlace.distances.rb +88 -0
  98. data/utils/enveomics/Scripts/JPlace.to_iToL.rb +320 -0
  99. data/utils/enveomics/Scripts/M5nr.getSequences.rb +81 -0
  100. data/utils/enveomics/Scripts/MeTaxa.distribution.pl +198 -0
  101. data/utils/enveomics/Scripts/MyTaxa.fragsByTax.pl +35 -0
  102. data/utils/enveomics/Scripts/MyTaxa.seq-taxrank.rb +49 -0
  103. data/utils/enveomics/Scripts/NCBIacc2tax.rb +92 -0
  104. data/utils/enveomics/Scripts/Newick.autoprune.R +27 -0
  105. data/utils/enveomics/Scripts/RAxML-EPA.to_iToL.pl +228 -0
  106. data/utils/enveomics/Scripts/RecPlot2.compareIdentities.R +32 -0
  107. data/utils/enveomics/Scripts/RefSeq.download.bash +48 -0
  108. data/utils/enveomics/Scripts/SRA.download.bash +67 -0
  109. data/utils/enveomics/Scripts/TRIBS.plot-test.R +36 -0
  110. data/utils/enveomics/Scripts/TRIBS.test.R +39 -0
  111. data/utils/enveomics/Scripts/Table.barplot.R +31 -0
  112. data/utils/enveomics/Scripts/Table.df2dist.R +30 -0
  113. data/utils/enveomics/Scripts/Table.filter.pl +61 -0
  114. data/utils/enveomics/Scripts/Table.merge.pl +77 -0
  115. data/utils/enveomics/Scripts/Table.prefScore.R +60 -0
  116. data/utils/enveomics/Scripts/Table.replace.rb +69 -0
  117. data/utils/enveomics/Scripts/Table.round.rb +63 -0
  118. data/utils/enveomics/Scripts/Table.split.pl +57 -0
  119. data/utils/enveomics/Scripts/Taxonomy.silva2ncbi.rb +227 -0
  120. data/utils/enveomics/Scripts/VCF.KaKs.rb +147 -0
  121. data/utils/enveomics/Scripts/VCF.SNPs.rb +88 -0
  122. data/utils/enveomics/Scripts/aai.rb +421 -0
  123. data/utils/enveomics/Scripts/ani.rb +362 -0
  124. data/utils/enveomics/Scripts/anir.rb +137 -0
  125. data/utils/enveomics/Scripts/clust.rand.rb +102 -0
  126. data/utils/enveomics/Scripts/gi2tax.rb +103 -0
  127. data/utils/enveomics/Scripts/in_silico_GA_GI.pl +96 -0
  128. data/utils/enveomics/Scripts/lib/data/dupont_2012_essential.hmm.gz +0 -0
  129. data/utils/enveomics/Scripts/lib/data/lee_2019_essential.hmm.gz +0 -0
  130. data/utils/enveomics/Scripts/lib/enveomics.R +1 -0
  131. data/utils/enveomics/Scripts/lib/enveomics_rb/anir.rb +293 -0
  132. data/utils/enveomics/Scripts/lib/enveomics_rb/bm_set.rb +175 -0
  133. data/utils/enveomics/Scripts/lib/enveomics_rb/enveomics.rb +24 -0
  134. data/utils/enveomics/Scripts/lib/enveomics_rb/errors.rb +17 -0
  135. data/utils/enveomics/Scripts/lib/enveomics_rb/gmm_em.rb +30 -0
  136. data/utils/enveomics/Scripts/lib/enveomics_rb/jplace.rb +253 -0
  137. data/utils/enveomics/Scripts/lib/enveomics_rb/match.rb +88 -0
  138. data/utils/enveomics/Scripts/lib/enveomics_rb/og.rb +182 -0
  139. data/utils/enveomics/Scripts/lib/enveomics_rb/rbm.rb +49 -0
  140. data/utils/enveomics/Scripts/lib/enveomics_rb/remote_data.rb +74 -0
  141. data/utils/enveomics/Scripts/lib/enveomics_rb/seq_range.rb +237 -0
  142. data/utils/enveomics/Scripts/lib/enveomics_rb/stats/rand.rb +31 -0
  143. data/utils/enveomics/Scripts/lib/enveomics_rb/stats/sample.rb +152 -0
  144. data/utils/enveomics/Scripts/lib/enveomics_rb/stats.rb +3 -0
  145. data/utils/enveomics/Scripts/lib/enveomics_rb/utils.rb +74 -0
  146. data/utils/enveomics/Scripts/lib/enveomics_rb/vcf.rb +135 -0
  147. data/utils/enveomics/Scripts/ogs.annotate.rb +88 -0
  148. data/utils/enveomics/Scripts/ogs.core-pan.rb +160 -0
  149. data/utils/enveomics/Scripts/ogs.extract.rb +125 -0
  150. data/utils/enveomics/Scripts/ogs.mcl.rb +186 -0
  151. data/utils/enveomics/Scripts/ogs.rb +104 -0
  152. data/utils/enveomics/Scripts/ogs.stats.rb +131 -0
  153. data/utils/enveomics/Scripts/rbm-legacy.rb +172 -0
  154. data/utils/enveomics/Scripts/rbm.rb +108 -0
  155. data/utils/enveomics/Scripts/sam.filter.rb +148 -0
  156. data/utils/enveomics/Tests/Makefile +10 -0
  157. data/utils/enveomics/Tests/Mgen_M2288.faa +3189 -0
  158. data/utils/enveomics/Tests/Mgen_M2288.fna +8282 -0
  159. data/utils/enveomics/Tests/Mgen_M2321.fna +8288 -0
  160. data/utils/enveomics/Tests/Nequ_Kin4M.faa +2970 -0
  161. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.tribs.Rdata +0 -0
  162. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.txt +7 -0
  163. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai-mat.tsv +17 -0
  164. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai.tsv +137 -0
  165. data/utils/enveomics/Tests/a_mg.cds-go.blast.tsv +123 -0
  166. data/utils/enveomics/Tests/a_mg.reads-cds.blast.tsv +200 -0
  167. data/utils/enveomics/Tests/a_mg.reads-cds.counts.tsv +55 -0
  168. data/utils/enveomics/Tests/alkB.nwk +1 -0
  169. data/utils/enveomics/Tests/anthrax-cansnp-data.tsv +13 -0
  170. data/utils/enveomics/Tests/anthrax-cansnp-key.tsv +17 -0
  171. data/utils/enveomics/Tests/hiv1.faa +59 -0
  172. data/utils/enveomics/Tests/hiv1.fna +134 -0
  173. data/utils/enveomics/Tests/hiv2.faa +70 -0
  174. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv +233 -0
  175. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.lim +1 -0
  176. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.rec +233 -0
  177. data/utils/enveomics/Tests/low-cov.bg.gz +0 -0
  178. data/utils/enveomics/Tests/phyla_counts.tsv +10 -0
  179. data/utils/enveomics/Tests/primate_lentivirus.ogs +11 -0
  180. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv1.rbm +9 -0
  181. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv2.rbm +8 -0
  182. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-siv.rbm +6 -0
  183. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-hiv2.rbm +9 -0
  184. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-siv.rbm +6 -0
  185. data/utils/enveomics/Tests/primate_lentivirus.rbm/siv-siv.rbm +6 -0
  186. data/utils/enveomics/build_enveomics_r.bash +45 -0
  187. data/utils/enveomics/enveomics.R/DESCRIPTION +31 -0
  188. data/utils/enveomics/enveomics.R/NAMESPACE +39 -0
  189. data/utils/enveomics/enveomics.R/R/autoprune.R +167 -0
  190. data/utils/enveomics/enveomics.R/R/barplot.R +203 -0
  191. data/utils/enveomics/enveomics.R/R/cliopts.R +141 -0
  192. data/utils/enveomics/enveomics.R/R/df2dist.R +192 -0
  193. data/utils/enveomics/enveomics.R/R/growthcurve.R +349 -0
  194. data/utils/enveomics/enveomics.R/R/prefscore.R +79 -0
  195. data/utils/enveomics/enveomics.R/R/recplot.R +419 -0
  196. data/utils/enveomics/enveomics.R/R/recplot2.R +1698 -0
  197. data/utils/enveomics/enveomics.R/R/tribs.R +638 -0
  198. data/utils/enveomics/enveomics.R/R/utils.R +90 -0
  199. data/utils/enveomics/enveomics.R/README.md +81 -0
  200. data/utils/enveomics/enveomics.R/data/growth.curves.rda +0 -0
  201. data/utils/enveomics/enveomics.R/data/phyla.counts.rda +0 -0
  202. data/utils/enveomics/enveomics.R/man/cash-enve.GrowthCurve-method.Rd +16 -0
  203. data/utils/enveomics/enveomics.R/man/cash-enve.RecPlot2-method.Rd +16 -0
  204. data/utils/enveomics/enveomics.R/man/cash-enve.RecPlot2.Peak-method.Rd +16 -0
  205. data/utils/enveomics/enveomics.R/man/enve.GrowthCurve-class.Rd +25 -0
  206. data/utils/enveomics/enveomics.R/man/enve.TRIBS-class.Rd +47 -0
  207. data/utils/enveomics/enveomics.R/man/enve.TRIBS.merge.Rd +23 -0
  208. data/utils/enveomics/enveomics.R/man/enve.TRIBStest-class.Rd +47 -0
  209. data/utils/enveomics/enveomics.R/man/enve.__prune.iter.Rd +26 -0
  210. data/utils/enveomics/enveomics.R/man/enve.__prune.reduce.Rd +26 -0
  211. data/utils/enveomics/enveomics.R/man/enve.__tribs.Rd +44 -0
  212. data/utils/enveomics/enveomics.R/man/enve.barplot.Rd +111 -0
  213. data/utils/enveomics/enveomics.R/man/enve.cliopts.Rd +67 -0
  214. data/utils/enveomics/enveomics.R/man/enve.col.alpha.Rd +34 -0
  215. data/utils/enveomics/enveomics.R/man/enve.col2alpha.Rd +25 -0
  216. data/utils/enveomics/enveomics.R/man/enve.df2dist.Rd +59 -0
  217. data/utils/enveomics/enveomics.R/man/enve.df2dist.group.Rd +63 -0
  218. data/utils/enveomics/enveomics.R/man/enve.df2dist.list.Rd +46 -0
  219. data/utils/enveomics/enveomics.R/man/enve.growthcurve.Rd +78 -0
  220. data/utils/enveomics/enveomics.R/man/enve.prefscore.Rd +50 -0
  221. data/utils/enveomics/enveomics.R/man/enve.prune.dist.Rd +44 -0
  222. data/utils/enveomics/enveomics.R/man/enve.recplot.Rd +147 -0
  223. data/utils/enveomics/enveomics.R/man/enve.recplot2-class.Rd +45 -0
  224. data/utils/enveomics/enveomics.R/man/enve.recplot2.ANIr.Rd +27 -0
  225. data/utils/enveomics/enveomics.R/man/enve.recplot2.Rd +77 -0
  226. data/utils/enveomics/enveomics.R/man/enve.recplot2.__counts.Rd +28 -0
  227. data/utils/enveomics/enveomics.R/man/enve.recplot2.__peakHist.Rd +24 -0
  228. data/utils/enveomics/enveomics.R/man/enve.recplot2.__whichClosestPeak.Rd +22 -0
  229. data/utils/enveomics/enveomics.R/man/enve.recplot2.changeCutoff.Rd +22 -0
  230. data/utils/enveomics/enveomics.R/man/enve.recplot2.compareIdentities.Rd +52 -0
  231. data/utils/enveomics/enveomics.R/man/enve.recplot2.coordinates.Rd +29 -0
  232. data/utils/enveomics/enveomics.R/man/enve.recplot2.corePeak.Rd +21 -0
  233. data/utils/enveomics/enveomics.R/man/enve.recplot2.extractWindows.Rd +45 -0
  234. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.Rd +34 -0
  235. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_e.Rd +23 -0
  236. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_m.Rd +24 -0
  237. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__emauto_one.Rd +31 -0
  238. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mow_one.Rd +56 -0
  239. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mower.Rd +20 -0
  240. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.em.Rd +51 -0
  241. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.emauto.Rd +43 -0
  242. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.mower.Rd +82 -0
  243. data/utils/enveomics/enveomics.R/man/enve.recplot2.peak-class.Rd +59 -0
  244. data/utils/enveomics/enveomics.R/man/enve.recplot2.seqdepth.Rd +27 -0
  245. data/utils/enveomics/enveomics.R/man/enve.recplot2.windowDepthThreshold.Rd +36 -0
  246. data/utils/enveomics/enveomics.R/man/enve.selvector.Rd +23 -0
  247. data/utils/enveomics/enveomics.R/man/enve.tribs.Rd +68 -0
  248. data/utils/enveomics/enveomics.R/man/enve.tribs.test.Rd +28 -0
  249. data/utils/enveomics/enveomics.R/man/enve.truncate.Rd +27 -0
  250. data/utils/enveomics/enveomics.R/man/growth.curves.Rd +14 -0
  251. data/utils/enveomics/enveomics.R/man/phyla.counts.Rd +13 -0
  252. data/utils/enveomics/enveomics.R/man/plot.enve.GrowthCurve.Rd +81 -0
  253. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBS.Rd +49 -0
  254. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBStest.Rd +48 -0
  255. data/utils/enveomics/enveomics.R/man/plot.enve.recplot2.Rd +125 -0
  256. data/utils/enveomics/enveomics.R/man/summary.enve.GrowthCurve.Rd +22 -0
  257. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBS.Rd +22 -0
  258. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBStest.Rd +22 -0
  259. data/utils/enveomics/globals.mk +8 -0
  260. data/utils/enveomics/manifest.json +9 -0
  261. data/utils/multitrim/Multitrim How-To.pdf +0 -0
  262. data/utils/multitrim/README.md +67 -0
  263. data/utils/multitrim/multitrim.py +1555 -0
  264. data/utils/multitrim/multitrim.yml +13 -0
  265. metadata +268 -6
@@ -0,0 +1,382 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "ogs.annotate.rb",
5
+ "description": ["Annotates Orthology Groups (OGs) using one or more",
6
+ "reference genomes."],
7
+ "see_also": ["ogs.mcl.rb"],
8
+ "help_arg": "--help",
9
+ "options": [
10
+ {
11
+ "name": "Input file",
12
+ "opt": "--in",
13
+ "arg": "in_file",
14
+ "mandatory": true,
15
+ "description": ["Input file containing the OGs (as generated by",
16
+ "ogs.mcl.rb)."]
17
+ },
18
+ {
19
+ "name": "Output file",
20
+ "opt": "--out",
21
+ "arg": "out_file",
22
+ "mandatory": true,
23
+ "description": "Output file containing the annotated OGs."
24
+ },
25
+ {
26
+ "name": "Annotations",
27
+ "opt": "-a",
28
+ "arg": "in_file",
29
+ "mandatory": true,
30
+ "multiple_sep": ",",
31
+ "description": ["Input file(s) containing the annotations. One or",
32
+ "more tab-delimited files with the gene names in the first column",
33
+ "and the annotation in the second."]
34
+ },
35
+ {
36
+ "opt": "--format",
37
+ "arg": "string",
38
+ "default": "(\\S+)\\.txt",
39
+ "description": ["Format of the filenames for the annotation files,",
40
+ "using regex syntax."]
41
+ },
42
+ {
43
+ "opt": "--quiet",
44
+ "description": "Run quietly (no STDERR output)."
45
+ }
46
+ ]
47
+ },
48
+ {
49
+ "task": "ogs.core-pan.rb",
50
+ "description": ["Subsamples the genomes in a set of Orthology Groups",
51
+ "(OGs) and estimates the trend of core genome and pangenome sizes."],
52
+ "help_arg": "--help",
53
+ "requires": [
54
+ {
55
+ "ruby_gem": "json"
56
+ }
57
+ ],
58
+ "see_also": ["ogs.mcl.rb"],
59
+ "options": [
60
+ {
61
+ "opt": "--ogs",
62
+ "arg": "in_file",
63
+ "mandatory": true,
64
+ "description": "Input file containing the precomputed OGs."
65
+ },
66
+ {
67
+ "opt": "--summary",
68
+ "arg": "out_file",
69
+ "description": ["Output file in tabular format with summary",
70
+ "statistics."]
71
+ },
72
+ {
73
+ "opt": "--tab",
74
+ "arg": "out_file",
75
+ "description": "Output file in tabular format."
76
+ },
77
+ {
78
+ "opt": "--json",
79
+ "arg": "out_file",
80
+ "description": "Output file in JSON format."
81
+ },
82
+ {
83
+ "opt": "--replicates",
84
+ "arg": "integer",
85
+ "description": "Number of replicates to estimate.",
86
+ "default": 100
87
+ },
88
+ {
89
+ "opt": "--threads",
90
+ "arg": "integer",
91
+ "description": "Children threads to spawn."
92
+ },
93
+ {
94
+ "opt": "--quiet",
95
+ "description": "Run quietly (no STDERR output)."
96
+ }
97
+ ]
98
+ },
99
+ {
100
+ "task": "ogs.extract.rb",
101
+ "description": ["Extracts sequences of Orthology Groups (OGs) from",
102
+ "genomes (proteomes)."],
103
+ "help_arg": "--help",
104
+ "see_also": ["ogs.mcl.rb"],
105
+ "options": [
106
+ {
107
+ "name": "Input file",
108
+ "opt": "--in",
109
+ "arg": "in_file",
110
+ "mandatory": true,
111
+ "description": ["Input file containing the OGs (as generated by",
112
+ "ogs.mcl.rb)."]
113
+ },
114
+ {
115
+ "name": "Output file",
116
+ "opt": "--out",
117
+ "arg": "out_file",
118
+ "mandatory": true,
119
+ "description": "Output directory where to place extracted sequences."
120
+ },
121
+ {
122
+ "name": "Sequences",
123
+ "opt": "--seqs",
124
+ "arg": "in_file",
125
+ "mandatory": true,
126
+ "description": ["Path to the proteomes in FastA format, using '%s'",
127
+ "to denote the genome. For example: /path/to/seqs/%s.faa."]
128
+ },
129
+ {
130
+ "opt": "--core",
131
+ "arg": "float",
132
+ "description": ["Use only OGs present in at least this fraction of",
133
+ "the genomes. To use only the strict core genome*, use --core 1."],
134
+ "note": ["* To use only the unus genome (OGs with exactly one gene",
135
+ "per genome), use: --core 1 --duplicates 1."]
136
+ },
137
+ {
138
+ "opt": "--duplicates",
139
+ "arg": "integer",
140
+ "description": ["Use only OGs with less than this number of",
141
+ "in-paralogs in a genome. To use only genes without in-paralogs*,",
142
+ "use --duplicates 1."],
143
+ "note": ["* To use only the unus genome (OGs with exactly one gene",
144
+ "per genome), use: --core 1 --duplicates 1."]
145
+ },
146
+ {
147
+ "opt": "--per-genome",
148
+ "description": ["If set, the output is generated per genome. By",
149
+ "default, the output is per OG."]
150
+ },
151
+ {
152
+ "opt": "--prefix",
153
+ "description": ["If set, each sequence is prefixed with the genome",
154
+ "name (or OG number, if --per-genome) and a dash."]
155
+ },
156
+ {
157
+ "opt": "--rand",
158
+ "description": ["Get only one gene per genome per OG (random)",
159
+ "regardless of in-paralogs. By default all genes are extracted."]
160
+ },
161
+ {
162
+ "opt": "--first",
163
+ "description": ["Get only one gene per genome per OG (first)",
164
+ "regardless of in-paralogs. By default all genes are extracted.",
165
+ "Takes precedence over --rand."]
166
+ },
167
+ {
168
+ "opt": "--quiet",
169
+ "description": "Run quietly (no STDERR output)."
170
+ }
171
+ ]
172
+ },
173
+ {
174
+ "task": "ogs.mcl.rb",
175
+ "description": ["Identifies Orthology Groups (OGs) in Reciprocal Best",
176
+ "Matches (RBM) between all pairs in a collection of genomes, using the",
177
+ "Markov Cluster Algorithm."],
178
+ "see_also": ["ogs.annotate.rb", "ogs.core-pan.rb", "ogs.extract.rb",
179
+ "ogs.stats.rb"],
180
+ "cite": [["Enright et al, 2002, NAR",
181
+ "http://dx.doi.org/10.1093/nar/30.7.1575"]],
182
+ "help_arg": "--help",
183
+ "options": [
184
+ {
185
+ "opt": "--out",
186
+ "arg": "out_file",
187
+ "mandatory": true,
188
+ "description": "Output file containing the detected OGs."
189
+ },
190
+ {
191
+ "opt": "--dir",
192
+ "arg": "in_dir",
193
+ "description": "Directory containing the RBM files.",
194
+ "note": "Mandatory, unless --abc is set to a non-empty file."
195
+ },
196
+ {
197
+ "opt": "--format",
198
+ "arg": "string",
199
+ "description": ["Format of the filenames for the RBM files (within",
200
+ "--dir), using regex syntax."],
201
+ "default": "(\\S+)-(\\S+)\\.rbm"
202
+ },
203
+ {
204
+ "opt": "--inflation",
205
+ "arg": "float",
206
+ "description": "Inflation parameter for MCL clustering.",
207
+ "default": 1.5
208
+ },
209
+ {
210
+ "opt": "--blind",
211
+ "description": ["If set, computes clusters without taking bitscore",
212
+ "into account."]
213
+ },
214
+ {
215
+ "opt": "--evalue",
216
+ "description": ["If set, uses the e-value to weight edges, instead",
217
+ "of the default Bit-Score."]
218
+ },
219
+ {
220
+ "opt": "--identity",
221
+ "description": ["If set, uses the identity to weight edges, instead",
222
+ "of the default Bit-Score."]
223
+ },
224
+ {
225
+ "opt": "--best-match",
226
+ "description": ["If set, it assumes best-matches instead reciprocal",
227
+ "best matches."]
228
+ },
229
+ {
230
+ "opt": "--mcl-bin",
231
+ "arg": "in_dir",
232
+ "description": ["Path to the directory containing the mcl binaries.",
233
+ "By default, assumed to be in the PATH."]
234
+ },
235
+ {
236
+ "name": "abc",
237
+ "arg": "out_file",
238
+ "opt": "--abc",
239
+ "description": "Use this abc file instead of a temporal file."
240
+ },
241
+ {
242
+ "opt": "--threads",
243
+ "arg": "integer",
244
+ "default": 2,
245
+ "description": "Number of threads to use."
246
+ },
247
+ {
248
+ "opt": "--quiet",
249
+ "description": "Run quietly (no STDERR output)."
250
+ }
251
+ ]
252
+ },
253
+ {
254
+ "task": "ogs.rb",
255
+ "description": ["Identifies Orthology Groups (OGs) in Reciprocal Best",
256
+ "Matches (RBM) between all pairs in a collection of genomes."],
257
+ "warn": ["This script suffers from chaining effect and is very",
258
+ "sensitive to spurious connections, because it applies a greedy",
259
+ "clustering algorithm. For most practical purposes, the use of this",
260
+ "script is discouraged and `ogs.mcl.rb` should be preferred."],
261
+ "help_arg": "--help",
262
+ "see_also": ["ogs.mcl.rb"],
263
+ "options": [
264
+ {
265
+ "opt": "--out",
266
+ "mandatory": true,
267
+ "arg": "out_file",
268
+ "description": "Output file containing the detected OGs."
269
+ },
270
+ {
271
+ "opt": "--dir",
272
+ "arg": "in_dir",
273
+ "description": "Directory containing the RBM files.",
274
+ "note": "Required unless --pre-ogs is passed."
275
+ },
276
+ {
277
+ "opt": "--pre-ogs",
278
+ "arg": "in_file",
279
+ "multiple_sep": ",",
280
+ "description": "Pre-computed OGs file(s), separated by commas."
281
+ },
282
+ {
283
+ "opt": "--unchecked",
284
+ "description": "Do not check internal redundancy in OGs."
285
+ },
286
+ {
287
+ "opt": "--format",
288
+ "arg": "string",
289
+ "default": "(\\S+)-(\\S+)\\.rbm",
290
+ "description": ["Format of the filenames for the RBM files (within",
291
+ "-d), using regex syntax."]
292
+ },
293
+ {
294
+ "opt": "--quiet",
295
+ "description": "Run quietly (no STDERR output)."
296
+ }
297
+ ]
298
+ },
299
+ {
300
+ "task": "ogs.stats.rb",
301
+ "description": ["Estimates some descriptive statistics on a set of",
302
+ "Orthology Groups (OGs)."],
303
+ "see_also": ["ogs.mcl.rb"],
304
+ "help_arg": "--help",
305
+ "requires": [ { "ruby_gem": "json" } ],
306
+ "options": [
307
+ {
308
+ "opt": "--ogs",
309
+ "arg": "in_file",
310
+ "mandatory": true,
311
+ "description": "Input file containing the precomputed OGs."
312
+ },
313
+ {
314
+ "opt": "--json",
315
+ "arg": "out_file",
316
+ "description": "Output file in JSON format."
317
+ },
318
+ {
319
+ "opt": "--tab",
320
+ "arg": "out_file",
321
+ "description": "Output file in tabular format."
322
+ },
323
+ {
324
+ "opt": "--transposed-tab",
325
+ "arg": "out_file",
326
+ "description": "Output file in transposed tabular format."
327
+ },
328
+ {
329
+ "opt": "--auto",
330
+ "description": "Run completely quiertly (no STDERR or STDOUT)."
331
+ },
332
+ {
333
+ "opt": "--quiet",
334
+ "description": "Run quietly (no STDERR output)."
335
+ }
336
+ ]
337
+ },
338
+ {
339
+ "task": "clust.rand.rb",
340
+ "description": ["Calculates the Rand Index and the Adjusted Rand Index",
341
+ "between two clusterings. The clustering format is a raw text file",
342
+ "with one cluster per line, each defined as comma-delimited members,",
343
+ "and a header line (ignored). Note that this is equivalent to the OGs",
344
+ "format for 1 genome."],
345
+ "see_also": ["ogs.mcl.rb"],
346
+ "help_arg": "--help",
347
+ "cite": [
348
+ ["Rand, 1971, J Am Stat Assoc",
349
+ "https://doi.org/10.2307%2F2284239"],
350
+ ["Hubert & Arabie, 1985, J Classif",
351
+ "https://doi.org/10.1007%2FBF01908075"]
352
+ ],
353
+ "options": [
354
+ {
355
+ "name": "Input file 1",
356
+ "opt": "--clust1",
357
+ "arg": "in_file",
358
+ "mandatory": true,
359
+ "description": "First input file."
360
+ },
361
+ {
362
+ "name": "Input file 2",
363
+ "opt": "--clust2",
364
+ "arg": "in_file",
365
+ "mandatory": true,
366
+ "description": "Second input file."
367
+ },
368
+ {
369
+ "name": "Precision",
370
+ "opt": "--prec",
371
+ "arg": "integer",
372
+ "description": "Precision to report.",
373
+ "default": 6
374
+ },
375
+ {
376
+ "opt": "--quiet",
377
+ "description": "Run quietly (no STDERR output)."
378
+ }
379
+ ]
380
+ }
381
+ ]
382
+ }