When using assembly: "GCF_005577135.1", species name is not parsed correctly, which is causing failure in create_bsgenome.R and design_guides.R about snakemake-crispr-guides HOT 4 CLOSED

When checking, it also did not become obvious to me if the downloaded meta data is anyways of importance for downstream analyses etc.
If it is not, a possible solution might be to just simply give some mock species ID or species name and ignore the meta data.

from snakemake-crispr-guides.

Checked again and "Picosynechococcus sp. PCC 11901" does not seem to be a valid species name either, so "Synechococcus sp. PCC 11901" would have to be entered instead, which is not part of the fasta headers or given in the .gff file which were downloaded. A possible solution would be to extract the species field given in the gff file - even though I am not sure if all gffs contain this field (in case of Synechococcus 11901 the following line: ##species https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?id=2579791)

from snakemake-crispr-guides.

good find!
will try to reproduce this and see how to fix. Ideally there is an automatic solution such that the user does not need to specify the species name manually

from snakemake-crispr-guides.

@ute-hoffmann check out dev branch, where this is fixed now.

from snakemake-crispr-guides.