烟台公司网站开发,如何开一家软件外包公司,做包装找灵感看什么网站,有没有免费的推广平台为什么介绍lotus2
因为快#xff0c;作者比较了lotus2流程和qiime2、dada2、vsearch等#xff0c;lotus2的速度最快、占用内存最小。 因为方便#xff0c;只需要一行代码#xff0c;即可完成全部分析。
lotus2 -i Example/ -m Example/miSeqMap.sm.txt -o myTestRun而且分…为什么介绍lotus2
因为快作者比较了lotus2流程和qiime2、dada2、vsearch等lotus2的速度最快、占用内存最小。 因为方便只需要一行代码即可完成全部分析。
lotus2 -i Example/ -m Example/miSeqMap.sm.txt -o myTestRun而且分类注释更加准确。
综合以上优点lotus2完全值得尝试接下来进入实战部分。
Lotus2的安装
推荐conda安装
conda install -c bioconda lotus2如果失败了可以尝试
conda create -c conda-forge -c bioconda --strict-channel-priority -n lotus2 lotus2安装成功后激活环境
conda activate lotus2细菌16s扩增子分析运行命令
lotus2 -i Example/ \ #序列所在位置的文件夹-m Example/miSeqMap.sm.txt \ #map文件示意文件参考后续内容-o myTestRun2 \ # 输出文件夹OTU table、tree、运行日志等都会放在这个文件夹中-s configs/sdm_miSeq2.txt \ #配置文件这个是作者提供的设定了运行命令所需的参数-p miSeq \ #测序平台比如我的真菌ITS测序用的是PacBio就填写PacBio即可-amplicon_type SSU \ #真菌ITS测序则填写ITS-forwardPrimer GTGYCAGCMGCCGCGGTAA \-reversePrimer GGACTACNVGGGTWTCTAAT \-CL dada2 \ #聚类的方法-refDB SLV \ #注释用的数据库-taxAligner lambda \ #注释的方法-t 14 #使用的线程数map文件的示意注意测序文件和map文件的所在目录。假如我们所在的工作目录为WD序列位置为WD/Example/seq.fastq.gz。
#SampleIDfastqFileSequencingRunANG1PID-0292-1_S1_L001_R1_001.fastq.sm.gz,PID-0292-1_S1_L001_R2_001.fastq.sm.gzCANG2PID-0292-1_S2_L001_R1_001.fastq.sm.gz,PID-0292-1_S2_L001_R2_001.fastq.sm.gzC
这里C的含义不是很清楚有了解的小伙伴可以留言一下。
真菌ITS扩增子分析运行命令
lotus2 -i YOURPATH \-p PacBio \-id 0.97 \-CL cdhit \-s configs/sdm_PacBio_ITS.txt \-refDB UNITE \-m ./map.txt \-o ./LotuS2 \-forwardPrimer GTACACACCGCCCGTCG \-reversePrimer CGCCTSCSCTTANTDATATGC \-t 14 \-amplicon_type ITS \-taxAligner vsearch \-buildPhylo 0 \ #一般认为ITS构建的遗传发育树不靠谱可以不用构建。#Recommendation from author: #We recommend the cautious usage of the #phylogenetic tree for ITS because high #variation of ITS sequences may lead to erroneous trees.
