单细胞原始SRA数据下载

拟下载10X基因组数据GSE134520,该数据集对应的SRA编号为SRP215370

这个数据集一共测了9个病人的单细胞数据(包括了3个浅表性胃炎、3个慢性萎缩性胃炎、6个肠上皮化生、1个早期胃癌)PS:由于数据太大,这里就不全部下载,只选择其中几个样本进行下载分析.

** 0、配置环境

conda create -n sra-toolkit
conda activate sra-toolkit

一、下载SRA数据前的准备

1、首先要先安装prefetch 和ascp

Prefetch是sratools中的一个小工具,可以直接安装:

conda install -c daler sratoolkit

可以显示帮助文档就说明安装成功

prefetch –h

如果要下载数据比如SRR文件,直接加ID号,指定输出目录就好

prefetch SRRxxxxxxx -O PATH

2、默认情况下,prefetch是利用https方式去下载原始数据,这个就像直接从网页下载一样,速度有一定的限制。因此这里需要先安装一款叫做"aspera"的下载工具,它是IBM旗下的商业高速文件传输软件,与NCBI和EBI有协作合同

** 二、下载aspera安装包

1.aspera下载
先进入https://www.ibm.com/products/aspera/downloads?list
然后找到如图所示的IBM Aspera Connect。因为数据一般是下载在服务器上,所以需要选择linux版本进行安装。

image.png

wget http://download.asperasoft.com/download/sw/connect/3.7.4/aspera-connect3.7.4.147727-linux-64.tar.gz
image.gif

将下载的安装包解压

tar zxvfibm-aspera-connect-3.9.6.173386-linux-g2.12-64.tar.gz

安装aspera

sh ibm-aspera-connect-3.9.6.173386-linux-g2.12-64.sh 

查看是否安装成功

cd /home/hej/.aspera/connect/bin

显示帮助文档就说明安装成功

./ascp –help

永久添加环境变量

echo 'export PATH=~/.aspera/connect/bin:$PATH'>> ~/.bashrc

最好写绝对路径

echo 'export PATH=/home/limf/.aspera/connect/bin:$PATH'>> ~/.bashrc

最后检查ascp是不是能用了

ascp –help

*ascp安装成功后,prefetch就会默认将下载方式从https转移到fasp,说明开启加速模式

二、开始下载SRA数据

1、将要下载的数据对应的SRR编号存在文件SRR_Acc_List.txt中

那么如何获取想要下载的数据集样本对应的SRR编号呢?

方法如下:

① 进入NCBI-SRA界面,选择SRA Run Selector

图片

② 输入SRP215370

图片

③下载Accession List

图片

Accession List具体内容如下:

图片

2、运用prefetch和aspera下载单个文件

对单个文件进行下载

prefetch SRR9713119-O `pwd` && echo "** SRR****.sra done **"

下载完成时会显示如下

图片

3、运用prefetch和aspera进行批量下载

对批量文件进行下载,将依次下载SRR_Acc_List.txt文件中SRR编号对应的数据

cat SRR_Acc_List.txt| while read i

4、有时候SRR_Acc_List.txt文件中列出来的数据,通过prefetch和aspera有可能下载不下来,所以这里介绍另一种补充的下载方式。

①、进入EBI官网https://www.ebi.ac.uk/ena

搜索SRA编号

图片

②、运用EBI下载的好处是可以选择直接下载fastq格式的文件,也可以选择下载SRA格式的文件。可以根据以下步骤获得fastq文件或者SRA文件所在的链接。

图片

③ 、可以直接在网页下载也可以将获得的链接粘贴到以下代码中进行下载

ascp -QT -l 300m-P33001 -i ~/.aspera/connect/etc/asperaweb_id_dsa.openssh era-fasp@fasp.sra.ebi.ac.uk:vol1/srr/SRR971/009/SRR9713119./

转载来自:
李呆呆的学习笔记

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,826评论 6 506
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,968评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,234评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,562评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,611评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,482评论 1 302
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,271评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,166评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,608评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,814评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,926评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,644评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,249评论 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,866评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,991评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,063评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,871评论 2 354

推荐阅读更多精彩内容