miga-base 0.3.0.0 → 0.3.0.1

Sign up to get free protection for your applications and to get access to all the features.
Files changed (260) hide show
  1. checksums.yaml +4 -4
  2. data/README.md +21 -4
  3. data/actions/init.rb +258 -0
  4. data/actions/run_local.rb +1 -2
  5. data/actions/test_taxonomy.rb +4 -1
  6. data/bin/miga +8 -1
  7. data/lib/miga/dataset.rb +4 -4
  8. data/lib/miga/dataset_result.rb +7 -4
  9. data/lib/miga/version.rb +2 -2
  10. data/scripts/_distances_noref_nomulti.bash +3 -1
  11. data/scripts/clade_finding.bash +1 -1
  12. data/scripts/init.bash +1 -1
  13. data/scripts/miga.bash +1 -1
  14. data/scripts/mytaxa.bash +78 -72
  15. data/scripts/mytaxa_scan.bash +67 -62
  16. data/scripts/ogs.bash +1 -1
  17. data/scripts/trimmed_fasta.bash +4 -3
  18. data/utils/enveomics/Examples/aai-matrix.bash +66 -0
  19. data/utils/enveomics/Examples/ani-matrix.bash +66 -0
  20. data/utils/enveomics/Examples/essential-phylogeny.bash +105 -0
  21. data/utils/enveomics/Examples/unus-genome-phylogeny.bash +100 -0
  22. data/utils/enveomics/LICENSE.txt +73 -0
  23. data/utils/enveomics/Makefile +52 -0
  24. data/utils/enveomics/Manifest/Tasks/aasubs.json +103 -0
  25. data/utils/enveomics/Manifest/Tasks/blasttab.json +703 -0
  26. data/utils/enveomics/Manifest/Tasks/distances.json +161 -0
  27. data/utils/enveomics/Manifest/Tasks/fasta.json +571 -0
  28. data/utils/enveomics/Manifest/Tasks/fastq.json +208 -0
  29. data/utils/enveomics/Manifest/Tasks/graphics.json +126 -0
  30. data/utils/enveomics/Manifest/Tasks/ogs.json +339 -0
  31. data/utils/enveomics/Manifest/Tasks/other.json +746 -0
  32. data/utils/enveomics/Manifest/Tasks/remote.json +355 -0
  33. data/utils/enveomics/Manifest/Tasks/sequence-identity.json +454 -0
  34. data/utils/enveomics/Manifest/Tasks/tables.json +308 -0
  35. data/utils/enveomics/Manifest/Tasks/trees.json +68 -0
  36. data/utils/enveomics/Manifest/Tasks/variants.json +111 -0
  37. data/utils/enveomics/Manifest/categories.json +132 -0
  38. data/utils/enveomics/Manifest/examples.json +154 -0
  39. data/utils/enveomics/Manifest/tasks.json +4 -0
  40. data/utils/enveomics/Pipelines/assembly.pbs/CONFIG.mock.bash +69 -0
  41. data/utils/enveomics/Pipelines/assembly.pbs/FastA.N50.pl +56 -0
  42. data/utils/enveomics/Pipelines/assembly.pbs/FastA.filterN.pl +60 -0
  43. data/utils/enveomics/Pipelines/assembly.pbs/FastA.length.pl +38 -0
  44. data/utils/enveomics/Pipelines/assembly.pbs/README.md +189 -0
  45. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-2.bash +112 -0
  46. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-3.bash +23 -0
  47. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-4.bash +44 -0
  48. data/utils/enveomics/Pipelines/assembly.pbs/RUNME.bash +50 -0
  49. data/utils/enveomics/Pipelines/assembly.pbs/kSelector.R +37 -0
  50. data/utils/enveomics/Pipelines/assembly.pbs/newbler.pbs +68 -0
  51. data/utils/enveomics/Pipelines/assembly.pbs/newbler_preparator.pl +49 -0
  52. data/utils/enveomics/Pipelines/assembly.pbs/soap.pbs +80 -0
  53. data/utils/enveomics/Pipelines/assembly.pbs/stats.pbs +57 -0
  54. data/utils/enveomics/Pipelines/assembly.pbs/velvet.pbs +63 -0
  55. data/utils/enveomics/Pipelines/blast.pbs/01.pbs.bash +38 -0
  56. data/utils/enveomics/Pipelines/blast.pbs/02.pbs.bash +73 -0
  57. data/utils/enveomics/Pipelines/blast.pbs/03.pbs.bash +21 -0
  58. data/utils/enveomics/Pipelines/blast.pbs/BlastTab.recover_job.pl +72 -0
  59. data/utils/enveomics/Pipelines/blast.pbs/CONFIG.mock.bash +98 -0
  60. data/utils/enveomics/Pipelines/blast.pbs/FastA.split.pl +55 -0
  61. data/utils/enveomics/Pipelines/blast.pbs/README.md +127 -0
  62. data/utils/enveomics/Pipelines/blast.pbs/RUNME.bash +109 -0
  63. data/utils/enveomics/Pipelines/blast.pbs/TASK.check.bash +128 -0
  64. data/utils/enveomics/Pipelines/blast.pbs/TASK.dry.bash +16 -0
  65. data/utils/enveomics/Pipelines/blast.pbs/TASK.eo.bash +22 -0
  66. data/utils/enveomics/Pipelines/blast.pbs/TASK.pause.bash +26 -0
  67. data/utils/enveomics/Pipelines/blast.pbs/TASK.run.bash +89 -0
  68. data/utils/enveomics/Pipelines/blast.pbs/sentinel.pbs.bash +29 -0
  69. data/utils/enveomics/Pipelines/idba.pbs/README.md +49 -0
  70. data/utils/enveomics/Pipelines/idba.pbs/RUNME.bash +95 -0
  71. data/utils/enveomics/Pipelines/idba.pbs/run.pbs +56 -0
  72. data/utils/enveomics/Pipelines/trim.pbs/README.md +54 -0
  73. data/utils/enveomics/Pipelines/trim.pbs/RUNME.bash +70 -0
  74. data/utils/enveomics/Pipelines/trim.pbs/run.pbs +130 -0
  75. data/utils/enveomics/README.md +40 -0
  76. data/utils/enveomics/Scripts/AAsubs.log2ratio.rb +171 -0
  77. data/utils/enveomics/Scripts/Aln.cat.rb +162 -0
  78. data/utils/enveomics/Scripts/Aln.convert.pl +35 -0
  79. data/utils/enveomics/Scripts/AlphaDiversity.pl +152 -0
  80. data/utils/enveomics/Scripts/BlastPairwise.AAsubs.pl +102 -0
  81. data/utils/enveomics/Scripts/BlastTab.addlen.rb +61 -0
  82. data/utils/enveomics/Scripts/BlastTab.advance.bash +48 -0
  83. data/utils/enveomics/Scripts/BlastTab.best_hit_sorted.pl +55 -0
  84. data/utils/enveomics/Scripts/BlastTab.catsbj.pl +106 -0
  85. data/utils/enveomics/Scripts/BlastTab.cogCat.rb +76 -0
  86. data/utils/enveomics/Scripts/BlastTab.filter.pl +47 -0
  87. data/utils/enveomics/Scripts/BlastTab.kegg_pep2path_rest.pl +194 -0
  88. data/utils/enveomics/Scripts/BlastTab.metaxaPrep.pl +104 -0
  89. data/utils/enveomics/Scripts/BlastTab.pairedHits.rb +157 -0
  90. data/utils/enveomics/Scripts/BlastTab.recplot2.R +40 -0
  91. data/utils/enveomics/Scripts/BlastTab.seqdepth.pl +86 -0
  92. data/utils/enveomics/Scripts/BlastTab.seqdepth_ZIP.pl +119 -0
  93. data/utils/enveomics/Scripts/BlastTab.seqdepth_nomedian.pl +86 -0
  94. data/utils/enveomics/Scripts/BlastTab.subsample.pl +47 -0
  95. data/utils/enveomics/Scripts/BlastTab.sumPerHit.pl +114 -0
  96. data/utils/enveomics/Scripts/BlastTab.taxid2taxrank.pl +90 -0
  97. data/utils/enveomics/Scripts/BlastTab.topHits_sorted.rb +101 -0
  98. data/utils/enveomics/Scripts/Chao1.pl +97 -0
  99. data/utils/enveomics/Scripts/CharTable.classify.rb +234 -0
  100. data/utils/enveomics/Scripts/EBIseq2tax.rb +83 -0
  101. data/utils/enveomics/Scripts/FastA.N50.pl +56 -0
  102. data/utils/enveomics/Scripts/FastA.filter.pl +52 -0
  103. data/utils/enveomics/Scripts/FastA.filterLen.pl +28 -0
  104. data/utils/enveomics/Scripts/FastA.filterN.pl +60 -0
  105. data/utils/enveomics/Scripts/FastA.fragment.rb +92 -0
  106. data/utils/enveomics/Scripts/FastA.gc.pl +42 -0
  107. data/utils/enveomics/Scripts/FastA.interpose.pl +87 -0
  108. data/utils/enveomics/Scripts/FastA.length.pl +38 -0
  109. data/utils/enveomics/Scripts/FastA.per_file.pl +36 -0
  110. data/utils/enveomics/Scripts/FastA.qlen.pl +57 -0
  111. data/utils/enveomics/Scripts/FastA.rename.pl +65 -0
  112. data/utils/enveomics/Scripts/FastA.revcom.pl +23 -0
  113. data/utils/enveomics/Scripts/FastA.slider.pl +85 -0
  114. data/utils/enveomics/Scripts/FastA.split.pl +55 -0
  115. data/utils/enveomics/Scripts/FastA.subsample.pl +131 -0
  116. data/utils/enveomics/Scripts/FastA.tag.rb +64 -0
  117. data/utils/enveomics/Scripts/FastA.wrap.rb +48 -0
  118. data/utils/enveomics/Scripts/FastQ.filter.pl +54 -0
  119. data/utils/enveomics/Scripts/FastQ.interpose.pl +90 -0
  120. data/utils/enveomics/Scripts/FastQ.offset.pl +90 -0
  121. data/utils/enveomics/Scripts/FastQ.split.pl +53 -0
  122. data/utils/enveomics/Scripts/FastQ.tag.rb +63 -0
  123. data/utils/enveomics/Scripts/FastQ.toFastA.awk +24 -0
  124. data/utils/enveomics/Scripts/GenBank.add_fields.rb +84 -0
  125. data/utils/enveomics/Scripts/HMM.essential.rb +254 -0
  126. data/utils/enveomics/Scripts/HMMsearch.extractIds.rb +83 -0
  127. data/utils/enveomics/Scripts/JPlace.distances.rb +88 -0
  128. data/utils/enveomics/Scripts/JPlace.to_iToL.rb +306 -0
  129. data/utils/enveomics/Scripts/M5nr.getSequences.rb +81 -0
  130. data/utils/enveomics/Scripts/MeTaxa.distribution.pl +198 -0
  131. data/utils/enveomics/Scripts/MyTaxa.fragsByTax.pl +35 -0
  132. data/utils/enveomics/Scripts/MyTaxa.seq-taxrank.rb +49 -0
  133. data/utils/enveomics/Scripts/NCBIacc2tax.rb +92 -0
  134. data/utils/enveomics/Scripts/Newick.autoprune.R +27 -0
  135. data/utils/enveomics/Scripts/RAxML-EPA.to_iToL.pl +228 -0
  136. data/utils/enveomics/Scripts/RefSeq.download.bash +48 -0
  137. data/utils/enveomics/Scripts/SRA.download.bash +50 -0
  138. data/utils/enveomics/Scripts/TRIBS.plot-test.R +36 -0
  139. data/utils/enveomics/Scripts/TRIBS.test.R +39 -0
  140. data/utils/enveomics/Scripts/Table.barplot.R +30 -0
  141. data/utils/enveomics/Scripts/Table.df2dist.R +30 -0
  142. data/utils/enveomics/Scripts/Table.filter.pl +61 -0
  143. data/utils/enveomics/Scripts/Table.merge.pl +77 -0
  144. data/utils/enveomics/Scripts/Table.replace.rb +69 -0
  145. data/utils/enveomics/Scripts/Table.round.rb +63 -0
  146. data/utils/enveomics/Scripts/Table.split.pl +57 -0
  147. data/utils/enveomics/Scripts/Taxonomy.silva2ncbi.rb +227 -0
  148. data/utils/enveomics/Scripts/VCF.KaKs.rb +147 -0
  149. data/utils/enveomics/Scripts/VCF.SNPs.rb +88 -0
  150. data/utils/enveomics/Scripts/aai.rb +373 -0
  151. data/utils/enveomics/Scripts/ani.rb +362 -0
  152. data/utils/enveomics/Scripts/gi2tax.rb +103 -0
  153. data/utils/enveomics/Scripts/in_silico_GA_GI.pl +96 -0
  154. data/utils/enveomics/Scripts/lib/data/essential.hmm.gz +0 -0
  155. data/utils/enveomics/Scripts/lib/enveomics_rb/enveomics.rb +26 -0
  156. data/utils/enveomics/Scripts/lib/enveomics_rb/jplace.rb +253 -0
  157. data/utils/enveomics/Scripts/lib/enveomics_rb/og.rb +182 -0
  158. data/utils/enveomics/Scripts/lib/enveomics_rb/remote_data.rb +74 -0
  159. data/utils/enveomics/Scripts/lib/enveomics_rb/seq_range.rb +237 -0
  160. data/utils/enveomics/Scripts/lib/enveomics_rb/stat.rb +30 -0
  161. data/utils/enveomics/Scripts/lib/enveomics_rb/vcf.rb +135 -0
  162. data/utils/enveomics/Scripts/ogs.annotate.rb +88 -0
  163. data/utils/enveomics/Scripts/ogs.core-pan.rb +160 -0
  164. data/utils/enveomics/Scripts/ogs.extract.rb +125 -0
  165. data/utils/enveomics/Scripts/ogs.mcl.rb +186 -0
  166. data/utils/enveomics/Scripts/ogs.rb +104 -0
  167. data/utils/enveomics/Scripts/ogs.stats.rb +131 -0
  168. data/utils/enveomics/Scripts/rbm.rb +137 -0
  169. data/utils/enveomics/Tests/Makefile +10 -0
  170. data/utils/enveomics/Tests/Mgen_M2288.faa +3189 -0
  171. data/utils/enveomics/Tests/Mgen_M2288.fna +8282 -0
  172. data/utils/enveomics/Tests/Mgen_M2321.fna +8288 -0
  173. data/utils/enveomics/Tests/Nequ_Kin4M.faa +2970 -0
  174. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.tribs.Rdata +0 -0
  175. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.txt +7 -0
  176. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai-mat.tsv +17 -0
  177. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai.tsv +137 -0
  178. data/utils/enveomics/Tests/a_mg.cds-go.blast.tsv +123 -0
  179. data/utils/enveomics/Tests/a_mg.reads-cds.blast.tsv +200 -0
  180. data/utils/enveomics/Tests/a_mg.reads-cds.counts.tsv +55 -0
  181. data/utils/enveomics/Tests/alkB.nwk +1 -0
  182. data/utils/enveomics/Tests/anthrax-cansnp-data.tsv +13 -0
  183. data/utils/enveomics/Tests/anthrax-cansnp-key.tsv +17 -0
  184. data/utils/enveomics/Tests/hiv1.faa +59 -0
  185. data/utils/enveomics/Tests/hiv1.fna +134 -0
  186. data/utils/enveomics/Tests/hiv2.faa +70 -0
  187. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv +233 -0
  188. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.lim +1 -0
  189. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.rec +233 -0
  190. data/utils/enveomics/Tests/phyla_counts.tsv +10 -0
  191. data/utils/enveomics/Tests/primate_lentivirus.ogs +11 -0
  192. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv1.rbm +9 -0
  193. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv2.rbm +8 -0
  194. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-siv.rbm +6 -0
  195. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-hiv2.rbm +9 -0
  196. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-siv.rbm +6 -0
  197. data/utils/enveomics/Tests/primate_lentivirus.rbm/siv-siv.rbm +6 -0
  198. data/utils/enveomics/build_enveomics_r.bash +44 -0
  199. data/utils/enveomics/enveomics.R/DESCRIPTION +31 -0
  200. data/utils/enveomics/enveomics.R/NAMESPACE +35 -0
  201. data/utils/enveomics/enveomics.R/R/autoprune.R +121 -0
  202. data/utils/enveomics/enveomics.R/R/barplot.R +165 -0
  203. data/utils/enveomics/enveomics.R/R/cliopts.R +119 -0
  204. data/utils/enveomics/enveomics.R/R/df2dist.R +117 -0
  205. data/utils/enveomics/enveomics.R/R/growthcurve.R +263 -0
  206. data/utils/enveomics/enveomics.R/R/recplot.R +320 -0
  207. data/utils/enveomics/enveomics.R/R/recplot2.R +745 -0
  208. data/utils/enveomics/enveomics.R/R/tribs.R +423 -0
  209. data/utils/enveomics/enveomics.R/R/utils.R +16 -0
  210. data/utils/enveomics/enveomics.R/README.md +52 -0
  211. data/utils/enveomics/enveomics.R/data/growth.curves.rda +0 -0
  212. data/utils/enveomics/enveomics.R/data/phyla.counts.rda +0 -0
  213. data/utils/enveomics/enveomics.R/man/enve.GrowthCurve-class.Rd +30 -0
  214. data/utils/enveomics/enveomics.R/man/enve.TRIBS-class.Rd +43 -0
  215. data/utils/enveomics/enveomics.R/man/enve.TRIBS.merge.Rd +19 -0
  216. data/utils/enveomics/enveomics.R/man/enve.TRIBStest-class.Rd +37 -0
  217. data/utils/enveomics/enveomics.R/man/enve.__prune.iter.Rd +24 -0
  218. data/utils/enveomics/enveomics.R/man/enve.__prune.reduce.Rd +24 -0
  219. data/utils/enveomics/enveomics.R/man/enve.__tribs.Rd +33 -0
  220. data/utils/enveomics/enveomics.R/man/enve.barplot.Rd +64 -0
  221. data/utils/enveomics/enveomics.R/man/enve.cliopts.Rd +37 -0
  222. data/utils/enveomics/enveomics.R/man/enve.col.alpha.Rd +19 -0
  223. data/utils/enveomics/enveomics.R/man/enve.col2alpha.Rd +18 -0
  224. data/utils/enveomics/enveomics.R/man/enve.df2dist.Rd +26 -0
  225. data/utils/enveomics/enveomics.R/man/enve.df2dist.group.Rd +25 -0
  226. data/utils/enveomics/enveomics.R/man/enve.df2dist.list.Rd +26 -0
  227. data/utils/enveomics/enveomics.R/man/enve.growthcurve.Rd +49 -0
  228. data/utils/enveomics/enveomics.R/man/enve.prune.dist.Rd +28 -0
  229. data/utils/enveomics/enveomics.R/man/enve.recplot.Rd +97 -0
  230. data/utils/enveomics/enveomics.R/man/enve.recplot2-class.Rd +40 -0
  231. data/utils/enveomics/enveomics.R/man/enve.recplot2.Rd +40 -0
  232. data/utils/enveomics/enveomics.R/man/enve.recplot2.__counts.Rd +24 -0
  233. data/utils/enveomics/enveomics.R/man/enve.recplot2.__findPeak.Rd +40 -0
  234. data/utils/enveomics/enveomics.R/man/enve.recplot2.__findPeaks.Rd +18 -0
  235. data/utils/enveomics/enveomics.R/man/enve.recplot2.__peakHist.Rd +22 -0
  236. data/utils/enveomics/enveomics.R/man/enve.recplot2.__whichClosestPeak.Rd +20 -0
  237. data/utils/enveomics/enveomics.R/man/enve.recplot2.changeCutoff.Rd +18 -0
  238. data/utils/enveomics/enveomics.R/man/enve.recplot2.corePeak.Rd +18 -0
  239. data/utils/enveomics/enveomics.R/man/enve.recplot2.extractWindows.Rd +27 -0
  240. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.Rd +53 -0
  241. data/utils/enveomics/enveomics.R/man/enve.recplot2.peak-class.Rd +44 -0
  242. data/utils/enveomics/enveomics.R/man/enve.tribs.Rd +44 -0
  243. data/utils/enveomics/enveomics.R/man/enve.tribs.test.Rd +21 -0
  244. data/utils/enveomics/enveomics.R/man/enveomics.R-package.Rd +15 -0
  245. data/utils/enveomics/enveomics.R/man/growth.curves.Rd +14 -0
  246. data/utils/enveomics/enveomics.R/man/phyla.counts.Rd +13 -0
  247. data/utils/enveomics/enveomics.R/man/plot.enve.GrowthCurve.Rd +43 -0
  248. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBS.Rd +29 -0
  249. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBStest.Rd +30 -0
  250. data/utils/enveomics/enveomics.R/man/plot.enve.recplot2.Rd +71 -0
  251. data/utils/enveomics/enveomics.R/man/summary.enve.GrowthCurve.Rd +18 -0
  252. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBS.Rd +18 -0
  253. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBStest.Rd +18 -0
  254. data/utils/enveomics/enveomics.R/man/z$-methods.Rd +27 -0
  255. data/utils/enveomics/globals.mk +8 -0
  256. data/utils/enveomics/manifest.json +9 -0
  257. data/utils/index_metadata.rb +0 -0
  258. data/utils/plot-taxdist.R +0 -0
  259. data/utils/requirements.txt +19 -19
  260. metadata +242 -2
@@ -0,0 +1,208 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "FastQ.filter.pl",
5
+ "description": "Extracts a subset of sequences from a FastQ file.",
6
+ "see_also": ["FastA.filter.pl"],
7
+ "help_arg": "-h",
8
+ "options": [
9
+ {
10
+ "name": "Reverse list",
11
+ "opt": "-r",
12
+ "description": "Extracts sequences NOT present in the list."
13
+ },
14
+ {
15
+ "name": "Quiet",
16
+ "opt": "-q",
17
+ "description": "Runs quietly."
18
+ },
19
+ {
20
+ "name": "List",
21
+ "arg": "in_file",
22
+ "mandatory": true,
23
+ "description": "List of sequences to extract."
24
+ },
25
+ {
26
+ "name": "Seqs.fq",
27
+ "arg": "in_file",
28
+ "mandatory": true,
29
+ "description": "FastQ file containing the superset of sequences."
30
+ },
31
+ ">",
32
+ {
33
+ "name": "Subset.fq",
34
+ "arg": "out_file",
35
+ "mandatory": true,
36
+ "description": "FastQ file to be created."
37
+ }
38
+ ]
39
+ },
40
+ {
41
+ "task": "FastQ.interpose.pl",
42
+ "description": ["Interposes sequences in FastQ format from two files",
43
+ "into one output file. If more than two files are provided, the script",
44
+ "will interpose all the input files."],
45
+ "warn": ["Note that this script will check for the consistency of the",
46
+ "names (assuming a pair of related reads contains the same name",
47
+ "varying only in a trailing slash (/) followed by a digit. If you want",
48
+ "to turn this feature off just set the -T option to zero. If you want",
49
+ "to decrease the sampling period (to speed the script up) or increase",
50
+ "it (to make it more sensitive to errors) just change -T option",
51
+ "accordingly."],
52
+ "see_also": ["FastQ.split.pl","FastA.interpose.pl"],
53
+ "help_arg": "",
54
+ "options": [
55
+ {
56
+ "name": "Checking period",
57
+ "opt": "-T",
58
+ "arg": "integer",
59
+ "default": 1000,
60
+ "description": "Sampling period for names evaluation."
61
+ },
62
+ {
63
+ "arg": "out_file",
64
+ "mandatory": true,
65
+ "description": "Output FastQ file."
66
+ },
67
+ {
68
+ "arg": "in_file",
69
+ "mandatory": true,
70
+ "description": "First input FastQ file."
71
+ },
72
+ {
73
+ "arg": "in_file",
74
+ "mandatory": true,
75
+ "description": "Second input FastQ file."
76
+ },
77
+ {
78
+ "arg": "in_file",
79
+ "multiple_sep": " ",
80
+ "description": "Any additional input FastQ files."
81
+ }
82
+ ]
83
+ },
84
+ {
85
+ "task": "FastQ.offset.pl",
86
+ "description": ["There are several FastQ formats. This script takes a",
87
+ "FastQ in any of them, identifies the type of FastQ (this is, the",
88
+ "offset), and generates a FastQ with the given offset."],
89
+ "warn": ["Note that Solexa+64 FastQ can cause problematic values when",
90
+ "using the offset 33, since there is no equivalent in Phred+33 for",
91
+ "negative values (the range of Solexa+64 is -5 to 40)."],
92
+ "help_arg": "",
93
+ "options": [
94
+ {
95
+ "arg": "in_file",
96
+ "mandatory": true,
97
+ "description": ["Input file in FastQ format (range is automatically",
98
+ "detected)."]
99
+ },
100
+ {
101
+ "name": "Offset",
102
+ "arg": "integer",
103
+ "default": 33,
104
+ "mandatory": true,
105
+ "description": ["Offset to use for the output. Use 0 (zero) to",
106
+ "detect the input format and exit."]
107
+ },
108
+ {
109
+ "opt": "force",
110
+ "description": ["If set, turns errors into warnings and continues.",
111
+ "Out-of-range values are set to the closest range limit."]
112
+ },
113
+ ">",
114
+ {
115
+ "arg": "out_file",
116
+ "mandatory": true,
117
+ "description": ["Output file in FastQ format with the specified",
118
+ "offset."]
119
+ }
120
+ ]
121
+ },
122
+ {
123
+ "task": "FastQ.split.pl",
124
+ "description": ["Splits a FastQ file into several FastQ files. This",
125
+ "script can be used to separate interposed sister reads using any even",
126
+ "number of output files."],
127
+ "help_arg": "",
128
+ "see_also": ["FastQ.interpose.pl","FastA.split.pl"],
129
+ "options": [
130
+ {
131
+ "name": "in_file.fq",
132
+ "arg": "in_file",
133
+ "mandatory": true,
134
+ "description": "Input file in FastQ format."
135
+ },
136
+ {
137
+ "name": "out_base",
138
+ "arg": "out_file",
139
+ "mandatory": true,
140
+ "description": ["Prefix for the name of the output files. It will be",
141
+ "appended with .<i>.fastq, where <i> is a consecutive number",
142
+ "starting in 1."]
143
+ },
144
+ {
145
+ "name": "no_files",
146
+ "arg": "integer",
147
+ "default": 2,
148
+ "description": "Number of files to generate."
149
+ }
150
+ ]
151
+ },
152
+ {
153
+ "task": "FastQ.tag.rb",
154
+ "description": "Generates easy-to-parse tagged reads from FastQ files.",
155
+ "see_also": ["FastA.tag.rb"],
156
+ "help_arg": "--help",
157
+ "options": [
158
+ {
159
+ "name": "Input file",
160
+ "opt": "--in",
161
+ "arg": "in_file",
162
+ "mandatory": true,
163
+ "description": "FastQ file containing the sequences."
164
+ },
165
+ {
166
+ "name": "Output file",
167
+ "opt": "--out",
168
+ "arg": "out_file",
169
+ "mandatory": true,
170
+ "description": "FastQ to create."
171
+ },
172
+ {
173
+ "opt": "--prefix",
174
+ "arg": "string",
175
+ "description": "Prefix to use in all IDs."
176
+ },
177
+ {
178
+ "opt": "--suffix",
179
+ "arg": "string",
180
+ "description": "Suffix to use in all IDs."
181
+ },
182
+ {
183
+ "opt": "--quiet",
184
+ "description": "Run quietly (no STDERR output)."
185
+ }
186
+ ]
187
+ },
188
+ {
189
+ "task": "FastQ.toFastA.awk",
190
+ "description": "Translates FastQ files into FastA.",
191
+ "help_arg": "'' --help",
192
+ "options": [
193
+ "<",
194
+ {
195
+ "arg": "in_file",
196
+ "mandatory": true,
197
+ "description": "Input FastQ file."
198
+ },
199
+ ">",
200
+ {
201
+ "arg": "out_file",
202
+ "mandatory": true,
203
+ "description": "Output FastA file."
204
+ }
205
+ ]
206
+ }
207
+ ]
208
+ }
@@ -0,0 +1,126 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "Table.barplot.R",
5
+ "description": "Creates nice barplots from tab-delimited tables.",
6
+ "requires": [ { "r_package": "optparse" } ],
7
+ "help_arg": "--help",
8
+ "options": [
9
+ {
10
+ "name": "Input file",
11
+ "opt": "--x",
12
+ "arg": "in_file",
13
+ "mandatory": true,
14
+ "description": ["A tab-delimited file containing header (first row)",
15
+ "and row names (first column)."]
16
+ },
17
+ {
18
+ "opt": "--sizes",
19
+ "arg": "string",
20
+ "description": ["A numeric vector containing the real size of the",
21
+ "samples (columns) in the same order of the input table. If set,",
22
+ "the values are assumed to be 100%, otherwise the sum of the",
23
+ "columns is used. Separate values by commas."]
24
+ },
25
+ {
26
+ "opt": "--top",
27
+ "arg": "integer",
28
+ "default": 25,
29
+ "description": ["Maximum number of categories to display. Any",
30
+ "additional categories will be listed as 'Others'."]
31
+ },
32
+ {
33
+ "opt": "--colors-per-group",
34
+ "arg": "integer",
35
+ "default": 9,
36
+ "description": ["Number of categories in the first two saturation",
37
+ "groups of colors. The third group contains the remaining",
38
+ "categories if needed."]
39
+ },
40
+ {
41
+ "opt": "--bars-width",
42
+ "arg": "integer",
43
+ "default": 4,
44
+ "description": "Width of the barplot with respect to the legend."
45
+ },
46
+ {
47
+ "opt": "--legend-ncol",
48
+ "arg": "integer",
49
+ "default": 1,
50
+ "description": "Number of columns in the legend."
51
+ },
52
+ {
53
+ "opt": "--other-col",
54
+ "arg": "string",
55
+ "default": "#000000",
56
+ "description": "Color of the 'Others' category."
57
+ },
58
+ {
59
+ "opt": "--add-trend",
60
+ "description": ["Controls if semi-transparent areas are to be",
61
+ "plotted between the bars to connect the regions (trend regions)."]
62
+ },
63
+ {
64
+ "opt": "--organic-trend",
65
+ "description": ["Controls if the trend regions are to be smoothed",
66
+ "(curves). By default, trend regions have straight edges. If TRUE,",
67
+ "forces add.trend=TRUE."]
68
+ },
69
+ {
70
+ "opt": "--sort-by",
71
+ "arg": "string",
72
+ "default": "median",
73
+ "description": ["Any function that takes a numeric vector and",
74
+ "returns a numeric scalar. This function is applied to each row,",
75
+ "and the resulting values are used to sort the rows",
76
+ "(decreasingly). Good options include: sd, min, max, mean, median."]
77
+ },
78
+ {
79
+ "opt": "--min-report",
80
+ "arg": "integer",
81
+ "default": 101,
82
+ "description": ["Minimum percentage to report the value in the plot.",
83
+ "Any value above 100 indicates that no values are to be reported."]
84
+ },
85
+ {
86
+ "opt": "--order",
87
+ "arg": "string",
88
+ "description": ["Controls how the rows should be ordered. If empty",
89
+ "(default), sort.by is applied per row and the results are sorted",
90
+ "decreasingly. If NA, no sorting is performed, i.e., the original",
91
+ "order is respected. If a vector is provided, it is assumed to be",
92
+ "the custom order to be used (either by numeric index or by row",
93
+ "names). Separate values by commas."]
94
+ },
95
+ {
96
+ "opt": "--col",
97
+ "arg": "string",
98
+ "description": ["Colors to use. If provided, overrides the variables",
99
+ "`top` and `colors.per.group`, but `other.col` is still used if",
100
+ "the vector is insufficient for all the rows. Separate values by",
101
+ "commas. An additional palette is available when using value",
102
+ "'coto' (palette contributed by Luis (Coto) Orellana)."]
103
+ },
104
+ {
105
+ "arg": "out_file",
106
+ "mandatory": true,
107
+ "description": "Output file in PDF format."
108
+ },
109
+ {
110
+ "name": "width",
111
+ "arg": "float",
112
+ "mandatory": true,
113
+ "default": 7,
114
+ "description": "Width of the plot (in inches)."
115
+ },
116
+ {
117
+ "name": "height",
118
+ "arg": "float",
119
+ "mandatory": true,
120
+ "default": 7,
121
+ "description": "Height of the plot (in inches)."
122
+ }
123
+ ]
124
+ }
125
+ ]
126
+ }
@@ -0,0 +1,339 @@
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
+ }