talk @ metabarcoding.org
ecotag with custom reference library - Printable Version

+- talk @ metabarcoding.org (http://talk.metabarcoding.org)
+-- Forum: Bioinformatic softwares (http://talk.metabarcoding.org/forumdisplay.php?fid=1)
+--- Forum: OBITools (http://talk.metabarcoding.org/forumdisplay.php?fid=2)
+---- Forum: Using OBITools (http://talk.metabarcoding.org/forumdisplay.php?fid=7)
+---- Thread: ecotag with custom reference library (/showthread.php?tid=110)



ecotag with custom reference library - davon - 08-25-2017

Hi,
I've built a custom reference library from in house sequencing and converted it to an ecopcr database with obiconvert. Now, I'm trying to use it to assign taxonomy to my data. I downloaded taxonomic data as a taxdumptar.gz from NCBI. I'm getting errors when I run ecotag, it seems maybe my obiconvert didn't work right, as apparently i'm missing a .pdx file?
Thanks very much in advance for help! 
Davon

Obiconvert code: obiconvert  --fasta --ecopcrdb-output=customdb -t ./TAXO ./custom/custom_ref_lib.fasta

output files:

  24328144 24 Aug 17:04 customdb.adx
146202025 24 Aug 17:04 customdb.ndx
           385 24 Aug 17:03 customdb.rdx
  75063901 24 Aug 17:03 customdb.tdx
     130312 24 Aug 17:04 customdb_001.sdx

ecotag code:
ecotag -d customdb -R ./custom/custom_ref_lib.fasta ./EKB_run1/ebl_data/ebl.fasta ebl.taxonomy.fasta

ecotag error:

Reading binary taxonomy database...
 [INFO : Taxon alias file found] 
 [INFO : Local taxon file not found] 
Taxonomical tree read
[Errno 2] No such file or directory: 'customdb.pdx'
 [INFO : Preferred taxon name file not found]
 ok
Reading reference DB ...  : 750
Traceback (most recent call last):
  File "/Users/DC/Code/OBITools-1.2.10/export/bin/ecotag", line 345, in <module>
    assert seqid not in taxonlink
AssertionError


RE: ecotag with custom reference library - celine.mercier - 09-14-2017

Hi!

The "[Errno 2] No such file or directory: 'customdb.pdx'" message is actually a misleading information message, not an error, it means that you don't have a file with preferred taxon names in your database, which is optional.

The error when you run ecotag is actually due to the fact that there are (at least) 2 identical sequence identifiers in your fasta reference sequence file, which is not good because ids are used to identify sequences in a unique way.

To fix it, you can run:

obiannotate --uniq-id my_ref_seqs.fasta > my_ref_seqs_with_unique_ids.fasta

Céline