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

Adding CGE-Predicted Phenotypes to staramr's Resfinder Output #167

Merged
merged 10 commits into from
Mar 14, 2023

Conversation

emarinier
Copy link
Member

@emarinier emarinier commented Mar 2, 2023

Adds CGE-predicted phenotypes to the output of staramr.

Example running with the README's tutorial data:

staramr search --pointfinder-organism salmonella -o out GCF_001478105.1.fasta

Example output:

Isolate ID	Gene	Predicted Phenotype	CGE Predicted Phenotype	%Identity	%Overlap	HSP Length/Total Length	Contig	Start	End	Accession	Sequence
GCF_001478105.1	blaCMY-2	ampicillin, amoxicillin/clavulanic acid, cefoxitin, ceftriaxone	Amoxicillin,  Amoxicillin+Clavulanic acid,  Ampicillin,  Ampicillin+Clavulanic acid,  Cefotaxime,  Cefoxitin,  Ceftazidime,  Piperacillin,  Piperacillin+Tazobactam,  Ticarcillin,  Ticarcillin+Clavulanic acid	100.00	100.00	1146/1146	ref|NZ_JYVD01000056.1|	25020	26165	X91840	ATGATGAAAAAATCGTTA[...]

@emarinier emarinier self-assigned this Mar 2, 2023
@emarinier emarinier marked this pull request as draft March 2, 2023 16:41
Copy link
Member

@apetkau apetkau left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What you have so far looks awesome. Great job getting this up and working so quickly 😄

The one comment I have so far is that you are storing a copy of the CGE gene drug table within the StarAMR data. Instead of doing this, I think it would be good to directly use the file installed as part of the ResFinder database (so that switching the ResFinder database used by StarAMR will also switch the gene-drug mapping to be consistent).

@emarinier
Copy link
Member Author

What you have so far looks awesome. Great job getting this up and working so quickly 😄

The one comment I have so far is that you are storing a copy of the CGE gene drug table within the StarAMR data. Instead of doing this, I think it would be good to directly use the file installed as part of the ResFinder database (so that switching the ResFinder database used by StarAMR will also switch the gene-drug mapping to be consistent).

Awesome. Okay, I'll look into using the installed version of the database.

@emarinier emarinier changed the title WIP: Adding CGE-predicted phenotypes to staramr's resfinder output Adding CGE-Predicted Phenotypes to staramr's Resfinder Output Mar 7, 2023
@emarinier emarinier marked this pull request as ready for review March 7, 2023 17:31
Copy link
Member

@apetkau apetkau left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks great. Thanks so much. Merging this in.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
No open projects
Status: Done
Development

Successfully merging this pull request may close these issues.

2 participants