我们进行获取数据并进行分析时，一个个的操作未免太过麻烦，因此可以利用shell脚本进行批量处理。在此利用之前老师上课讲过的生物软件来利用脚本做一个流程化处理。
首先我们明确一下操作流程
1、利用Aspera批量下载文件
2、利用FastQC语需要.fastq.gz文件，所以我们需要利用SRA-toolkits做格式转换
3、利用FastQC做质控
我的需要批量下载的序列放在sra_list.txt文件中
sra_list.txt
下面是一些结果和脚本代码
第一步结果
第一步结果
第二步结果
第二步结果
第三步结果
第三步结果
脚本运行结果
脚本运行结果
因为是在sra_list.txt中放了三个文件位置，却只下载了两条，有一条错误，我在下面也试了很多次，代码应该没有什么问题，估计可能是服务器不稳定，如果有其他什么发现，欢迎同学们告知
关于后面的序列组装SPAdes和velvet以及序列拼接评价QUAST的脚本，大家有兴趣的可以在后面添加，与此如出一辙。
今天看了俊辉的脚本发现可以对PE还是SE进行区别，大家也可以看看
我们进行获取数据并进行分析时，一个个的操作未免太过麻烦，因此可以利用shell脚本进行批量处理。在此利用之前老师上...
数据下载 数据转换 数据质控 数据过滤
下载单个文件脚本代码 上传单个文件脚本代码 批量下载脚本代码 批量上传脚本代码 命令解释 登录FTP << 是使用...
最近要下载fnl再分析数据，网站原有的压缩包下载功能没有了，只能通过脚本批量下载。在研究利用python，...
批量下载安装： 小工具：比对、cricos图、下载工具： 基因功能分析预测 基因组拼接 测序数据质控
上篇文章写了如何将protoco转成我们iOS 开发用的.pb.swiftSwift5 与 Protobuf由于p...
安装依赖 request 创建储存视频的文件夹 请求本地文件(txt) 下载视频
利用python自带工具2to3.py完成批量转换！
利用python脚本，批量打包碎图。会先更新svn，然后执行打包操作
创建项目 下载原始FastQ数据 1.1 使用循环创建下载地址文件 1.2 编写批量下载脚本aspera.sh 1...
本文标题：关于利用脚本批量下载、转换、质控
本文链接：https://www.haomeiwen.com/subject/pwgotqtx.html
网友评论