mkdir 20220923_CXQ
cd 20220923_CXQ
conda install -c hcc aspera-cli -y
ascp -h
which ascp
# /public/home/djs/miniconda3/bin/ascp
# 查找密钥
ls /public/home/djs/miniconda3/etc/asperaweb_id_dsa.openssh
#
cat download.txt |while read id ;do ascp -l 100M -P 33001 -QT -k 2 -i /public/home/djs/miniconda3/etc/asperaweb_id_dsa.openssh era-fasp@fasp.sra.ebi.ac.uk:/vol1/srr/${id:0:6}/00${id:9}/${id} ${id};done
更新
当我再次使用时发现上述口令不好使了,就是那个 URL构建的方式变了
可以直接下载 EBI提供的 filereport 文件,文件打开后格式如下,提取sra_aspera 那一列即可:
下载口令的构建如下:
cat filereport_read_run_PRJNA751490_tsv.txt | \
cut -f4 | \
while read id ; \
do ascp -l 100M -P 33001 -QT -k 2 -i /public/home/djs/miniconda3/etc/asperaweb_id_dsa.openssh \
era-fasp@${id} $(basename ${id}) ; \
done