-
Notifications
You must be signed in to change notification settings - Fork 36
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
error standardized_DIAMOND_analysis_counter.py #77
Comments
Hi Arian, Sorry about that, let's see if we can figure out an answer! A couple things to check: First, are you using the standard bacterial RefSeq database that is available for distribution with SAMSA2, or is this a custom-built database? Second, it looks from this error like there's at least one entry where there's an ID but no description. Could you try re-running with the following code added to replace the existing line 138 of the standardized_DIAMOND_analysis_counter.py script?
The idea here is to see the offending line, as well as where in your reference database this line occurs (for easy fixing, if there's an issue with it). This would require you to re-run the script, although you don't need to re-run the whole master_script; you could just run the command:
(You'll need to correct the pathways and replace the variables with your output file and the RefSeq_db path.) Essentially, this looks to be an issue with the reference database; happy to help identify where it's occurring and what, specifically, is causing it. |
Hi Sam! By the way, I mentioned earlier that I downloaded Fungi protein sequences from Ensembl. Do you have any suggestions for other good sources or databases to use? I used NCBI refseq as well but I believe the sequence files were corrupted/weird, because the sequences had numbers or "_" in between which resulted an error with DIAMOND as well. Thank you again for your help. |
Hello again Sam, After replacing my old Fungi database with the one I downloaded from Ensembl, I've got list index out of range error in line 159. :/
here is the the python script from line 144 to 159:
I think the fungi dataset causes this issue, perhaps the annotations in the fasta files are not compatible with your script. Do you know where I could find a fungi dataset compatible with DIAMOND/your script? Thanks. BTW, here is an example the header/description of Fungi data which causes the error (I think?):
|
Any solution for this? |
Hello, I am having issues with standardized_DIAMOND_analysis_counter.py script
I am getting an IndexError in line 138
I've modified your master_script_for_sample_files.bash file and the error comes after STEP4 is DONE
command from your script:
STEP 5: AGGREGATING WITH ANALYSIS_COUNTER
error:
Now reading through the m8 results infile.
Analysis of /projects/bact.fun.unmapped.RefSeq_annotated complete.
Number of total lines: 574668
Number of unique sequences: 574668
Time elapsed: 1.8101940155 seconds.
then "Starting database analysis now." message pops and goes until
198M lines processed so far in 2025.08801007 seconds.
Then I get this error:
Traceback (most recent call last):
File "/projects/tools/samsa2/python_scripts/standardized_DIAMOND_analysis_counter.py", line 138, in
db_entry = db_entry[1][:-1]
IndexError: list index out of range
Here is an snapshot of your script from line 127 to 138
Thanks in advance
The text was updated successfully, but these errors were encountered: