aboutsummaryrefslogtreecommitdiff
path: root/workflows/phylogeny/README.md
diff options
context:
space:
mode:
authorPjotr Prins2021-01-26 16:22:19 +0000
committerGitHub2021-01-26 16:22:19 +0000
commita54122023ef2bde4ef4fc39355b3491678bbf773 (patch)
tree150bef4a574486db9e7c16d87be89ad219d84bb1 /workflows/phylogeny/README.md
parentb112b3203e34ea61dfdf802bce5036f938eaa774 (diff)
parent4842117f530aba08b4a253aee533a8d4802b1c12 (diff)
downloadbh20-seq-resource-a54122023ef2bde4ef4fc39355b3491678bbf773.tar.gz
bh20-seq-resource-a54122023ef2bde4ef4fc39355b3491678bbf773.tar.lz
bh20-seq-resource-a54122023ef2bde4ef4fc39355b3491678bbf773.zip
Merge pull request #124 from urbanslug/master
Add phylogeny workflow
Diffstat (limited to 'workflows/phylogeny/README.md')
-rw-r--r--workflows/phylogeny/README.md38
1 files changed, 38 insertions, 0 deletions
diff --git a/workflows/phylogeny/README.md b/workflows/phylogeny/README.md
new file mode 100644
index 0000000..6ecc047
--- /dev/null
+++ b/workflows/phylogeny/README.md
@@ -0,0 +1,38 @@
+A workflow to generate a phylogeny that can be visualized using [auspice](https://github.com/urbanslug/auspice).
+Expects a multi-fasta file path at [pggb_fasta][1] and generates a tree in `json` format.
+
+#### Dependencies
+
+Depends on:
+ - [pggb](https://github.com/pangenome/pggb/blob/master/pggb)
+ * [wfmash](https://github.com/ekg/wfmash)
+ * [seqwish](https://github.com/ekg/seqwish)
+ * [smoothxg](https://github.com/pangenome/smoothxg)
+ * [odgi](https://github.com/vgteam/odgi)
+
+ - [taxophages](https://github.com/urbanslug/taxophages/)
+ * Clone and run with `python main.py ...`
+
+ - [augur](https://github.com/nextstrain/augur)
+
+
+#### Running
+
+Expects that taxophages is cloned in a previous dir but you can update the path [main_py_script][2] to wherever it is.
+
+Run the phylogeny workflow with the bleow after specifying your path to [pggb_fasta][1].
+
+```bash
+R_PACKAGES="${HOME}/RLibraries" \ # a directory holding R packages. Needed if R packages installed using install.packages on server e.g https://github.com/urbanslug/taxophages/blob/master/scripts/deps.R
+TAXOPHAGES_ENV=server \ # helps taxophages figure out where it is being ran
+AUGUR_RECURSION_LIMIT=30000 \ # augur isn't used to working with so many nested values
+cwltool --preserve-entire-environment --no-container phylogeny.cwl clado-job.yml
+```
+
+Alternatively run any workflow with
+```
+cwltool --no-container <workflow>.cwl clado-job.yml
+```
+
+[1]: clado-job.yml#L8
+[2]: clado-job.yml#L28