sh 从sra #fastq获取东西

Posted 2021-05-07

esearch -db sra -query PRJNA40075  | efetch --format runinfo | cut -d ',' -f 1 | grep SRR > seqs.tsv

while read p; do
  fastq-dump --split-files --gzip $p
done <../seqs.tsv