| 查看: 1966 | 回复: 1 | ||
[求助]
QIIME做OTU分析去除嵌合体出错
|
|
我用Virtualbox 5.0.12，装的QIIME-1.9.1-amd64，做OTU分析。 在做第一步去嵌合体时，用了以下命令：identify_chimeric_seqs.py -i trim.fasta -m usearch61 -o check_chimeras -r /home/qiime/diversity/analysis/silva.gold.ng.fasta 注：我已经将silva.gold.ng.fasta放在了analysis文件夹了。 结果运行后出现下列错误： Traceback (most recent call last): File "/usr/local/bin/identify_chimeric_seqs.py", line 354, in <module> main() File "/usr/local/bin/identify_chimeric_seqs.py", line 350, in main threads=threads) File "/usr/local/lib/python2.7/dist-packages/qiime/identify_chimeric_seqs.py", line 774, in usearch61_chimera_check log_lines, verbose, threads) File "/usr/local/lib/python2.7/dist-packages/qiime/identify_chimeric_seqs.py", line 891, in identify_chimeras_usearch61 consout_filepath=output_consensus_fp) File "/usr/local/lib/python2.7/dist-packages/bfillings/usearch.py", line 2280, in usearch61_smallmem_cluster app_result = app() File "/usr/local/lib/python2.7/dist-packages/burrito/util.py", line 285, in __call__ 'StdErr:\n%s\n' % open(errfile).read()) burrito.util.ApplicationError: Unacceptable application exit status: 139 Command: cd "/home/qiime/diversity/analysis/check_chimeras/"; usearch61 --maxaccepts 1 --consout "/home/qiime/diversity/analysis/check_chimeras/trim.fasta_consensus_with_abundance.fasta" --usersort --id 0.97 --sizeout --minseqlength 64 --wordlength 8 --uc "/home/qiime/diversity/analysis/check_chimeras/trim.fasta_consensus_with_abundance.uc" --cluster_smallmem "/home/qiime/diversity/analysis/trim.fasta" --maxrejects 8 --strand plus --log "/home/qiime/diversity/analysis/check_chimeras/trim.fasta_smallmem_clustered.log" > "/tmp/tmpHO5kHDhuaxiCdEO1ZTPT.txt" 2> "/tmp/tmpuuaVrSi6Xn7UsQa60jl0.txt" StdOut: StdErr: Segmentation fault 结果只在analysis文件夹下面生成了check_chimeras文件夹，但是check_chimeras文件夹里面没能生成chimeras.txt文件。 |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
kshdd
金虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 2819.8
- 散金: 763
- 红花: 12
- 沙发: 2
- 帖子: 1795
- 在线: 432.2小时
- 虫号: 3384976
- 注册: 2014-08-27
- 性别: GG
- 专业: 环境微生物学
2楼2018-09-02 00:54:55