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

hundred XM genomes called as 'Unassigned' or as BA.1 #46

Closed
akifoss opened this issue May 3, 2022 · 5 comments
Closed

hundred XM genomes called as 'Unassigned' or as BA.1 #46

akifoss opened this issue May 3, 2022 · 5 comments

Comments

@akifoss
Copy link

akifoss commented May 3, 2022

Using the mutation profile mentioned in the pango-designation issue that assigned lineage XM (C15240T C2470T A18163G C19955T A20055G), the attached list of genomes did not get an XM lineage designation although they seem to belong to XM when looking at them with sc2rf. Would model training improve the lineage assignments or is this not in place anymore due to the upgrade to pangolin version 4?

Here's a sc2rf screenshot:
image

Many thanks!

XM-notCalled-2022-05-02-GISAID.csv

@corneliusroemer
Copy link

corneliusroemer commented May 3, 2022

Thanks for reporting this. I ran the sequences through Nextclade and Usher and indeed, most of these are good XM.

Only dozen or so are close but not quite XM (according to Usher).

https://nextstrain.org/fetch/genome.ucsc.edu/trash/ct/subtreeAuspice1_genome_ba86_149cd0.json?c=pango_lineage_usher

image

image

@corneliusroemer
Copy link

In particular all the German sequences are very likely real XM (XM is mostly German)
image

@JosetteSchoenma
Copy link

JosetteSchoenma commented May 3, 2022

I have noticed the same thing last week. Only 4 out of 11 Dutch samples that were called XM by Usher were called XM by GISAID/pangolearn. They were either unassigned or n/a. The ones that were assigned correctly are from just after designation.

@JosetteSchoenma
Copy link

Screenshot_20220503-230816_Chrome.jpg
And indeed, also in the world, most of these are unassigned.

@SVN-PhD
Copy link

SVN-PhD commented May 4, 2022

I ran your sequences through pangolin (v4.0.6) on the command line with pangolin-data v1.8. I think this is another scorpio overwriting Usher calls issue.

Screenshot from 2022-05-04 07-35-56

@corneliusroemer corneliusroemer transferred this issue from cov-lineages/pango-designation Aug 21, 2022
@rmcolq rmcolq closed this as completed Jun 7, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

5 participants