miga-base 1.2.17.0 → 1.2.17.2

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
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,291 @@
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 checking period to zero. If you",
49
+ "want to decrease the sampling period (to speed the script up) or",
50
+ "increase it (to make it more sensitive to errors) just change the",
51
+ "checking period 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.maskQual.rb",
86
+ "description": "Masks low-quality bases in a FastQ file.",
87
+ "help_arg": "--help",
88
+ "options": [
89
+ {
90
+ "opt": "--input",
91
+ "arg": "in_file",
92
+ "mandatory": true,
93
+ "description": ["Path to the FastQ file containing the sequences.",
94
+ "Supports compression with .gz extension."]
95
+ },
96
+ {
97
+ "opt": "--output",
98
+ "arg": "out_file",
99
+ "mandatory": true,
100
+ "description": ["Path to the output FastQ file.",
101
+ "Supports compression with .gz extension."]
102
+ },
103
+ {
104
+ "opt": "--qual",
105
+ "arg": "integer",
106
+ "default": 15,
107
+ "description": "Minimum quality score to allow a base."
108
+ },
109
+ {
110
+ "opt": "--offset",
111
+ "arg": "integer",
112
+ "default": 33,
113
+ "description": "Q-score offset."
114
+ },
115
+ {
116
+ "opt": "--fasta",
117
+ "description": "Output sequences in FastA format."
118
+ },
119
+ {
120
+ "opt": "--quiet",
121
+ "description": "Run quietly."
122
+ }
123
+ ]
124
+ },
125
+ {
126
+ "task": "FastQ.offset.pl",
127
+ "description": ["There are several FastQ formats. This script takes a",
128
+ "FastQ in any of them, identifies the type of FastQ (this is, the",
129
+ "offset), and generates a FastQ with the given offset."],
130
+ "warn": ["Note that Solexa+64 FastQ can cause problematic values when",
131
+ "using the offset 33, since there is no equivalent in Phred+33 for",
132
+ "negative values (the range of Solexa+64 is -5 to 40)."],
133
+ "help_arg": "",
134
+ "options": [
135
+ {
136
+ "arg": "in_file",
137
+ "mandatory": true,
138
+ "description": ["Input file in FastQ format (range is automatically",
139
+ "detected)."]
140
+ },
141
+ {
142
+ "name": "Offset",
143
+ "arg": "integer",
144
+ "default": 33,
145
+ "mandatory": true,
146
+ "description": ["Offset to use for the output. Use 0 (zero) to",
147
+ "detect the input format and exit."]
148
+ },
149
+ {
150
+ "opt": "force",
151
+ "description": ["If set, turns errors into warnings and continues.",
152
+ "Out-of-range values are set to the closest range limit."]
153
+ },
154
+ ">",
155
+ {
156
+ "arg": "out_file",
157
+ "mandatory": true,
158
+ "description": ["Output file in FastQ format with the specified",
159
+ "offset."]
160
+ }
161
+ ]
162
+ },
163
+ {
164
+ "task": "FastQ.split.pl",
165
+ "description": ["Splits a FastQ file into several FastQ files. This",
166
+ "script can be used to separate interposed sister reads using any even",
167
+ "number of output files."],
168
+ "help_arg": "",
169
+ "see_also": ["FastQ.interpose.pl","FastA.split.pl"],
170
+ "options": [
171
+ {
172
+ "name": "in_file.fq",
173
+ "arg": "in_file",
174
+ "mandatory": true,
175
+ "description": "Input file in FastQ format."
176
+ },
177
+ {
178
+ "name": "out_base",
179
+ "arg": "out_file",
180
+ "mandatory": true,
181
+ "description": ["Prefix for the name of the output files. It will be",
182
+ "appended with .<i>.fastq, where <i> is a consecutive number",
183
+ "starting in 1."]
184
+ },
185
+ {
186
+ "name": "no_files",
187
+ "arg": "integer",
188
+ "default": 2,
189
+ "description": "Number of files to generate."
190
+ }
191
+ ]
192
+ },
193
+ {
194
+ "task": "FastQ.tag.rb",
195
+ "description": "Generates easy-to-parse tagged reads from FastQ files.",
196
+ "see_also": ["FastA.tag.rb"],
197
+ "help_arg": "--help",
198
+ "options": [
199
+ {
200
+ "name": "Input file",
201
+ "opt": "--in",
202
+ "arg": "in_file",
203
+ "mandatory": true,
204
+ "description": [
205
+ "FastQ file containing the sequences.",
206
+ "Supports compression with .gz extension."
207
+ ]
208
+ },
209
+ {
210
+ "name": "Output file",
211
+ "opt": "--out",
212
+ "arg": "out_file",
213
+ "mandatory": true,
214
+ "description": [
215
+ "FastQ to create.",
216
+ "Supports compression with .gz extension."
217
+ ]
218
+ },
219
+ {
220
+ "opt": "--prefix",
221
+ "arg": "string",
222
+ "description": "Prefix to use in all IDs."
223
+ },
224
+ {
225
+ "opt": "--suffix",
226
+ "arg": "string",
227
+ "description": "Suffix to use in all IDs."
228
+ },
229
+ {
230
+ "opt": "--quiet",
231
+ "description": "Run quietly (no STDERR output)."
232
+ }
233
+ ]
234
+ },
235
+ {
236
+ "task": "FastQ.toFastA.awk",
237
+ "description": "Translates FastQ files into FastA.",
238
+ "see_also": "FastA.toFastQ.rb",
239
+ "help_arg": "'' --help",
240
+ "options": [
241
+ "<",
242
+ {
243
+ "arg": "in_file",
244
+ "mandatory": true,
245
+ "description": "Input FastQ file."
246
+ },
247
+ ">",
248
+ {
249
+ "arg": "out_file",
250
+ "mandatory": true,
251
+ "description": "Output FastA file."
252
+ }
253
+ ]
254
+ },
255
+ {
256
+ "task": "FastQ.test-error.rb",
257
+ "description": ["Compares the estimated error of sequencing reads",
258
+ "(Q-score) with observed mismatches (identity against a know",
259
+ "reference sequence)."],
260
+ "help_arg": "--help",
261
+ "options": [
262
+ {
263
+ "name": "FastQ",
264
+ "opt": "--fastq",
265
+ "arg": "in_file",
266
+ "mandatory": true,
267
+ "description": "FastQ file containing the sequences."
268
+ },
269
+ {
270
+ "name": "Tabular BLAST",
271
+ "opt": "--blast",
272
+ "arg": "in_file",
273
+ "mandatory": true,
274
+ "description": ["Tabular BLAST file mapping reads to reference",
275
+ "sequences."]
276
+ },
277
+ {
278
+ "name": "Output",
279
+ "opt": "--out",
280
+ "arg": "out_file",
281
+ "mandatory": true,
282
+ "description": "Output tab-delimited file to create."
283
+ },
284
+ {
285
+ "opt": "--quiet",
286
+ "description": "Run quietly (no STDERR output)."
287
+ }
288
+ ]
289
+ }
290
+ ]
291
+ }
@@ -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,165 @@
1
+ {
2
+ "tasks": [
3
+ {
4
+ "task": "BedGraph.tad.rb",
5
+ "description": [
6
+ "Estimates the truncated average sequencing depth (TAD)",
7
+ "from a BedGraph file."
8
+ ],
9
+ "warn": [
10
+ "This script doesn't consider zero-coverage positions if",
11
+ "missing from the file. If you produce your BedGraph file with",
12
+ "bedtools genomecov and want to consider zero-coverage position, be",
13
+ "sure to use -bga (not -bg)."
14
+ ],
15
+ "see_also": [
16
+ "BedGraph.window.rb", "BlastTab.seqdepth.pl", "BlastTab.seqdepth_ZIP.pl"
17
+ ],
18
+ "help_arg": "--help",
19
+ "options": [
20
+ {
21
+ "opt": "--input",
22
+ "arg": "in_file",
23
+ "mandatory": true,
24
+ "description": [
25
+ "Input BedGraph file.",
26
+ "Supports compression with .gz extension, use - for STDIN."
27
+ ]
28
+ },
29
+ {
30
+ "opt": "--output",
31
+ "arg": "out_file",
32
+ "default": "-",
33
+ "description": [
34
+ "Output tab-delimited values (by default, STDOUT).",
35
+ "Supports compression with .gz extension, use - for STDOUT."
36
+ ]
37
+ },
38
+ {
39
+ "opt": "--range",
40
+ "arg": "float",
41
+ "default": 0.5,
42
+ "description": [
43
+ "Central range to consider, between 0 and 1.",
44
+ "By default: inter-quartile range (0.5)."
45
+ ]
46
+ },
47
+ {
48
+ "opt": "--name",
49
+ "arg": "string",
50
+ "description": [
51
+ "Name (ID) of the sequence, added as first column."
52
+ ]
53
+ },
54
+ {
55
+ "opt": "--per-seq",
56
+ "description": [
57
+ "Calculate averages per reference sequence, not total.",
58
+ "Assumes a sorted BedGraph file."
59
+ ]
60
+ },
61
+ {
62
+ "opt": "--length",
63
+ "description": "Add sequence length to the output."
64
+ }
65
+ ]
66
+ },
67
+ {
68
+ "task": "BedGraph.window.rb",
69
+ "description": ["Estimates the sequencing depth per windows from a",
70
+ "BedGraph file."],
71
+ "warn": ["This script doesn't consider zero-coverage positions if",
72
+ "missing from the file. If you produce your BedGraph file with",
73
+ "bedtools genomecov and want to consider zero-coverage position, be",
74
+ "sure to use -bga (not -bg)."],
75
+ "see_also": ["BedGraph.tad.rb",
76
+ "BlastTab.seqdepth.pl", "BlastTab.seqdepth_ZIP.pl"],
77
+ "help_arg": "--help",
78
+ "options": [
79
+ {
80
+ "opt": "--input",
81
+ "arg": "in_file",
82
+ "mandatory": true,
83
+ "description": "Input BedGraph file."
84
+ },
85
+ {
86
+ "name": "Window size",
87
+ "opt": "--win",
88
+ "arg": "float",
89
+ "default": 1000,
90
+ "description": "Window size, in base pairs."
91
+ }
92
+ ]
93
+ },
94
+ {
95
+ "task": "sam.filter.rb",
96
+ "description": ["Filters a SAM or BAM file by target sequences and/or",
97
+ "identity."],
98
+ "see_also": ["anir.rb"],
99
+ "help_arg": "--help",
100
+ "options": [
101
+ {
102
+ "opt": "--genome",
103
+ "arg": "in_file",
104
+ "mandatory": true,
105
+ "description": ["Genome assembly.",
106
+ "Supports compression with .gz extension."]
107
+ },
108
+ {
109
+ "opt": "--mapping",
110
+ "arg": "in_file",
111
+ "mandatory": true,
112
+ "description": ["Mapping file.",
113
+ "Supports compression with .gz extension."]
114
+ },
115
+ {
116
+ "opt": "--out-sam",
117
+ "arg": "out_file",
118
+ "mandatory": true,
119
+ "description": ["Output filtered file in SAM format.",
120
+ "Supports compression with .gz extension."]
121
+ },
122
+ {
123
+ "opt": "--g-format",
124
+ "arg": "select",
125
+ "values": ["fasta", "list"],
126
+ "default": "fasta",
127
+ "description": ["Genome assembly format."]
128
+ },
129
+ {
130
+ "opt": "--m-format",
131
+ "arg": "select",
132
+ "values": ["sam", "bam"],
133
+ "default": "sam",
134
+ "description": ["Mapping file format. SAM supports compression with",
135
+ ".gz file extension."]
136
+ },
137
+ {
138
+ "opt": "--identity",
139
+ "arg": "float",
140
+ "description": "Set a fixed threshold of percent identity.",
141
+ "default": 95.0
142
+ },
143
+ {
144
+ "opt": "--no-header",
145
+ "description": "Do not include the headers."
146
+ },
147
+ {
148
+ "opt": "--threads",
149
+ "arg": "integer",
150
+ "description": "Threads to use.",
151
+ "default": 2
152
+ },
153
+ {
154
+ "opt": "--log",
155
+ "arg": "out_file",
156
+ "description": "Log file to save output."
157
+ },
158
+ {
159
+ "opt": "--quiet",
160
+ "description": "Run quietly."
161
+ }
162
+ ]
163
+ }
164
+ ]
165
+ }