QIIME2导入数据-fq数据转换成qza数据-使用方法心得

启动QIIME2运行环境

conda activate qiime2-2019.4

###新建并定位设置到存在fq数据的文件夹

mkdir qiime2-importing-tutorial ##建立新的文件夹
cd */*/qiime2-importing-tutorial

在这里插入图片描述

将fq数据转换成qza数据

1.清单文件的建立

由于我们拿到的数据多为FASTQ数据即fq格式。如果你不是EMP或CASAVA格式的数据,则需要先创建一个“清单文件”,然后使用导入命令qiime工具导入命令,手动将数据导入到qiime 2

我们大都不是,所以必须建立清单来促进数据的导入。
清单文件格式如上清单文件是逗号分隔(即.csv)的文本文件。每行的第一个字段是qiime应该使用的样本名sample-id (注意:必须加-,不然报错,这是双端33格式的要求),第二个字段是绝对文件路径absolute-filepath(注意:必须加-,不然报错,这是双端33格式的要求),第三个字段是读取方向direction。以#开头的行和空行将被忽略。sample-id里是你自己的样品名,absolute-filepath里$PWD是目前绝对路径,后面加上你的样品名(双端的话,会有1,2两组数据,记住要在压缩状态,即*.fq.gz)。direction主要包括:forward和reverse,*_1.fq.gz为forward。

2.代码运行

在这之前,先讲一下,数据存放问题。fq数据在qiime2-importing-tutorial 下,同时新建文件夹SampleData,将清单文件放入。
在这里插入图片描述mkdir SampleData
上传清单文件
在这里插入图片描述运行代码

qiime tools import \
  --type 'SampleData[PairedEndSequencesWithQuality]' \
  --input-path SampleData/manifest.csv \
  --output-path  paired-end-demux.qza \
  --input-format PairedEndFastqManifestPhred33

在这里插入图片描述fq数据转换成qza数据,大功告成。

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值