|
1 | 1 | <img src="https://github.com/sunbeam-labs/sunbeam/blob/stable/docs/images/sunbeam_logo.gif" width=120, height=120 align="left" /> |
2 | 2 |
|
3 | | -# sbx_virus_id |
| 3 | +# sbx_cenote_taker |
4 | 4 |
|
5 | 5 | <!-- badges: start --> |
6 | | -[](https://github.com/sunbeam-labs/sbx_virus_id/actions/workflows/tests.yml) |
7 | | -[](https://hub.docker.com/repository/docker/sunbeamlabs/sbx_virus_id/) |
| 6 | +[](https://github.com/sunbeam-labs/sbx_cenote_taker/actions/workflows/tests.yml) |
| 7 | +[](https://hub.docker.com/repository/docker/sunbeamlabs/sbx_cenote_taker/) |
8 | 8 | <!-- badges: end --> |
9 | 9 |
|
10 | 10 | ## Introduction |
11 | 11 |
|
12 | | -sbx_virus_id is a [sunbeam](https://github.com/sunbeam-labs/sunbeam) extension for identifying viruses in samples. This pipeline uses [MEGAHIT](https://github.com/voutcn/megahit) or [SPAdes](https://github.com/ablab/spades) for assembly of contigs and [Cenote-Taker2](https://github.com/mtisza1/Cenote-Taker2) or [Virsorter2](https://github.com/jiarong/VirSorter2) for viral identification. |
| 12 | +sbx_cenote_taker is a [sunbeam](https://github.com/sunbeam-labs/sunbeam) extension for identifying viruses in samples with [Cenote-Taker3](https://github.com/jedvachey/Cenote-Taker3). This pipeline uses [MEGAHIT](https://github.com/voutcn/megahit) for assembly of contigs and then processes assemblies with Cenote-Taker3. |
13 | 13 |
|
14 | | -N.B. If using Megahit for assembly, this extension requires also having sbx_assembly installed. |
| 14 | +N.B. This extension requires also having sbx_assembly installed. |
15 | 15 |
|
16 | 16 | ### Installation |
17 | 17 |
|
18 | 18 | ``` |
19 | | -sunbeam extend https://github.com/sunbeam-labs/sbx_virus_id.git |
| 19 | +sunbeam extend https://github.com/sunbeam-labs/sbx_assembly.git |
| 20 | +sunbeam extend https://github.com/sunbeam-labs/sbx_cenote_taker.git |
20 | 21 | ``` |
21 | 22 |
|
22 | | -# Installing blast dbs |
| 23 | +### Cenote-Taker database |
23 | 24 |
|
24 | | -Install blast db: |
| 25 | +sbx_cenote_taker expects the Cenote-Taker3 reference database to be available locally. Download the database following the official instructions, for example: |
25 | 26 |
|
26 | 27 | ``` |
27 | | -conda create -n blast |
28 | | -conda activate blast |
29 | | -conda install -c bioconda blast |
30 | | -mkdir refseq_select_prot/ |
31 | | -cd refseq_select_prot/ |
32 | | -perl `which update_blastdb.pl` --decompress refseq_select_prot |
| 28 | +conda activate cenote-taker |
| 29 | +get_ct3_dbs -o /path/to/ct3_db --hmm T --hallmark_tax T --refseq_tax T --mmseqs_cdd T --domain_list T --hhCDD T --hhPFAM T --hhPDB T |
33 | 30 | ``` |
34 | 31 |
|
35 | | -Install viral blast db: |
| 32 | +Update the `cenote_taker_db` entry in your Sunbeam configuration to point at the resulting directory. |
36 | 33 |
|
37 | | -``` |
38 | | -conda stuff from above ^^^ |
39 | | -mkdir viral_prot/ && cd viral_prot/ |
40 | | -wget https://ftp.ncbi.nlm.nih.gov/refseq/release/viral/viral.1.protein.faa.gz && gzip -d viral.1.protein.faa.gz |
41 | | -makeblastdb -in viral.1.protein.faa -parse_seqids -title "viral" -dbtype prot |
42 | | -``` |
| 34 | +### Running |
43 | 35 |
|
44 | | -## Running |
45 | | - |
46 | | -Run with sunbeam on the target `all_virus_id`, |
| 36 | +Run with sunbeam on the target `all_cenote_taker`: |
47 | 37 |
|
48 | 38 | ``` |
49 | | -sunbeam run --profile /path/to/project/ all_virus_id |
| 39 | +sunbeam run --profile /path/to/project/ all_cenote_taker |
50 | 40 | ``` |
51 | 41 |
|
52 | 42 | ### Options for config.yml |
53 | 43 |
|
54 | 44 | - blast_db: path to blast db (default: "") (NOTE: this should be the database file not just the directory it's in) |
55 | 45 | - blastx_threads: number of threads for running blastx (default: 4) |
56 | 46 | - bowtie2_build_threads: number of threads for running bowtie2-build (default: 4) |
57 | | - - cenote_taker2_db: path to cenote-taker2 db (default: "") (NOTE: this should be a directory) |
58 | | - - virsorter_db: path to virsorter2 db (default: "") (NOTE: this should be a directory) |
| 47 | + - cenote_taker_db: path to cenote-taker3 db (default: "") (NOTE: this should be a directory) |
59 | 48 | - include_phages: Whether to include phages in the output (default: False) |
60 | | - - use_spades: Whether to use SPAdes instead of MEGAHIT (default: False) |
61 | | - - use_virsorter: Whether to use Virsorter2 instead of Cenote-Taker2 (default: False) |
62 | | - |
63 | | -## Legacy Installation |
64 | | - |
65 | | -``` |
66 | | -git clone https://github.com/sunbeam-labs/sbx_virus_id.git extensions/sbx_virus_id |
67 | | -cd extensions/sbx_virus_id |
68 | | -cat config.yml >> /path/to/sunbeam_config.yml |
69 | | -``` |
0 commit comments