Phylogenetic analysis and docking study of neuraminidase gene of influenza A/H1N1 viruses circulating in Iran from 2010 to 2019

•In 2015–2016 flu season, A/H1N1/pdm09 was predominant all around the world and in our region.•Software including MEGA-X, MODELLER, UCSF ChimeraX, auto-dock 4.2, and other online tools were used to analyze the phylogenetic relationship, antiviral resistance and vaccine efficiency of selected NA proteins.•Iranian isolates were located in clade 6B similar to A/Michigan/45/2015 reference sequence.•The phylogenetic study showed that most Iranian NA sequences (between 2015 and 2016) were located in a single clade and following years were located in its subclade by 3 major mutations (G77R/K, V81A, and J188T).•Resistant mutations in drug targets of NA including I117M, D151E, I223V, and S247N were ascertained in 10 isolates during the 2015–2016 flu seasons.