scottdaniel / sbx_gene_clusters

sbx_gene_clusters

Reads-level based alignment to gene clusters of interest, e.g. bai operon or butyrate producing genes. Please refer to sunbeam_database for details.

Installation

To install, activate your conda environment (using the name of your environment) and use sunbeam extend:

conda activate <i>sunbeamX.X.X</i>
sunbeam extend https://github.com/sunbeam-labs/sbx_gene_clusters.git

Now take UniRef50 database as an example. First download the uniref50.fasta into your current sunbeam_output/mapping/sbx_gene_family/databases/.

mkdir -p sunbeam_output/mapping/sbx_gene_family/database/
wget ftp://ftp.uniprot.org/pub/databases/uniprot/uniref/uniref50/uniref50.fasta.gz -P sunbeam_output/mapping/sbx_gene_family/database/

Be sure to update the config.yml with the proper path.

Usage

To generate alignments,

sunbeam run --profile /path/to/project all_gene_family

Configuration

• threads: Is the number of threads to run parallel processes with
• genes_fp: Is the path to the downloaded database
• evalue:
• alnLen:
• mismatch:

Legacy Installation

For sunbeam versions <3 or if sunbeam extend isn't working, you can use git directly to install an extension:

git clone https://github.com/sunbeam-labs/sbx_gene_clusters.git extensions/sbx_gene_clusters

and then include it in the config for any given project with:

cat extensions/sbx_gene_clusters/config.yml >> /path/to/project/sunbeam_config.yml