sra-toolkit使用

参考网址:https://github.com/ncbi/sra-tools/wiki/08.-prefetch-and-fasterq-dump

(Sequence Read Archive, SRA)数据可通过多个云提供商和NCBI服务器获得,是最大的高通量测序数据公开可用存储库。该档案接受来自生命的所有分支以及宏基因组和环境调查的数据。SRA存储原始测序数据和比对信息,以提高可重复性,并通过数据分析促进新发现。

prefetch + fasterq-dump的组合是从sra -access中提取fastq -文件的最快方法。prefetch工具下载所有必要的文件到您的计算机。如果下载不成功,可以多次调用prefetch工具。它不会每次都从头开始;相反,它将从上次调用失败的地方开始。下载之后,您可以选择使用vdb-validate工具测试下载的数据。下载成功后,不需要网络连接。您可以将通过prefetch创建的文件夹移动到不同的位置,以在其他地方执行到fastq格式的转换(例如到没有互联网访问的计算集群)。

检查“prefetch”工具的最大大小限制

prefetch工具的默认最大下载大小为20G。如果请求的增量大于20G,则需要增加该限制。无论请求的访问有多大,您都可以指定一个极高的限制。你也可以使用vdb-dump-tool和--info选项查询access -size。例如,vdb-dump SRR000001--info告诉您这个登录有多大(以及其他信息)。接入SRR000001有932,308,473字节,低于默认限制,因此不需要进一步操作。接入SRR1951777有410,112,373,995字节。要下载此条目,您必须解除该大小以上的限制:

prefetch SRR1951777 --max-size 420000000000

从SRA 文件中提取fastq文件

在执行提取之前,您应该快速估计所需的硬盘空间。最终的fastq文件大小大约是sra文件的7倍。在转换过程中,fasterq-dump工具需要大约1.5倍于最终fastq文件量的临时空间(临时空间)。总的来说,您在转换过程中需要的空间大约是sra大小的17倍。如果超过限制,'fasterq-dump'-工具将失败,并将显示一条消息。

fasterq-dump SRR000001

运行fasterq-dump最简单的方法

运行该命令行的前提是,您之前已经将要下载的序列prefetched到当前工作目录。如果目录SRR000001不存在(即未进行prefetch获得文件),fasterq-dump工具将尝试通过访问网络下载该数据。这将慢得多,并可能最终由于网络超时而失败。

需要注意的事项

fasterq-dump tool默认执行split-3操作。在命令行选项方面,fasterq-dump-tool与以前的fastq-dump-tool并不相同。fastq-dump和fasterq-dump的对比如下:

split-3

fastq-dump SRR000001 --split-3 --skip-technical

fasterq-dump SRR000001

split-spot

fastq-dump SRR000001 --split-spot --skip-technical

fasterq-dump SRR000001 --split-spot

split-files

fastq-dump SRR000001 --split-files --skip-technical

fasterq-dump SRR000001 --split-files

concatenated

fastq-dump SRR000001

fasterq-dump SRR000001 --concatenate-reads --include-technical

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

推荐阅读更多精彩内容