Skip to content

Instantly share code, notes, and snippets.

Last active April 2, 2020 12:15
  • Star 1 You must be signed in to star a gist
  • Fork 1 You must be signed in to fork a gist
Star You must be signed in to star a gist
What would you like to do?
For retrieving fastq file from NCBI GEO using GSM number.
if ! type "efetch" > /dev/null; then
print "Please install E-utilitie."
GSMs=`cat $GSM_LIST|cut -f1`
#echo $SRRs
! type "foo" > /dev/null 2>&1;
for GSM in $GSMs; do
echo $GSM retrieves from NCBI GEO.....
SRR=`esearch -db sra -query $GSM |efetch -format docsum |xtract -pattern DocumentSummary -element Run@acc`
fastq-dump -A $SRR -O sra_fastq
echo fastq-dump done.....
Copy link

abmmki commented Apr 2, 2020

not working

/ GSM.List
GSM4432378 retrieves from NCBI GEO.....
unrecognized option: '-f'
fastq-dump done.....

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment