Fastq文件的获取

news2024/9/22 23:24:42

Fastq文件的获取

  • Fastq文件的获取
      • linux 下安装SRA Toolkit工具
      • 使用SRA Toolkit工具下载SRA数据

Fastq文件的获取

author:CYH-BI

date:2023.10.25

Fastq文件得获取 ,可以从NCBI获取,也可以自己测序得到结果,下面介绍从网上获取数据。

linux 下安装SRA Toolkit工具

2023.3.24

1、NCBI 官网查看工具地址,并下载安装包

具体操作步骤如下:

  • 打开NCBI 官网

    在这里插入图片描述

  • 点击Download Tools并进入

    在这里插入图片描述

  • 找到SRA Toolskit

在这里插入图片描述

  • 会转跳到github ,复制下载连接,后续使用Linux wget 命令下载

    在这里插入图片描述

  • 打开Linux 总端,使用wget 命令下载

    wget https://ftp-trace.ncbi.nlm.nih.gov/sra/sdk/3.0.7/sratoolkit.3.0.7-centos_linux64.tar.gz
    

2、解压sratoolkit.3.0.7 安装包.。注意包名称,我们的不一定一样(不同版本不一样)

tar xvf sratoolkit.3.0.7-centos_linux64.tar.gz

3、查看有无安装成功(如果软件含有bin目录,大多情况下配置环境时路径需要精确到bin)

cd /home/cyh/sratoolkit/sratoolkit.3.0.7-centos_linux64/bin #看看能不能打开bin

4、配置SRA Toolkit工具环境

nano ~/.bsahrc  #编辑.bashrc文件 或者是 vim /etc/.bashrc

export PATH="/home/cyh/sratoolkit/sratoolkit.3.0.7-linux64/bin:$PATH"  #最后一行输入
#编辑结束后按CTRL+C保存,按提示后按回车
# 格式:export PATH="<路径>:$PATH" ;路径要精确到bin
#注意:你的全路径。是你的,不是我的。
#注意:你的全路径。是你的,不是我的。
#注意:你的全路径。是你的,不是我的。
source ~/.bashrc  #重启文件

SRA Toolkit 工具的目的是下载NCBI未处理的SRA数据,还能将SRA文件处理成fastq文件

使用SRA Toolkit工具下载SRA数据

1、使用命令下载感兴趣的数据,SRA Toolkit 工具有个命令可以稳定下载数据,就是有点慢,但是稳定。

[root@localhost sra_data]# prefetch SRR12415658  ##以SRR12415658为例

以SRR12415658数据为例:

建议使用后台下载:(如果是服务器,后台下载不会因为连接中断而中断,除非服务器网断了)

nohup prefetch SRR12415658 &   

对于多个数据下载:(方法多种多样,可以自己写循环)

2、下载后,使用fastq-dump拆分sra数据,转变为fastq格式数据

fastq-dump --split-3 SRR12415658.sra

--split-3表示双端测序

加以加上 -o 参数,不加就会下到sra 里面(好像在SRA Toolkit里,忘记了,上次找了好久)。

fastq-dump --split-3 SRR12415658.sra -o /home/cyh/Desktop/sra

将结果输出到 /home/cyh/Desktop/fastq ,(文件你得先有)

3、还需了解其他参数,可以使用 -h参数。比如,如果你需要将文件输出到其他地方,加个 -o参数。

fastq-dump -h

后续可以得到两个fastq文件 ,用于后续FastQC 质控。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1131984.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

MySQL3:MySQL中一条更新SQL是如何执行的?

MySQL3&#xff1a;MySQL中一条更新SQL是如何执行的&#xff1f; MySQL中一条更新SQL是如何执行的&#xff1f;1.Buffer Pool缓冲池2.Redo logredo log作用Redo log文件位置redo log为什么是2个&#xff1f; 3.Undo log4.更新过程5.InnoDB官网架构InnoDB架构-内存结构①Buffer …

【uniapp+云函数调用】人脸识别,实人认证,适用于app,具体思路解析,已实现

2023.10.8 需求: uniapp开发的app项目中使用人脸识别 app项目都是第一次搞,更别提人脸识别了。目前已有的就是Dcloud账号已申请,实现需求的时间没那么紧迫 此篇会详细记录从0到1的过程 2023.10.24 今天开始探究实现的过程 可能会记录的有些冗余 效果图如下: uniapp开发指南…

GoLong的学习之路(九)语法之结构体(非常重要,不看就等于不会Go语言)

书接上回&#xff0c;上回书说到&#xff0c;map等数据结构&#xff0c;接下来说结构体 文章目录 类型别名和自定义类型自定义类型别名类型定义和类型别名的区别 结构体结构体的定义结构体实例化基本实例化 匿名结构体创建指针类型结构体并取用获取结构体指针取用结构体 结构体…

NewStarCTF2023week4-More Fast(GC回收)

打开链接&#xff0c;存在很多个类&#xff0c;很明显是php反序列化漏洞利用&#xff0c;需要构造pop链 &#xff0c; 关于pop链构造的详细步骤教学&#xff0c;请参考我之前的博客&#xff0c;真的讲得很详细也容易理解&#xff1a; http://t.csdnimg.cn/wMYNB 如果你是刚接…

降本增效神器?Share Creators 智能数字资产管理系统真香!

降本增效似乎是一个持续又永久的话题。尤其在今年&#xff0c; 显得格外的重要~ 在各行各业都受到了疫情所带来巨大冲击的背景下&#xff0c;降本增效对很多企业来说不再是锦上添花&#xff0c;而可能是一条唯一的出路。 随着市场的收缩和竞争的加剧&#xff0c;在更“卷”的…

树形数据增删改查

功能描述&#xff1a; 默认展示所有项目点击项目展示当前项目下的所有区域点击区域展示当前区域下的所有工位以上以树形图格式展示项目&#xff0c;区域&#xff0c;和工位都可进行增加 修改 和删除&#xff0c;每个图标hover时显示对应提示信息项目&#xff0c;区域&#xff…

数据库数据恢复—Oracle数据库报错ORA-01110错误的数据恢复案例

Oracle数据库故障&#xff1a; 北京某公司一台运行oracle数据库的服务器&#xff0c;机房意外断电导致该服务器重启&#xff0c;重启后发现oracle数据库报错。该Oracle数据库没有备份。 Oracle数据库数据恢复过程&#xff1a; 1、北亚企安数据恢复工程师检查该oracle数据库的数…

外卖跑腿小程序开发是否需要定期更新和维护?

外卖跑腿小程序已成为现代生活的一部分&#xff0c;但它们的成功和可靠性取决于持续的定期更新和维护。本文将探讨为什么外卖跑腿小程序需要定期更新和维护&#xff0c;以及如何实施这些关键的技术措施。 为什么需要定期更新和维护&#xff1f; 1. 改进性能和稳定性 外卖跑…

Android Studio模拟器/虚拟设备连接互联网的方法

如图&#xff0c;无线、网络都无法联网 找到本机的DNS 找到emu-launch-params.txt&#xff0c;添加DNS -dns-server 192.168.124.1 重启虚拟机&#xff0c;关闭无线

Python使用psycopg2读取PostgreSQL的geometry字段出现二进制乱码

1、问题 读取geometry字段出现二进制乱码 查询语句&#xff1a; sql "select * from public"Note: 这种写法在PostgreSQL中直接查询, 没有问题&#xff0c;不会报错。 但是在Python中查询&#xff0c;如果导出的geom还是一长串的geometry 格式的话&#xff0c; …

Anaconda创建环境出现一堆DEBUG

问题描述如图所示&#xff1a; 试了好几种方法&#xff0c;比如以为是清华的源不能用了嘛&#xff0c;然后改了之后发现不行&#xff1b;然后又以为是anaconda出啥问题了&#xff0c;卸了重装发现都没能解决。 但是好在没有放弃啊哈哈&#xff1a;【试试下面这个】 解决办法&a…

clion安装C++远程linux开发并调试 从装centos虚拟机到完美开发调试

下载镜像 阿里云镜像 从vmware上安装虚拟机并提权开放ssh 更新编译环境 一般 gcc gdb 版本都是比较低的&#xff0c;适配不了clion的最低要求。 升级gdb参考博客 升级gcc参考博客 安装CMAKE 官方的源使用wget即可下载 未找到openssl的解决办法 注意版本&#xff0c;又踩个…

C语言--有 n 个整数,使前面各数顺序向后移 m 个位置,最后 m 个数变成最前面的 m个数,见下图,写一个函数实现该功能

图解&#xff1a; 代码&#xff1a; void Move(int* arr, int n, int m) {if (m < 0 || m > n)return;//创建m个长度的int数组int* brr (int*)malloc(m * sizeof(int));assert(brr ! NULL);if (brr NULL)return;//把后面的m个数字移到新数组int i;for (i 0; i <…

【Axure高保真原型】图片手电筒效果

今天和大家分享图片手电筒效果的原型模板&#xff0c;鼠标移入图片区域后&#xff0c;会显示一个光圈&#xff0c;光圈会跟随鼠标移动&#xff0c;照亮对应的区域&#xff1b;鼠标拖动时可以移动地图图片&#xff0c;查看更多区域的内容&#xff0c;具体效果可以打开下方原型地…

【linux】安装rpmrebuild

rpmrebuild是一种从已经安装的包中构建RPM文件的工具。它可以用于轻松构建修改后的包&#xff0c;并适用于任何使用RPM的Linux发行版。 访问地址 rpm rebuild download | SourceForge.net 选择版本 版本地址&#xff1a;版本地址 下载安装包 安装 rpm -ivh rpmrebuild-2.15…

探索光模块的MSA多源协议

在当今高度互联的世界中&#xff0c;光模块作为网络设备的重要部分&#xff0c;其性能和质量直接影响到整个网络系统的运行。其中光模块由于其灵活性和高效性&#xff0c;已经成为数据中心和云计算领域的主要选择。本文易天光通信将深入探讨光模块的MSA协议&#xff0c;揭示其重…

6. Cesium中的Entity

1. Entity类简介 Entity类是Cesium中描述和呈现地球上实体对象的核心类。它具有丰富的属性和方法&#xff0c;用于控制和定制地理实体的外观和行为。Entity对象可以表示各种地理实体&#xff0c;如点、线、面等&#xff0c;并具有位置、方向、模型、标牌、折线、多边形等属性&…

小主机折腾记18

这个月冲动消费了小两千块钱…… 1.880g5twr 由于四根2400t的内存条没有用出去&#xff0c;我又把它们装回了惠普的800g5twr&#xff1b; 看到pdd有400块钱的9350K&#xff0c;于是想着给他上一个9350k 在参考了pdd、咸鱼以及淘宝的价格后&#xff0c;我发现400块钱的9350k都…

易点易动设备管理系统:提升企业备件管理和维修效率的智能解决方案

在当今竞争激烈的商业环境中&#xff0c;大型生产制造企业面临着日益复杂的设备管理挑战。为了提高生产效率和降低成本&#xff0c;企业需要一种高效的设备管理系统。易点易动设备管理系统是一款智能化的解决方案&#xff0c;旨在帮助企业提升备件管理和维修效率。本文将详细介…

【外贸干货】广交会后9大细节推进客户下单

10月的展会正如火如荼的进行&#xff0c;外贸人也到了真正忙碌的时刻。我们在展会中挖掘到的客户是否已经得到妥善跟进&#xff1f;如何进行跟进&#xff0c;以吸引客户的注意&#xff1f;如何撰写高回复率的开发信&#xff1f;不同的展后跟进方式将直接影响我们的参展成果&…