Skip to content

Commit

Permalink
very preliminary credits and tutorial
Browse files Browse the repository at this point in the history
  • Loading branch information
Katherine Eaton authored and ktmeaton committed Apr 27, 2022
1 parent e9c41e6 commit d4e0aa8
Showing 1 changed file with 16 additions and 22 deletions.
38 changes: 16 additions & 22 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -54,48 +54,35 @@ SARS-CoV-2 recombinant sequence detection inspired by [nextstrain/ncov](https://
conda activate ncov-recombinant
```

## Test

- Testing uses the "controls" dataset `data/controls`.
- These are publicly available recombinant sequences from [pango-designation](https://github.com/cov-lineages/pango-designation) issues.
- Testing uses the system profile `profiles/laptop`.
- This specifies system resources of 1 cpus and 4 GB of memory.

> Tip: Remember to run `conda activate ncov-recombinant` before running the commands below!

1. Preview the run configuration.

```bash
snakemake --profile profiles/laptop print_config
```
## Tutorial

1. Preview the commands that are going to be run.

```bash
snakemake --profile profiles/laptop --dryrun --printshellcmds
snakemake --profile profiles/tutorial --dryrun --printshellcmds
```

1. Run the workflow.

```bash
snakemake --profile profiles/laptop
snakemake --profile profiles/tutorial
```

1. Inspect the output

- Sample line list:

```bash
csvtk pretty -t results/controls/summary.tsv
csvtk pretty -t results/tutorial/linelist.tsv
```

- Breakpoints: `results/controls/sc2rf.txt`
- Breakpoints:

```bash
less -S results/controls/sc2rf.txt
less -S results/tutorial/sc2rf.recombinants.ansi.txt
```

- Trees: PNG and Auspice JSON in `results/controls/subtrees_collapse`
- Trees: Aupoad JSON in `results/tutorials/subtrees_collapse/` to <https://auspice.us/>

## SLURM

Expand All @@ -110,8 +97,15 @@ bash scripts/slurm.sh --profile profiles/hpc --partition MyPartition
- Unlock the workflow after a failed run:

```bash
snakemake --profile profiles/laptop --unlock
snakemake --profile profiles/tutorial --unlock
```

- Logs for each rule will be stored at: `logs/<rule>/<build>_<date>.txt`
- If using slurm (`profiles/hpc`) the master log will be at: `logs/slurm/ncov-recombinant_<date>_<jobid>.log`
- If using slurm (`profiles/hpc`) the master log will be at: `logs/ncov-recombinant/ncov-recombinant_<date>_<jobid>.log`

## Credits

- nextstrain/ncov: https://github.com/nextstrain/ncov
- nextstrain/nextclade: https://github.com/nextstrain/nextclade
- lenaschimmel/sc2rf: https://github.com/lenaschimmel/sc2rf
- yatisht/usher: https://github.com/yatisht/usher

0 comments on commit d4e0aa8

Please sign in to comment.