从ENA中下载fastq数据

#!/bin/bash
#注意，输入的参数全部为绝对路径！！！

SRR=$1
sra_dir=/media/whq/282A932A2A92F3D2/WHQ/mysradata/sra/srafile/sra
genome=/media/whq/282A932A2A92F3D2/WHQ/reference_genome/Homo_sapiens.GRCh38.dna.toplevel.fa
gtf=/media/whq/282A932A2A92F3D2/WHQ/reference_index/ensembl_GRCH38_hisat2_index/Homo_sapiens.GRCh38.99.gtf
index=/media/whq/282A932A2A92F3D2/WHQ/reference_index/ensembl_GRCH38_hisat2_index/ref.fa
get_count_trans=/media/whq/282A932A2A92F3D2/WHQ/myRscript/get_count_trans.R
get_fastq_url_from_ena=/media/whq/282A932A2A92F3D2/WHQ/myRscript/get_fastq_url_from_ena.R
#获取ftp链接
cd $sra_dir
Rscript $get_fastq_url_from_ena $SRR
#从ENA中下载fastq序列
rm ascp_com2
cat url.list.txt|while read id;do id=${id#*vol1/} ; echo "ascp -QT -l 300m -P33001 -i ~/.aspera/connect/etc/asperaweb_id_dsa.openssh era-fasp@fasp.sra.ebi.ac.uk:/vol1/${id} ./" ; done>ascp_com2
echo "已经建立好下载通道"
./ascp_com2