Dear colleague,
I am working on the analysis of prokaryotic genomes from NCBI genome database.
After opening the file, we could see one column named "Chromosome RefSeq". (e.g., NZ_AQXM00000000)
Each file has a name like "GCF_000834735.1_ASM83473v1_protein.faa.gz".
It is odd that the two datasets use different accession number system. In this case, how to identify if the genome of the proteins is annotated in the file "prok_representative_genomes.txt"？
My aim is to retrieve all the protein sequences for the genomes listed in the file "prok_representative_genomes.txt".
Thanks a lot,
Kind regards
Tom
In the same FTP path, there is another file prokaryotes.txt that can be helpful. Join the two files on the organism name (column 3 of prok_representative_genomes.txt and column 1 of prokaryotes.txt). Then, extract the path to the assembly folder on the FTP site (column 21 of prokaryotes.txt file) and download the *.protein.faa.gz files.
Accessions with the NZ_ prefix are RefSeq Chromosome accessions whereas the ones with GCF_ prefix are assembly accessions. An assembly includes both chromosome and plasmid(s).
As vkkodali wrote, use the prokayotes.txt. However, there's no need to join, one can filter on the Reference column. The NCBI discriminates between manually selected reference (exist for common lab organisms) and computationally selected representative genomes (exists for each species, almost).
Log in to answer this question.
Thanks,
Does anybody know why there are some 0 values in the "Genomes" column of the file prok_representative_genomes.txt? For example, "Acetobacter aceti".
The README says it's manually updated. I tend to believe you'll see a more complete picture when using the
prokaryotes.txt(computationally updated...)