Skip to content

Instantly share code, notes, and snippets.

@jungbluth
Last active September 20, 2023 03:12
Show Gist options
  • Save jungbluth/ffbf49b8fcf3a17b7344f44761a198f7 to your computer and use it in GitHub Desktop.
Save jungbluth/ffbf49b8fcf3a17b7344f44761a198f7 to your computer and use it in GitHub Desktop.
run_sourmash.sh
#!/usr/bin/env bash
for line in $(cat list); do
FILE=/global/cfs/cdirs/img/web-data/sandbox.blast.data/${line}/${line}.a.fna
if [ -f "$FILE" ]; then
mamba activate sourmash
sourmash sketch dna -p k=31,scaled=1000 -o ${line}.fna.sig "$FILE"
fi
done
sourmash gather 3300009984.fna.sig gtdb-rs207.genomic-reps.dna.k31.zip --save-matches matches.zip
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment