Skip to content

Conversation

@j23414
Copy link
Contributor

@j23414 j23414 commented Aug 18, 2025

Description of proposed changes

Similar to the fix in the Nextclade workflow

  1. Align the nextclade genome reference and the phylogenetic genome reference
  2. Compare P and I CDS and gene annotations
  3. Glance at the annotated translations and the automatic translations

Basically changing from original:

Screenshot 2025-08-18 at 11 32 08 AM

to:

Screenshot 2025-08-18 at 11 37 04 AM

Related issue(s)

Checklist

  • Checks pass
  • Update changelog

@j23414 j23414 marked this pull request as ready for review August 18, 2025 19:00
@j23414
Copy link
Contributor Author

j23414 commented Oct 7, 2025

Handles slippage in P gene related to nextstrain/augur#1901 (comment)

Comment on lines 83 to 85
gene 1916..2426
/gene="I"
CDS 1916..2431
CDS join(1916..2380,2377..2426)
Copy link
Member

Choose a reason for hiding this comment

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

Tested with nextstrain/augur#1901. P gene looks good now, but there's a new error with I gene:

$ nextstrain build --augur=/path/to/local/nextstrain/augur phylogenetic results/global/aa_muts.json --rerun-incomplete
ERROR: Reference file 'defaults/reference.gb' failed validation:
    'I' has length 515 which is not a multiple of 3.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Thanks, fixup in 4018a3c based on what I can see in Geneious and genbank annotation

Screenshot 2025-10-10 at 12 06 46 PM

Copy link
Member

Choose a reason for hiding this comment

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

I'm not familiar with the reasoning, but augur translate now succeeds in the test command so this looks good to me!

Similar to the fix in the Nextclade workflow.

* https://github.com/nextstrain/mumps/blame/41654ec8e601f4669146d4c3ac1834ae7f74c12f/nextclade/defaults/genome/reference.gb#L54

1. Align the nextclade genome reference and the phylogenetic genome reference
2. Compare P and I CDS and gene annotations
3. Glance at the annotated translations and the automatic translations
@j23414 j23414 force-pushed the fix-p-gene-annotation branch from 4018a3c to 573491f Compare October 10, 2025 20:23
@j23414 j23414 merged commit 61ccca3 into main Oct 10, 2025
4 checks passed
@j23414 j23414 deleted the fix-p-gene-annotation branch October 10, 2025 20:29
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants