enable protein sketches for pangenome hash correlations

[AnneliektH]

calc-hash-presence.py from pangenome-hash-cor, does not yet allow me to compare protein sketches instead of nucleotide sketches.

I'd like to compare protein-pangenomes instead of nucleotide ones

Trying in /group/ctbrowngrp2/scratch/annie/2023-swine-sra/sourmash/pangenomics/test_virpan

python ../../2024-pangenome-hash-corr/calc-hash-presence.py \
cluster1371.ranktable.csv \
cluster1371.zip \
-o cluster1371.dump \
--scaled=1

loaded 204371 hashvals... downsampling soon.
found 0 metagenomes
Traceback (most recent call last):
  File "/group/ctbrowngrp2/scratch/annie/2023-swine-sra/sourmash/pangenomics/test_virpan/sigs/../../2024-pangenome-hash-corr/calc-hash-presence.py", line 80, in <module>
    sys.exit(main())
  File "/group/ctbrowngrp2/scratch/annie/2023-swine-sra/sourmash/pangenomics/test_virpan/sigs/../../2024-pangenome-hash-corr/calc-hash-presence.py", line 37, in main
    query_minhash = next(iter(idx.signatures())).minhash.copy_and_clear()
StopIteration

[ctb]

just added it now - ctb/2024-pangenome-hash-corr@7440f8f

you'll need to update your sourmash_plugin_pangenomics first, either to the latest released version or the latest dev version.

[ctb]

(let me know if it breaks ;)

[ctb]

hmm, some more things need to be fixed before it will work. sorry!

[ctb]

ok fixed in 2024-pangenome-hash-corr ctb/2024-pangenome-hash-corr@f2433f2 and v0.2.2 of sourmash_plugin_pangenomics, I think.

you may need to use --protein --no-dna with the pangenome plugin scripts. You'll also need to specify the k-mer size and scaled value far too many times. Sorry, work in progress ;).