I am trying to use my list of protein IDs (e.,g VUX63899.1,QDO61010.1, QDO50771.1) to retrieve the assembly accession number.
Using the follow commands ..
epost -input teste.txt -db protein | elink -target nuccore -db protein | elink -target assembly |esummary | xtract -pattern AssemblyAccession -element AssemblyAccession > assemblyAccession.txt
However, I need two columns with assembly accession with its corresponding protein IF. I could not figure out yet what the command I need to add.
Thanks in advance.