Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Save golobor/fea8b74fce5e1b8119e6c8dc96671f07 to your computer and use it in GitHub Desktop.
Save golobor/fea8b74fce5e1b8119e6c8dc96671f07 to your computer and use it in GitHub Desktop.
For retrieving fastq file from NCBI GEO using GSM number.
#!/bin/sh
if ! type "efetch" > /dev/null; then
print "Please install E-utilitie."
fi
GSM_LIST=$1
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`
echo downloading $SRR
fastq-dump -A $SRR -O sra_fastq
echo fastq-dump done.....
done
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment