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,161 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "TRIBS.test.R",
5
+ "description": ["Estimates the empirical difference between all the",
6
+ "distances in a set of objects and a subset, together with its",
7
+ "statistical significance."],
8
+ "help_arg": "--help",
9
+ "requires": [
10
+ { "r_package": "optparse" },
11
+ { "r_package": "enveomics.R" }
12
+ ],
13
+ "see_also": ["TRIBS.plot-test.R"],
14
+ "options": [
15
+ {
16
+ "opt": "--dist",
17
+ "arg": "in_file",
18
+ "mandatory": true,
19
+ "description": "A tab-delimited matrix of distances."
20
+ },
21
+ {
22
+ "opt": "--selection",
23
+ "arg": "in_file",
24
+ "mandatory": true,
25
+ "description": "A list of names with the selection to evaluate."
26
+ },
27
+ {
28
+ "opt": "--replicates",
29
+ "arg": "integer",
30
+ "default": 1000,
31
+ "description": "Number of replications per point."
32
+ },
33
+ {
34
+ "opt": "--summary-fx",
35
+ "arg": "string",
36
+ "description": ["Function to summarize the distance distributions in",
37
+ "a given replicate. By default, the median distance is estimated."]
38
+ },
39
+ {
40
+ "opt": "--dist-method",
41
+ "arg": "string",
42
+ "default": "euclidean",
43
+ "description": ["Distance method between random points and samples",
44
+ "in the transformed space. See the R function 'dist'."]
45
+ },
46
+ {
47
+ "opt": "--dimensions",
48
+ "arg": "integer",
49
+ "description": ["Dimensions to use in the NMDS. By default, 5% of",
50
+ "the selection length."]
51
+ },
52
+ {
53
+ "opt": "--threads",
54
+ "arg": "integer",
55
+ "default": 2,
56
+ "description": "Number of threads to use."
57
+ },
58
+ {
59
+ "opt": "--verbosity",
60
+ "arg": "integer",
61
+ "default": 1,
62
+ "description": ["Verbosity. Use 0 to run quietly, increase for",
63
+ "additional information."]
64
+ },
65
+ {
66
+ "arg": "out_file",
67
+ "mandatory": true,
68
+ "description": "Output file with the TRIBS test R object."
69
+ },
70
+ {
71
+ "name": "Bins",
72
+ "arg": "integer",
73
+ "default": 50,
74
+ "description": "Number of bins to evaluate in the range of distances."
75
+ }
76
+ ]
77
+ },
78
+ {
79
+ "task": "TRIBS.plot-test.R",
80
+ "description": "Plots an `enve.TRIBStest` object.",
81
+ "requires": [
82
+ { "r_package": "optparse" },
83
+ { "r_package": "enveomics.R" }
84
+ ],
85
+ "see_also": ["TRIBS.test.R"],
86
+ "help_arg": "--help",
87
+ "options": [
88
+ {
89
+ "name": "Input TRIBS Test",
90
+ "opt": "--x",
91
+ "arg": "in_file",
92
+ "mandatory": true,
93
+ "description": "'enve.TRIBStest' object to plot."
94
+ },
95
+ {
96
+ "opt": "--type",
97
+ "arg": "select",
98
+ "values": ["overlap","difference"],
99
+ "default": "overlap",
100
+ "description": ["What to plot. 'overlap' generates a plot of the two",
101
+ "contrasting empirical PDFs (to compare against each other),",
102
+ "'difference' produces a plot of the differences between the",
103
+ "empirical PDFs (to compare against zero)."]
104
+ },
105
+ {
106
+ "opt": "--col",
107
+ "arg": "string",
108
+ "default": "#00000044",
109
+ "description": "Main color of the plot if type='difference'."
110
+ },
111
+ {
112
+ "opt": "--col1",
113
+ "arg": "string",
114
+ "description": ["First color of the plot if 'Type' is overlap. 'Col'",
115
+ "by default."]
116
+ },
117
+ {
118
+ "opt": "--col2",
119
+ "arg": "string",
120
+ "default": "#44001144",
121
+ "description": "Second color of the plot if 'Type' is overlap."
122
+ },
123
+ {
124
+ "opt": "--ylab",
125
+ "arg": "string",
126
+ "default": "Probability",
127
+ "description": "Y-axis label."
128
+ },
129
+ {
130
+ "opt": "--xlim",
131
+ "arg": "string",
132
+ "description": "X-axis limits, separate values by comma."
133
+ },
134
+ {
135
+ "opt": "--ylim",
136
+ "arg": "string",
137
+ "description": "Y-axis limits, separate values by comma."
138
+ },
139
+ {
140
+ "arg": "out_file",
141
+ "mandatory": true,
142
+ "description": "Output file in PDF format."
143
+ },
144
+ {
145
+ "name": "Width",
146
+ "arg": "float",
147
+ "default": 7.0,
148
+ "mandatory": true,
149
+ "description": "Width of the plot in inches."
150
+ },
151
+ {
152
+ "name": "Height",
153
+ "arg": "float",
154
+ "default": 7.0,
155
+ "mandatory": true,
156
+ "description": "Height of the plot in inches."
157
+ }
158
+ ]
159
+ }
160
+ ]
161
+ }
@@ -0,0 +1,571 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "FastA.N50.pl",
5
+ "description": ["Calculates the N50 value of a set of sequences.",
6
+ "Alternatively, it can calculate other N** values. It also calculates",
7
+ "the total number of sequences and the total added length."],
8
+ "help_arg": "",
9
+ "see_also": ["FastA.length.pl"],
10
+ "options": [
11
+ {
12
+ "name": "Sequences",
13
+ "arg": "in_file",
14
+ "mandatory": true,
15
+ "description": "A FastA file containing the sequences."
16
+ },
17
+ {
18
+ "name": "Minimum length",
19
+ "arg": "integer",
20
+ "description": "Minimum length to take into consideration.",
21
+ "default": 0
22
+ },
23
+ {
24
+ "name": "NXX",
25
+ "arg": "integer",
26
+ "description": "Value NXX to calculate.",
27
+ "default": 50
28
+ }
29
+ ]
30
+ },
31
+ {
32
+ "task": "FastA.filter.pl",
33
+ "description": "Extracts a subset of sequences from a FastA file.",
34
+ "help_arg": "-h",
35
+ "see_also": ["FastQ.filter.pl"],
36
+ "options": [
37
+ {
38
+ "opt": "-r",
39
+ "description": ["Reverse list. Extracts sequences NOT present in the",
40
+ "list."]
41
+ },
42
+ {
43
+ "opt": "-q",
44
+ "description": "Runs quietly."
45
+ },
46
+ {
47
+ "name": "List",
48
+ "arg": "in_file",
49
+ "mandatory": true,
50
+ "description": "List of sequences to extract."
51
+ },
52
+ {
53
+ "name": "Sequences",
54
+ "arg": "in_file",
55
+ "mandatory": true,
56
+ "description": "FastA file containing the superset of sequences."
57
+ },
58
+ ">",
59
+ {
60
+ "name": "Subset",
61
+ "arg": "out_file",
62
+ "mandatory": true,
63
+ "description": "FastA file to be created."
64
+ }
65
+ ]
66
+ },
67
+ {
68
+ "task": "FastA.filterLen.pl",
69
+ "description": "Filters a multi-FastA file by length.",
70
+ "requires": [ { "perl_lib": "Bio::SeqIO" } ],
71
+ "cite": [["Stajich et al, 2002, GRes",
72
+ "http://dx.doi.org/10.1101/gr.361602"]],
73
+ "help_arg": "",
74
+ "options": [
75
+ {
76
+ "arg": "in_file",
77
+ "mandatory": true,
78
+ "description": "Input FastA file."
79
+ },
80
+ {
81
+ "name": "Minimum length",
82
+ "arg": "integer",
83
+ "mandatory": true,
84
+ "description": "Minimum length to report a sequence."
85
+ },
86
+ ">",
87
+ {
88
+ "arg": "out_file",
89
+ "mandatory": true,
90
+ "description": "Filtered FastA file."
91
+ }
92
+ ]
93
+ },
94
+ {
95
+ "task": "FastA.filterN.pl",
96
+ "description": ["Filter sequences by N-content and presence of long",
97
+ "homopolymers."],
98
+ "help_arg": "",
99
+ "options": [
100
+ {
101
+ "name": "Sequences",
102
+ "arg": "in_file",
103
+ "mandatory": true,
104
+ "description": "Input file in FastA format."
105
+ },
106
+ {
107
+ "name": "Content",
108
+ "arg": "float",
109
+ "default": 0.5,
110
+ "mandatory": true,
111
+ "description": ["A number between 0 and 1 indicating the maximum",
112
+ "proportion of Ns (1 to turn off, 0.5 by default)."]
113
+ },
114
+ {
115
+ "name": "Stretch",
116
+ "arg": "integer",
117
+ "default": 100,
118
+ "description": ["A number indicating the maximum number of",
119
+ "consecutive identical nucleotides allowed (0 to turn off, 100 by",
120
+ "default)."]
121
+ },
122
+ ">",
123
+ {
124
+ "name": "Filtered",
125
+ "arg": "out_file",
126
+ "mandatory": true,
127
+ "description": "Filtered set of sequences."
128
+ }
129
+ ]
130
+ },
131
+ {
132
+ "task": "FastA.gc.pl",
133
+ "description": "Estimates the G+C content of sequences.",
134
+ "help_arg": "",
135
+ "options": [
136
+ {
137
+ "name": "seqs.fa",
138
+ "arg": "in_file",
139
+ "multiple_sep": " ",
140
+ "mandatory": true,
141
+ "description": "One or more FastA files."
142
+ },
143
+ ">",
144
+ {
145
+ "name": "gc.txt",
146
+ "arg": "out_file",
147
+ "mandatory": true,
148
+ "description": "A table with the G+C content of the sequences."
149
+ }
150
+ ]
151
+ },
152
+ {
153
+ "task": "FastA.interpose.pl",
154
+ "description": ["Interpose sequences in FastA format from two files into",
155
+ "one output file. If more than two files are provided, the script will",
156
+ "interpose all the input files."],
157
+ "see_also": ["FastA.split.pl","FastQ.interpose.pl"],
158
+ "warn": ["Please note that this script will check for the consistency of",
159
+ "the names (assuming a pair of related reads contains the same name",
160
+ "varying only in a trailing slash (/) followed by a digit. If you want",
161
+ "to turn this feature off just set the $eval_T variable to zero. If",
162
+ "you want to decrease the sampling period (to speed the script up) or",
163
+ "increase it (to make it more sensitive to errors) just change $eval_T",
164
+ "accordingly."],
165
+ "help_arg": "",
166
+ "options": [
167
+ {
168
+ "arg": "out_file",
169
+ "mandatory": true,
170
+ "description": "Output FastA file."
171
+ },
172
+ {
173
+ "name": "Input FastA 1",
174
+ "arg": "in_file",
175
+ "mandatory": true,
176
+ "description": "First FastA file."
177
+ },
178
+ {
179
+ "name": "Input FastA 2",
180
+ "arg": "in_file",
181
+ "mandatory": true,
182
+ "description": "Second FastA file."
183
+ },
184
+ {
185
+ "name": "Additional FastA",
186
+ "arg": "in_file",
187
+ "multiple_sep": " ",
188
+ "description": "Any additional FastA files (or none)."
189
+ }
190
+ ]
191
+ },
192
+ {
193
+ "task": "FastA.length.pl",
194
+ "description": "Returns the length of sequences in (multi-)FastA.",
195
+ "see_also": ["FastA.N50.pl"],
196
+ "help_arg": "",
197
+ "options": [
198
+ {
199
+ "arg": "in_file",
200
+ "multiple_sep": " ",
201
+ "mandatory": true,
202
+ "description": "One or more FastA files."
203
+ },
204
+ ">",
205
+ {
206
+ "arg": "out_file",
207
+ "mandatory": true,
208
+ "description": "A table with the lengths of the sequences."
209
+ }
210
+ ]
211
+ },
212
+ {
213
+ "task": "FastA.qlen.pl",
214
+ "description": ["Calculates the quartiles of the length in a set of",
215
+ "sequences. The Q2 is also known as the median. Q0 is the minimum",
216
+ "length, and Q4 is the maximum length. It also calculates TOTAL, the",
217
+ "added length of the sequences in the file, and AVG, the average",
218
+ "length."],
219
+ "help_arg": "",
220
+ "options": [
221
+ {
222
+ "name": "Sequences",
223
+ "arg": "in_file",
224
+ "mandatory": true,
225
+ "description": "A FastA file containing the sequences."
226
+ },
227
+ {
228
+ "name": "Minimum Length",
229
+ "arg": "integer",
230
+ "default": 0,
231
+ "description": "The minimum length to take into consideration."
232
+ }
233
+ ]
234
+ },
235
+ {
236
+ "task": "FastA.rename.pl",
237
+ "description": "Renames a set of sequences in FastA format.",
238
+ "help_arg": "-h",
239
+ "options": [
240
+ {
241
+ "name": "Filter list",
242
+ "opt": "-f",
243
+ "description": "Ignore sequences NOT present in the list."
244
+ },
245
+ {
246
+ "name": "Quiet",
247
+ "opt": "-q",
248
+ "description": "Runs quietly."
249
+ },
250
+ {
251
+ "name": "List.txt",
252
+ "arg": "in_file",
253
+ "mandatory": true,
254
+ "description": ["Tab-delimited list of sequences, with the original",
255
+ "ID in the first column and the ID to use in the second."]
256
+ },
257
+ {
258
+ "name": "Seqs.fa",
259
+ "arg": "in_file",
260
+ "mandatory": true,
261
+ "description": "FastA file containing the superset of sequences."
262
+ },
263
+ ">",
264
+ {
265
+ "name": "Renamed.fa",
266
+ "arg": "out_file",
267
+ "mandatory": true,
268
+ "description": "FastA file to be created."
269
+ }
270
+ ]
271
+ },
272
+ {
273
+ "task": "FastA.revcom.pl",
274
+ "description": "Reverse-complement sequences in FastA format.",
275
+ "requires": [ { "perl_lib": "Bio::SeqIO" } ],
276
+ "cite": [["Stajich et al, 2002, GRes",
277
+ "http://dx.doi.org/10.1101/gr.361602"]],
278
+ "help_arg": "--help",
279
+ "options": [
280
+ "<",
281
+ {
282
+ "arg": "in_file",
283
+ "mandatory": true,
284
+ "description": "Input file in FastA format."
285
+ },
286
+ ">",
287
+ {
288
+ "arg": "out_file",
289
+ "mandatory": true,
290
+ "description": "Output file in FastA format."
291
+ }
292
+ ]
293
+ },
294
+ {
295
+ "task": "FastA.slider.pl",
296
+ "description": "Slices sequences in fixed- or variable-length windows.",
297
+ "help_arg": "-help",
298
+ "options": [
299
+ {
300
+ "opt": "-seq",
301
+ "arg": "in_file",
302
+ "mandatory": true,
303
+ "description": "Input file in FastA format."
304
+ },
305
+ {
306
+ "opt": "-out",
307
+ "arg": "out_file",
308
+ "mandatory": true,
309
+ "description": "Output file in FastA format."
310
+ },
311
+ {
312
+ "opt": "-win",
313
+ "arg": "integer",
314
+ "default": 18,
315
+ "description": "Window size."
316
+ },
317
+ {
318
+ "opt": "-step",
319
+ "arg": "integer",
320
+ "default": 1,
321
+ "description": "Step size."
322
+ },
323
+ {
324
+ "name": "Length error",
325
+ "opt": "-lerr",
326
+ "arg": "integer",
327
+ "default": 2,
328
+ "description": "Expected error in chunks length."
329
+ },
330
+ {
331
+ "opt": "-comm",
332
+ "arg": "select",
333
+ "values": [0,1],
334
+ "default": 0,
335
+ "description": ["Generate FastA comments (leaded by semi-colon) to",
336
+ "separate input sequences. Set to 0 for no-comments, to 1 for",
337
+ "comments."]
338
+ },
339
+ {
340
+ "opt": "-short",
341
+ "arg": "select",
342
+ "values": [0,1],
343
+ "default": 0,
344
+ "description": ["Use chunks shorter than the window size 'as is'.",
345
+ "Set to 0 to discard those chunks, to 1 to use them."]
346
+ }
347
+ ]
348
+ },
349
+ {
350
+ "task": "FastA.split.pl",
351
+ "description": "Splits a FastA file into two or more files.",
352
+ "see_also": ["FastA.interpose.pl","FastQ.split.pl"],
353
+ "help_arg": "",
354
+ "options": [
355
+ {
356
+ "arg": "in_file",
357
+ "mandatory": true,
358
+ "description": "Input file in FastA format."
359
+ },
360
+ {
361
+ "name": "Out base",
362
+ "arg": "out_file",
363
+ "mandatory": true,
364
+ "description": ["Prefix for the name of the output files. It will be",
365
+ "appended with .<i>.fa, where <i> is a consecutive number starting",
366
+ "in 1."]
367
+ },
368
+ {
369
+ "name": "Number of files",
370
+ "arg": "integer",
371
+ "default": 12,
372
+ "description": "Number of files to generate."
373
+ }
374
+ ]
375
+ },
376
+ {
377
+ "task": "FastA.subsample.pl",
378
+ "description": "Subsamples a set of sequences.",
379
+ "help_arg": "-h",
380
+ "options": [
381
+ {
382
+ "name": "Fraction",
383
+ "opt": "-f",
384
+ "arg": "string",
385
+ "default": "10",
386
+ "description": ["Fraction of the library to be sampled (as",
387
+ "percentage). It can include several values (separated by comma),",
388
+ "as well as ranges of values in the form 'from-to/by'. For",
389
+ "example, the -f value 1-5/1,10-50/10,75,99 will produce 12",
390
+ "subsamples with expected fractions 1%, 2%, 3%, 4%, 5%, 10%, 20%,",
391
+ "30%, 40%, 50%, 75%, and 99%."]
392
+ },
393
+ {
394
+ "name": "Replicates",
395
+ "opt": "-r",
396
+ "arg": "integer",
397
+ "default": 1,
398
+ "description": "Number of replicates per fraction."
399
+ },
400
+ {
401
+ "name": "Out base",
402
+ "opt": "-o",
403
+ "arg": "out_file",
404
+ "description": ["Prefix of the output files to be created. The",
405
+ "output files will have a suffix of the form",
406
+ "'.fraction-replicate.fa', where 'fraction' is the percentage",
407
+ "sampled and 'replicate' is an increasing integer for replicates",
408
+ "of the same fraction. By default: Path to the input file."]
409
+ },
410
+ {
411
+ "name": "Force",
412
+ "opt": "-F",
413
+ "description": "Force overwriting output file(s)."
414
+ },
415
+ {
416
+ "name": "Zeroes",
417
+ "opt": "-z",
418
+ "description": ["Include leading zeroes in the numeric parts of the",
419
+ "output files (e.g., file.002.50-01.fa instead of file.2.50-1.fa),",
420
+ "so that alphabetic sorting of files reflects the sampled",
421
+ "fraction."]
422
+ },
423
+ {
424
+ "name": "Quiet",
425
+ "opt": "-q",
426
+ "description": "Run quietly."
427
+ },
428
+ {
429
+ "arg": "in_file",
430
+ "mandatory": true,
431
+ "multiple_sep": " ",
432
+ "description": "Input multi-FastA file(s)."
433
+ }
434
+ ]
435
+ },
436
+ {
437
+ "task": "FastA.tag.rb",
438
+ "description": "Generates easy-to-parse tagged reads from FastA files.",
439
+ "see_also": ["FastQ.tag.rb"],
440
+ "help_arg": "--help",
441
+ "options": [
442
+ {
443
+ "name": "Input file",
444
+ "opt": "--in",
445
+ "arg": "in_file",
446
+ "mandatory": true,
447
+ "description": "FastA file containing the sequences."
448
+ },
449
+ {
450
+ "name": "Output file",
451
+ "opt": "--out",
452
+ "arg": "out_file",
453
+ "mandatory": true,
454
+ "description": "FastA file to create."
455
+ },
456
+ {
457
+ "opt": "--prefix",
458
+ "arg": "string",
459
+ "description": "Prefix to use in all IDs."
460
+ },
461
+ {
462
+ "opt": "--suffix",
463
+ "arg": "string",
464
+ "description": "Suffix to use in all IDs."
465
+ },
466
+ {
467
+ "opt": "--defline",
468
+ "description": "Keep the original defline after a space."
469
+ },
470
+ {
471
+ "opt": "--quiet",
472
+ "description": "Run quietly (no STDERR output)."
473
+ }
474
+ ]
475
+ },
476
+ {
477
+ "task": "FastA.per_file.pl",
478
+ "description": ["Extracts all the sequences in a multi-FastA into",
479
+ "multiple single-FastA files."],
480
+ "see_also": "FastA.split.pl",
481
+ "help_arg": "",
482
+ "options": [
483
+ {
484
+ "arg": "out_dir",
485
+ "mandatory": true,
486
+ "description": "Output directory for the individual files."
487
+ },
488
+ {
489
+ "arg": "in_file",
490
+ "mandatory": true,
491
+ "multiple_sep": " ",
492
+ "description": "Input multi-FastA file(s)."
493
+ }
494
+ ]
495
+ },
496
+ {
497
+ "task": "FastA.fragment.rb",
498
+ "description": ["Simulates incomplete (fragmented) drafts from complete",
499
+ "genomes."],
500
+ "help_arg": "--help",
501
+ "options": [
502
+ {
503
+ "name": "Input file",
504
+ "opt": "--in",
505
+ "arg": "in_file",
506
+ "mandatory": true,
507
+ "description": "FastA file containing the complete sequences."
508
+ },
509
+ {
510
+ "name": "Output file",
511
+ "opt": "--out",
512
+ "arg": "out_file",
513
+ "mandatory": true,
514
+ "description": "FastA to create."
515
+ },
516
+ {
517
+ "opt": "--completeness",
518
+ "arg": "float",
519
+ "mandatory": true,
520
+ "description": ["Fraction of genome completeness to simulate from 0",
521
+ "to 1."]
522
+ },
523
+ {
524
+ "name": "Minimum length",
525
+ "opt": "--minlen",
526
+ "arg": "integer",
527
+ "description": "Minimum fragment length to report.",
528
+ "default": 500
529
+ },
530
+ {
531
+ "opt": "--sorted",
532
+ "description": ["Keep fragments sorted as in the input file. By",
533
+ "default, fragments are shuffled."]
534
+ },
535
+ {
536
+ "opt": "--quiet",
537
+ "description": "Run quietly (no STDERR output)."
538
+ }
539
+ ]
540
+ },
541
+ {
542
+ "task": "FastA.wrap.rb",
543
+ "description": "Wraps sequences in a FastA to a given line length.",
544
+ "help_arg": "--help",
545
+ "options": [
546
+ {
547
+ "name": "Input FastA",
548
+ "opt": "--in",
549
+ "arg": "in_file",
550
+ "mandatory": true,
551
+ "description": "Input FastA file."
552
+ },
553
+ {
554
+ "name": "Output FastA",
555
+ "opt": "--out",
556
+ "arg": "out_file",
557
+ "mandatory": true,
558
+ "description": "Output FastA file."
559
+ },
560
+ {
561
+ "name": "Wrap length",
562
+ "opt": "--wrap",
563
+ "arg": "integer",
564
+ "description": ["Line length to wrap sequences. Use 0 to generate",
565
+ "1-line sequences."],
566
+ "default": 70
567
+ }
568
+ ]
569
+ }
570
+ ]
571
+ }