This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Use esearch/efetch to output relationship table of GSM to SRR (bascially SRA file names)

Hello all,

I want to use NCBI command line utilities (esearch, efetch, etc) to achieve the following: given the GSE ID, I want a simple tab-separated output, giving you GSM experiment in column 1, and SRR file name in column 2. Ideally, it would also be collapsed by GSM ID (in cases when you have more than 1 sra file per GSM).

I'm reading through numerous pages of manuals that are relatively obscure for people without much experience with databases, so if you can help me to figure this out, I'd be most grateful :)

Thank you in advance!

geo geo-omnibus ncbi
• 5,387 views • link
written by predeus predeus

1 answer

Ok I've figured it out. To get GSM to SRR relationship, you can use this:

esearch -db sra -query "GSM123456" | efetch -format docsum | xtract -pattern DocumentSummary -element Runs | perl -ne '@mt = ($_ =~ /SRR\d+/g); print "@mt\n"'

note that it will also print multiple SRR IDs in one line (if you have more than one sra file per GSM).

• 0 views • link
updated by Ram Ram • written by predeus predeus

Can you give more information on

  1. if we should have a perl script to run the command?
  2. if it can be run in cygwin
  3. a url to the perl script if one exists.

Thank you.

• 0 views • link
updated by Ram Ram • written by an • an

Log in to answer this question.

More posts like this