picrust2 run error

44 views
Skip to first unread message

er qw

unread,
May 12, 2025, 4:41:23 PMMay 12
to picrust-users
picrust2_pipeline.py -s rep-seqs_single/dna-sequences.fasta -i table_single/otu_table_single.tsv -o qiime_realsc -p 180

Warning - 15 input sequences aligned poorly to reference sequences (--min_align option specified a minimum proportion of 0.8 aligning to reference sequences). These input sequences will not be placed and will be excluded from downstream steps.

This is the set of poorly aligned input sequences to be excluded: 4474759, 338381, 887744, 4307094, 335132, 187469, 523483, 509621, 516366, 4294655, 4397579, 267005, 3910247, 184175, 312973




Traceback (most recent call last):
  File "/home/tools/anaconda3/envs/picrust2-sc/bin/picrust2_pipeline.py", line 312, in <module>
    main()
  File "/home/tools/anaconda3/envs/picrust2-sc/bin/picrust2_pipeline.py", line 272, in main
    func_outfiles, pathway_outfiles = full_pipeline_split(study_fasta=args.study_fasta,
  File "/home/tools/picrust2-2.6.2/picrust2/pipeline.py", line 306, in full_pipeline_split
    prune_tree(list(nsti_dom1_df.index.values)+ref_d1_seqs, out_tree_ref1, out_tree_ref1_red)
  File "/home/tools/picrust2-2.6.2/picrust2/util.py", line 496, in prune_tree
    tree.prune(names)
  File "/home/tools/anaconda3/envs/picrust2-sc/lib/python3.10/site-packages/ete3/coretype/tree.py", line 531, in prune
    to_keep = set(_translate_nodes(self, *nodes))
  File "/home/tools/anaconda3/envs/picrust2-sc/lib/python3.10/site-packages/ete3/coretype/tree.py", line 2607, in _translate_nodes
    raise TreeError("Invalid target node: "+str(n))
ete3.coretype.tree.TreeError: 'Invalid target node: 180227'

this is error message,but when i run  picrust2_pipeline_oldIMG.py -s rep-seqs_single/dna-sequences.fasta -i table_single/otu_table_single.tsv -o qiime_realsc -p 180  is OK, that's why?

Robyn Wright

unread,
May 12, 2025, 4:43:25 PMMay 12
to picrust-users
Hi there,

I think this could be an issue with your sequence names being only numbers, and one of the packages used by the new version therefore confusing them for internal nodes of the tree. Try adding "Sequence" to the beginning of each of your sequence names and then re-running the new version. 

Robyn

Reply all
Reply to author
Forward
0 new messages