GEO的各种英文缩写
GPL GEO Platform 平台
GSE GEO Series
SRA-----Sequence Read Archive
存档来自高通量测序平台的原始测序数据和比对信息，包括Roche 454GSSystem®，Illumina GenomeAnalyzer®，Applied BiosystemsSOLiDSystem®，HelicosHeliscope®，CompleteGenomics®和PacificBiosciencesSMRT®。
SRA - SRA submission accession
SRP - SRA study accession
SRX - SRA experiments accession
SRR - SRA run accession
SRS - SRA sample accession
SRZ - SRA analysis accession
下机数据分析常规流程
image.png
安装aspera
####aspera支持NCBI断点续传，同样可以提高下载速度？
wget -c https://download.asperasoft.com/download/sw/connect/3.8.1/ibm-aspera-connect-3.8.1.161274-linux-g2.12-64.tar.gz
tar zxvf ibm-aspera-connect-3.8.1.161274-linux-g2.12-64.tar.gz
cd /home/xjqi/.aspera/connect/bin
./ascp --help
数据下载
可以下载sra文件或fastq文件
sra文件为fastq.z又经压缩后的内容
用prefetch命令下载文件会直接调用ascp（添加到环境变量后，这里我还没搞通透）
我们需要的信息中的变量其实是GSEXXXXXX，根据链接格式进行更改
EBI数据库下载fastq.gz 列表----
fastq.gz文件. https://www.ebi.ac.uk/ena/data/search?query=GSE111229
image.png
image.png
image.png
GEO下载sra文件列表----
sra文件. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111229
image.png
image.png
image.png
####PRJNA436229.txt为fastq.gz的列表，prefetch.txt为SRR的list，是我们需要下载内容的名称列表
cat PRJNA436229.txt |grep 'gz'|cut -f 10|head -n 2|while read line;do wget $line;done
cat prefetch.txt |while read id;do prefetch $id -O ~/;done
####可以看到下载下来的文件有很多NC、KC、DQ开头的，看了一下是二进制文件，但具体什么情况，尚不了解。
image.png
SRR转为fastq数据
####Single-end 43 bp long reads were aligned at ESCG to mm10 mouse genome.这是数据搜索时看到的测序策略描述。
####single end和mate pair（pair end）数据转换时参数是不同的？
fastq-dump --gzip -O ./ ./SRR1039508.sra
