This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to change GCA_ to NC_ annotations (RefSeq)

I downloaded the data_summary.tsv file via NCBI Datasets, but I need to change the GCA_ style annotations in the file to NC_ annotations in order to create a RefSeq (NC_) reference list. Is there any way to convert RefSeq Assembly to RefSeq accession number?

ncbi refseq
• 1,710 views • link
updated by Istvan Albert Istvan Albert • written by beginner123 • beginner123

it would be helpful if you can post what you have downloaded and what you want.

• 0 views • link
written by cpad0112 cpad0112

GCA accessions are GenBank assemblies where as corresponding GCF accessions (if they exist) would be RefSeq.

One way to convert these would be using EntrezDirect:

$ esearch -db assembly -query GCF_000266945 | elink -target nuccore | efetch -format acc
NC_018026.1
NC_018025.1
CP003361.1
CP003360.1
• 0 views • link
written by GenoMax GenoMax

looks like we wrote the exact same thing at the same time :-) - just took a detour in investigating the assembly stats

• 0 views • link
written by Istvan Albert Istvan Albert

1 answer

Interesting question I assumed the assembly summary would have that - turns out it does not. Upon some trial and error it seems you can link this up in the following way using Entrez Direct:

esearch -db assembly -query GCA_009858895 | elink -target nuccore | efetch -format acc

prints:

NC_045512.2
MN908947.3

it shows both the RefSeq and the GenBank entries of the same data.

• 0 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this