如何使用BLAST将分类法分配给trnL序列

时间:2015-09-28 23:19:45

标签: sequence bioinformatics taxonomy blast qiime

我正在使用trnL叶绿体基因从食草动物粪中鉴定植物,目前我正试图将分类学分配给我的Illumina输出中的trnL序列。以下是我想要运行的QIIME脚本和选项:

assign_taxonomy.py -i rep_set_numbered.fa -r sequence.fasta -t id_to_taxonomy.txt -e 0.01 -m blast

我有来自数据管道的输入文件,以及来自NCBI GenBank的参考文件(205,703个序列)。但是,我没有制表符分隔的分类法文本文件。通常我会从Excel生成一个,但由于FASTA文件太大(超过500 MB),因此无法在Excel中完全查看,因此无法可靠地编辑。

我的问题是,是否有一个命令行方法用于从我的参考FASTA文件生成我自己的制表符分隔的分类文件,如果有,我该怎么做?如果没有,在" assign_taxonomy.py"处理这个必需选项的其他选择是什么? QIIME脚本?

0 个答案:

没有答案