Question: Extracting the assemblies from wgs accessions
0
gravatar for genomes_and_MGEs
8 months ago by
genomes_and_MGEs0 wrote:

Hey everyone,

I just blasted a sequence against the wgs DB on NCBI. When I blast against the nt DB, I can import the accessions on Batch Entrez and directly fetch the assemblies. However, I can't do that when uploading the hits against the wgs DB. For example, I have these accession hits:

PJQP01000010.1 PJCO01000020.1 PJCN01000016.1 PESJ01000016.1 LKKS01000134.1 JFBC01000520.1 NEJN01000012.1 LQXW01000009.1 JWLN01000008.1 RJAI01000031.1 NTME01000046.1 MJVE01002323.1 LDPF01000101.1 JOJY01000001.1 AXUP01000411.1 ...

How can I extract the assemblies for all these accessions? Thank you guys! Cheers

sequence assembly genome • 158 views
ADD COMMENTlink modified 8 months ago by vkkodali1.3k • written 8 months ago by genomes_and_MGEs0
1
gravatar for vkkodali
8 months ago by
vkkodali1.3k
United States
vkkodali1.3k wrote:

Batch Entrez is expecting a list of GIs. If you can get those somehow, you can use Batch Entrez. Alternately, assuming you want to download the sequences in FASTA format, you can use Entrez Direct for this as follows:

epost -input <accs_file> -db nuccore -format acc \
  | efetch -db nuccore -format fasta

where <accs_file> is a text file with accessions, one accession per line. You have options other than -format fasta; check them out in the documentation.

ADD COMMENTlink written 8 months ago by vkkodali1.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1775 users visited in the last hour