生物信息分析人员一般会接触到从NCBI等网站下载的SRA数据,之前也介绍了下载SRA数据的几种方式。下面,我就简单介绍一下如何将下载的sra格式数据转换成为常用的fastq等格式。

1、fastq-dump命令

sratoolkit的下载,该部分详见上一篇文章(SRA数据库及linux本地下载)

单端测序:

fastq-dump SRR14306907.sra -O ./ (结果生成:SRR14306907.fastq)
fastq-dump --fasta SRR14306907.sra -O ./ (结果生成:SRR14306907.fasta)

双端测序:

fastq-dump SRR14306907.sra --split-3 -O ./ (结果生成:SRR14306907_1.fastq,SRR14306907_2.fastq)
fastq-dump SRR14306907.sra --split-3 --gzip -O ./ (结果生成:SRR14306907_1.fastq.gz, SRR14306907_2.fastq.gz)

2、pfastq-dump

1)下载pfastq-dump

git clone https://github.com/inutano/pfastq-dump
cd pfastq-dump/bin/
chmod a+x pfastq-dump

复制到自己的文件夹

cp pfastq-dump /home/xuyang/SRAToolkit/sratoolkit.2.10.9-ubuntu64/bin

2)转换格式

单端测序:

pfastq-dump SRR14306907.sra -O . -t 10

循环

for id in *sra; do pfastq-dump --threads 10 ./$id --gzip; done

双端测序:

 pfastq-dump SRR14306907.sra --split-3 --gzip -O ./ -t 10

循环

for id in *sra; do pfastq-dump --threads 8 ./$id --split-3 --gzip; done

for i in `tail -n+1 sra_new.ids|cut -f1`;do
pfastq-dump ${i}.sra --split-3 --gzip -O ./ -t 10
done 

其中sra_new.ids为所有srafile的ID号,t为线程数

3.批量转换

for循环函数

for i in `tail -n+1 sra_new.ids|cut -f1`;do
pfastq-dump ${i}.sra --split-3 --gzip -O ./ -t 10
done 

其中sra_new.ids是sra文件的列表

ref:inutano/pfastq-dump: parallel-fastq-dump implementation in bash script (github.com)

NGS小技能(2):如何进行SRA到fastq格式的快速转换 - 简书 (jianshu.com)

SRA到fastq格式的批量转换相关推荐

  1. 视频文件的格式可以批量转换吗

    在我们的工作生活中,如果遇到视频格式不适用而导致视频无法正常播放或者使用情况,有什么办法可以快速解决吗?今天小编就来分享一下如何批量转换多个视频格式,快来看看具体的操作方法吧. 首先我们需要先下载安装 ...

  2. 前端 Git CRLF 与 LF格式的批量转换

    问题出现的原因: 从远程下代码,发现有很多的变更,但文件内容没变,最后发现是换行符不同导致的,遂需要批量解决. 实现方法: 1.点击下载 UNIX2DOS工具.下载后是一个ZIP的压缩包,解压到某目录 ...

  3. Mac下.CR2格式图片批量转换

    我的相机是佳能X7i型号,之前拍的照片都是.cr2格式,想上传照片到Flickr,因为格式原因上传不了,开始一张张手动格式转换,太傻瓜.直到今天,发现了这个,感觉太好了. No need for sl ...

  4. excel排版成印刷版格式:批量转换成pdf

    下载地址:http://blog.csdn.net/bq_cui/article/details/49685983 插入页码后的excel,只需要1键即可生成pdf文件. 点击"全部转PDF ...

  5. html5 mp4转换ogv格式,FFmpeg 批量转换HTML5视频(mp4/webm/ogv)

    现在Html5也逐渐普及了,设备对于mp4/webm/ogv的支持也逐渐广泛,但是Html5视频的转换软件却不是太多,这个Ffmpeg确实不错.而且功能也强大. REM "H.264 / A ...

  6. 怎么批量转换图片格式?

    怎么批量转换图片格式?在日常的工作过程中,小伙们时常需要将上百上千的图片进行格式转换,如果你一张一张图片进行转换,就需要消耗大量的时间和精力,时间就是金钱,我们不应该浪费大量时间在这种简单又重复的事情 ...

  7. c#如何wmf图片转换成png图片_C# 批量转换图片格式 支持/JPG/TIF/PNG/WMF/GIF等

    C# 批量转换图片格式实例源码,可以在BMP/JPG/TIF/PNG/WMF/GIF这几种格式之间批量转换,是否支持逆转换,这个未测试,有兴趣的自己下载源码编译一下吧.本源码实现的原理:将转换过程写入 ...

  8. bmp格式如何转换成jpg?可批量转换

    bmp对于很多小伙伴们来说可能比较陌生,因为接触使用的比较少,但它也是Windows操作系统中的标准图像文件格式,小编在百度里经常下载到,由于这种图片格式不进行压缩,所以比较占用内存空间,并且很多网站 ...

  9. RNA-seq流程学习笔记(4)-使用FastQC软件对fastq格式的数据进行质量控制

    今天开始学习使用FastQC软件对范例SRA测序文件的质量进行分析. 主要参考文章: RNA-seq(3):sra到fastq格式转换并进行质量控制 转录组入门(3):了解fastq测序数据 用Fas ...

最新文章

  1. 《css世界》的那些实用技巧
  2. angular中集中页面传参(我只是知识的搬运工)
  3. 关于水晶报表的一些错误
  4. I/O多路转接之poll——基于TCP协议
  5. matlab三角形分割,MATLAB 2014b及以上版本中带有画家渲染器的三角形拆分补丁
  6. linux服务器下降,linux - 远程升级Ubuntu:如何最大程度地降低丢失服务器的风险? - Ubuntu问答...
  7. 数据仓库之电商数仓-- 2、业务数据采集平台
  8. 3层b+树索引访问磁盘次数_从B+树到LSM树,及LSM树在HBase中的应用
  9. UITablView索引列表
  10. 学习云计算有什么用?企业怎么才能“上云”?
  11. xutils使用流程
  12. java ee是编程语言_什么是Java编程语言?
  13. docker 重启gitlab_gitlab从入门到绝望
  14. android 微信小程序原理,Android开发微信小程序页面的图文教程
  15. 微信公众号群发图文API插入视频音频解决方案
  16. 雅思备考-口语词组积累-第一章
  17. hse是什么意思的缩写_什么是HSE
  18. W2019SRV and W10 (1809) Office2019 CN
  19. opencv图像处理学习(六十八)——肤色检测
  20. 基于C++11的muduo网络库

热门文章

  1. 微信消灭病毒哪个服务器好,消灭病毒120和140哪个好用
  2. 求两个数最大公因数问题,辗转相除法
  3. 51单片机静态动态数码管显示
  4. 计算机维修工作案例分析,西门子工控设备维修案例分析
  5. 6G中的卫星通信高效天基计算技术
  6. 基于Python的旅游管理系统微信小程序设计与实现毕业论文+项目源码及数据库
  7. 世界首款骨传导蓝牙立体声耳机问世
  8. 1399元起!荣耀9X发布:全系搭载7nm麒麟810,升降全面屏,主流手游满帧体验
  9. 小米路由器R2D刷机
  10. [图文教程]阿里云万网域名一键签发SSL证书 快捷申请免费SSL证书