Assign-Taxonomy-with-BLAST icon indicating copy to clipboard operation
Assign-Taxonomy-with-BLAST copied to clipboard

Published 20 hours ago • verified publisher icon Joseph7e

using script on custom database

Open javiercnav opened this issue 6 years ago • 1 comments

Hi Joe, I am trying to use the script to assign taxonomy to some 18S sequences. For this, I have generated my own blast database from a blast search of first 10 hits for each of my sequences and tax file that connects my blast database sequences with their taxonomy. My tax database has been curated to have 7 levels only k,p,c,o,f,g,s...

I tried to run the script and I get this message: Traceback (most recent call last): File "taxonomy_assignment_BLAST_V2.py", line 340, in best_level_taxonomy, blast_percent = Assign_Taxonomy(current_query, current_best_hits) File "taxonomy_assignment_BLAST_V2.py", line 288, in Assign_Taxonomy s.add(j[i]) IndexError: list index out of range

Could it be that the script is not recognizing the categories in my tax file?

blast_output_custom_format.txt log_file.txt

javiercnav avatar Feb 06 '19 23:02 javiercnav