Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

sketch,build sbts w sourmash 4.0, multiple alphabets #15

Closed
wants to merge 6 commits into from

Conversation

bluegenes
Copy link
Contributor

@bluegenes bluegenes commented Oct 22, 2020

Here's a Snakefile for building gtdb (or other) sbts using sourmash 4.0 sketching & indexing. Also added an environment.yml with environment than can be used to run this snakefile.

Configfile requires (default index_config.yml):

  • genome and/or protein csv with at minimum, accession, filename columns
    • (provide signame for better signature naming)
  • data-dir: where to look for the genome/proteome files
  • specify only desired alphabet-ksize-scaled combinations for building sbts.

To run on different data or with different parameters, modify index_config.yml or provide a different configfile.

Alphabet defaults are taken from sourmash-bio/sourmash#1168.
Workflow does not generate translated sketches, as we don't want to build databases from 6-frame translations.

I also modified the farm sbatch files, and have corresponding snakemake profiles to use, but not sure we want to add those to this repo.

Likely needs modification (e.g. batching?) for even more genomes or additional parameters (bloom filter size?), but here's a starting point for 4.x.

@bluegenes
Copy link
Contributor Author

Closing in favor of #16

@bluegenes bluegenes closed this May 11, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant