24小时热门版块排行榜    

查看: 2003  |  回复: 1

tjneuron

捐助贵宾 (小有名气)

[求助] 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文件。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kshdd

金虫 (著名写手)

请问,这个silva.gold.ng.fasta在哪找到的,具体的链接能不能发我
2楼2018-09-02 00:54:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tjneuron 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602 289分求调剂 +5 WWW西西弗斯 2026-03-24 5/250 2026-03-24 18:51 by jhhcooi
[考研] 340求调剂 +3 话梅糖111 2026-03-24 3/150 2026-03-24 17:47 by hyzs6688
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +14 吃吃吃才有意义 2026-03-19 15/750 2026-03-24 17:28 by dick_runner
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 316求调剂 +7 梁茜雯 2026-03-19 7/350 2026-03-23 16:21 by lingjue
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
请填处理意见