|
1 | | -This folder is meant to contain all resources necessary for running the workflow, for example reference sequences or databases. |
| 1 | +# Resources |
2 | 2 |
|
3 | | -Wherever feasible, they can also be downloaded programmatically via rules defined in the pipeline. |
4 | | -Please create a directory, or multiple directories, here for your library, or libraries, that contains your raw fastq files. |
| 3 | +## Purpose of this directory |
5 | 4 |
|
6 | | -Each database should be added to its respective folder in the database directory. |
| 5 | +This directory contains all external reference data and input resources required to run MMCAW, including sequencing inputs, reference genomes, and taxonomic databases. Correct placement and formatting of these resources is essential for reproducible execution of the workflow. |
| 6 | + |
| 7 | +## Sequencing inputs |
| 8 | + |
| 9 | +### FASTQ files |
| 10 | +- **Description:** Basecalled and demultiplexed ONT reads generated using Guppy v6.0.6 (high-accuracy model). |
| 11 | +- **Format:** `.fastq.gz` |
| 12 | +- **Expected location:** User-defined; specified in the metadata file referenced in `config/config.yaml`. |
| 13 | + |
| 14 | +### Guppy sequence summaries |
| 15 | +- **Description:** Guppy summary files corresponding to each sequencing run. |
| 16 | +- **Format:** `.txt` |
| 17 | +- **Usage:** Used for run-level QC and tracking. |
| 18 | + |
| 19 | +### Unblocked read ID lists |
| 20 | +- **Description:** Lists of read IDs used to exclude blocked reads from downstream analysis. |
| 21 | +- **Format:** Plain text files. |
| 22 | + |
| 23 | +## Reference data |
| 24 | + |
| 25 | +### Human reference genome (GRCh38.p14) |
| 26 | +- **Description:** Used for optional host read filtering. |
| 27 | +- **Default path in config:** resources/databases/human_reference/GCF_000001405.40_GRCh38.p14_genomic.fna |
| 28 | + |
| 29 | +## Taxonomic databases |
| 30 | + |
| 31 | +### Kraken2 database |
| 32 | +- **Description:** Kraken2 standard database including bacterial, archaeal, viral, and human genomes. |
| 33 | +- **Default path in config:** ~/Kraken2_Simple_Workflow/resources/databases/krakenstd_06_2023/kraken2_std_database |
| 34 | + |
| 35 | +### BLAST (NT) database |
| 36 | +- **Description:** NCBI nucleotide (NT) database used for BLAST-based taxonomic assignment. |
| 37 | +- **Default path in config:** resources/databases/NCBI_blast_database/nt |
| 38 | + |
| 39 | +### CAT databases |
| 40 | +- **Description:** NCBI taxonomy and protein databases used by CAT for contig classification. |
| 41 | +- **Default paths in config:** resources/databases/20240422_CAT_nr/db, resources/databases/20240422_CAT_nr/tax |
| 42 | + |
| 43 | +### NCBI Taxonomy (taxdump) |
| 44 | +- **Description:** Required for taxonomy name resolution and LCA assignment. |
| 45 | +- **Default path in config:** resources/databases/taxdump |
| 46 | + |
| 47 | +## Optional resources |
| 48 | + |
| 49 | +Additional databases may be generated using the included database creation subworkflow by setting: include_db_creation: True in `config/config.yaml`. |
| 50 | + |
| 51 | +## How to obtain and prepare resources |
| 52 | + |
| 53 | +- Human reference genome: Download from NCBI (GRCh38.p14). |
| 54 | +- Kraken2 database: Build using Kraken2 standard database scripts or provide an existing installation. |
| 55 | +- BLAST NT database: Download via NCBI BLAST+ `update_blastdb.pl`. |
| 56 | +- CAT databases: Download from the official CAT repository or build using provided scripts. |
| 57 | + |
| 58 | +Ensure all paths in `config/config.yaml` correctly point to the corresponding files. |
| 59 | + |
| 60 | +## Directory structure (recommended) |
| 61 | + |
| 62 | +resources/ |
| 63 | + |
| 64 | +└── databases/ |
| 65 | + |
| 66 | +├── human_reference/ |
| 67 | + |
| 68 | +├── krakenstd_06_2023/ |
| 69 | + |
| 70 | +├── 20240422_CAT_nr/ |
| 71 | + |
| 72 | +├── NCBI_blast_database/ |
| 73 | + |
| 74 | +└── taxdump/ |
0 commit comments