ITADN

is there a way to give protein sites non-sequential numbering?

#1968Openjbloom 创建于 2026-03-02
enhancement
J
jbloomcommented
### Context In Nextstrain trees, is there a way to give protein sites non-sequential numbering for things like labeling branches and coloring by protein site? Right now I understand how you can pass a GFF to say define HA1 and HA2 to start at some position within the full coding sequence. But people often use H3 numbering for other subtypes (like H7 and H5). Here the numbers of sites in the H7 or H5 HA are not always strictly sequential: if H7 has a deletion relative to H3 then the H3 numbering skips a site, and if H7 has an indel then there is a duplicate number (eg, site 133a). PDB files are often numbered this way. Is there any established way to do this within Nextstrain trees? ### Description As an example, we currently use [H3 numbering](https://dms-vep.org/Flu_H5_American-Wigeon_South-Carolina_2021-H5N1_DMS/numbering.html) for mutations in a H5 HA. This mapping of numbering is shown [here](https://github.com/dms-vep/Flu_H5_American-Wigeon_South-Carolina_2021-H5N1_DMS/blob/main/data/site_numbering_map.csv) (see *reference_site* and *sequential_site* columns). ### Possible solution Perhaps one possible solution is to define a very complex GFF that defines a non-sequential feature, although this could not capture all forms of indels between the numbering schemes. Not sure if @lmoncla @kistlerk @Haddox @jordan-ort have dealt with this?
0 条评论