24小时热门版块排行榜    

查看: 2002  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
请填处理意见