Problems with choosing reference sequences

Skip to first unread message

Jul 5, 2017, 10:34:35 PM7/5/17
to Qiime 1 Forum

Dear all. 

I'm new in this field, having some problems with my first analysis.
 Could you help me to solve the problems of reference? 

I'm doing my works with 16S rRNA sequences. 

My problems are, 

1. Reference sequences: RDP

  - I want to use RDP reference sequence, but I have no idea which RDP data should I use for analysis; training set, unaligned, aligned.
  - I want to know the meaning of "RDP Training set". What is RDP Training set? Is this using when the taxonomy assignment method is 'RDP classifyer'? 

2. Error in Taxonomy assignment with RDP reference data

  I got error when '' with unaligned RDP data as reference sequence.

Clustering was: -i seqs_chimera_filtered.fna -o w_rdp_otus/ -r unaligned_rdp.fa -f -a --suppress_taxonomy_assignment

Then I tried to assignment with unaligned rdp reference sequence but I failed.

Here is the command: -i rep_set.fna -r Trimmed_reads/unaligned_rdp.fa -t Trimmed_reads/trainset16_022016.rdp/ -o unaligned_rdp_assignment

Traceback (most recent call last):
  File "/qiime/qiime/build/scripts-2.7/", line 417, in <module>
  File "/qiime/qiime/build/scripts-2.7/", line 394, in main
  File "/qiime/python/python2/lib/python2.7/site-packages/qiime-1.9.1.dev0-py2.7.egg/qiime/", line 1306, in __call__
    result = self._uc_to_assignments(app_result['ClusterFile'])
  File "/qiime/python/python2/lib/python2.7/site-packages/qiime-1.9.1.dev0-py2.7.egg/qiime/", line 1364, in _uc_to_assignments
    tax = self.id_to_taxonomy[subject_id].split(';')
KeyError: 'S003990863'

But, when I use greengenes as reference, there's no error. What should I do for solve this problem? 

I need some help to understand how to use the right reference sequence in each process. 
Thank you!

2017. 07. 06. THU.
- Hyojung Kim 

Se Jin Song

Jul 7, 2017, 2:32:11 PM7/7/17
to Qiime 1 Forum
Reply all
Reply to author
0 new messages