Unable to retrieve the "Assembly method" of an assembly using Edirect in the terminal
1
0
Entering edit mode
5 months ago

Dear all , I am having a list of assembly ID and would like to retrieve specific info on the assembly method (eg. SPADes) from the docsum of that assembly. I have been trying these codes that worked for some other retrievable info, but I could not find the Assembly method although it is shown on the ncbi website. I also inspected the docsum itself where I can not find assembly method. I still see some tutorial using some commands as shown below but it did not work. Here is my code:

printing the docsum

esearch -db assembly -query "GCF_037000505.1" | efetch -format docsum 

I looked into this, there is no assembly method

Some tutorials suggested

efetch -db assembly -id GCF_037000505.1 -format docsum | xtract -pattern DocumentSummary -element AssemblyMethod

This did not produce any output.

Can you advise on a way in linux using esearch to retrieve the assembly method ? Why it is not retrieved via docsum although present on ncbi ?

Thanks

Assembly Sequencing Genomics • 615 views
ADD COMMENT
0
Entering edit mode

may be this is the reason:

https://ncbiinsights.ncbi.nlm.nih.gov/2023/10/18/ncbi-datasets-access-sequence-data/

Effective June 2024, NCBI Datasets will replace legacy Genome and Assembly web resources

ADD REPLY
0
Entering edit mode

I not sure, but when I used the same code to retrieve other info (e.g. AssemblyStatus), it came out correctly as shown in the docsum !

ADD REPLY
1
Entering edit mode

You can only extract things shown in docsum output which does not contain this info.

ADD REPLY
2
Entering edit mode
5 months ago
GenoMax 147k

Using datasets and dataformat tools from NCBI:

$ datasets summary genome accession GCF_037000505 --as-json-lines | dataformat tsv genome --fields organism-name,assminfo-assembly-method
Organism Name   Assembly Assembly Method
Klebsiella pneumoniae   SPAdes v. 3.7
ADD COMMENT
0
Entering edit mode

Thanks a lot that worked nicely.

ADD REPLY

Login before adding your answer.

Traffic: 1855 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6