生物信息分析人员一般会接触到从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. 皮一皮:男生的求生欲,以及神回复...
  2. H3C PPP MP配置示例二(续)
  3. github建立静态网站,域名解析和跳转
  4. oracle的服务关闭,oracle 各种服务关闭启动
  5. [HDU 6643] Ridiculous Netizens(点分治+根号分治+dp)
  6. 2012年3月51CTO壁纸点评活动获奖名单【已结束】
  7. Leetcode每日一题:147.insertion-sort-list(对链表进行插入排序)
  8. 安装deepin后会删除win10吗_不管是先装Win10还是先装Deepin,启动问题轻松搞定
  9. 人民日报发推欢迎Google重返大陆,FB上长文阐述详细立场
  10. spring 4.x下让http请求返回json串
  11. Julia : 用@async提升循环性能
  12. 2022第四届中国制造业数字化转型国际峰会
  13. 43岁,转行当了大学老师
  14. 随手查_python
  15. python 数据库框架peewee_python peewee模块使用
  16. 【优秀易用的无损音乐播放器】Amarra Luxe for Mac 4.3
  17. ffmpeg 音频解码二
  18. js简单实现一个json格式化功能
  19. java 与 c gzip_zlib和gzip生成不同的数据
  20. VScode受难记 - 0

热门文章

  1. w8ndows 秒表,Win8开机时间:机械硬盘堪比SSD_固态硬盘评测-中关村在线
  2. leetcode系列-70. 爬楼梯
  3. 商品审核网页界面_网页设计与制作竞赛通知
  4. 从一道笔试题谈算法优化(上)
  5. 【工具】audacity提示声音设备出错
  6. DTMF双音频信令检测,基于Goertzel(戈泽尔算法),FPGA实现
  7. libdhnetsdk.so: 错误 ELF 类: ELFCLASS64
  8. 为什么深度学习需要与边缘计算进行结合?如何结合?
  9. GIS—DEM与数字地形分析
  10. win2003 iis运行php环境,win2003 IIS6.0配置php环境如何加载php5isapi.dll