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

Delimeters for Genomic Addres #34

Open
sgsutcliffe opened this issue Jan 8, 2025 · 0 comments
Open

Delimeters for Genomic Addres #34

sgsutcliffe opened this issue Jan 8, 2025 · 0 comments
Labels
bug Something isn't working

Comments

@sgsutcliffe
Copy link
Contributor

Description of the bug

Right now, gasnomenclature does not work as expected, or has some issues if a delimeter for the genomic address besides the default ".".

GAS call is supposed to be flexible in the output, in that a delimeter other than the default can be set by using --gm_delimiter but the expected behavior of a command like:

nextflow run main.nf -profile docker,test --outdir results --gm_delimiter -

Still gives an output with the delimeter as the default "." which might be a GAS issue.

In addition, the tool is supposed to be able to accept address in the samplesheet with none "." characters but this is currently not the case. For example a samplesheet using:

sample,mlst_alleles,address
sampleQ,https://raw.githubusercontent.com/phac-nml/gasnomenclature/dev/tests/data/reports/sampleQ.mlst.json,
sample1,https://raw.githubusercontent.com/phac-nml/gasnomenclature/dev/tests/data/reports/sample1.mlst.json,1-1-1
sample2,https://raw.githubusercontent.com/phac-nml/gasnomenclature/dev/tests/data/reports/sample2.mlst.json,1-1-1
sample3,https://raw.githubusercontent.com/phac-nml/gasnomenclature/dev/tests/data/reports/sample3.mlst.json,1-1-2

Causes an error. Though it was due to the schema for address not including anything besides "."

"pattern": "^\\d+(\\.\\d+)*$",

However when I changed the pattern to "pattern": "^\\d+([_.-]\\d+)*$", it still fails. I think this is due to an issue with GAS 0.1.3
as the error looks like:

Error something went wrong with cluster assignment. check error messages ["Genomic address too short for samples: ['sample1', 'sample2', 'sample3'] based on {0: 10.0, 1: 5.0, 2: 0.0}"]

And when I run GAS with the file files I get the same thing.

In the PR 31 we also made it so the addresses in the database and samplesheet must have the same delimeter output by GAS but that does not break anything because right now gasnomenclature can only handle the delimeter "."

Command used and terminal output

No response

Relevant files

No response

System information

No response

@sgsutcliffe sgsutcliffe added the bug Something isn't working label Jan 8, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant