Error in DIAMOND_specific_organism_retriever.py

31 views
Skip to first unread message

Indah Oktaviani

unread,
Nov 12, 2022, 6:30:53 PM11/12/22
to SAMSA bioinformatics group
Hi, Sam. 

I've been running agregation step in samsa2 and want to run  DIAMOND_specific_organism_retriever.py but I meet this error message:

Database is read and set up, moving on to the infile...
Traceback (most recent call last):
  File "DIAMOND_specific_organism_retriever.py", line 148, in <module>
    infile = open (infile_name, "r")
IOError: [Errno 2] No such file or directory: 'view_RMTA1a.m8'


Here I attach the command I used: python DIAMOND_specific_organism_retriever.py -I view_RMTA1a.m8 -D nr1_1653.faa -SO 

Could you help me to solve this? any suggestion?

Best,
IO

Sam Westreich

unread,
Nov 12, 2022, 6:53:29 PM11/12/22
to Indah Oktaviani, SAMSA bioinformatics group
Hi Indah,

This looks pretty straightforward; can you confirm that the infile, 'view_RMTA1a.m8', is in the current directory where you're running this command?  If not, you'll need to supply a path to the file with the -I command, instead of just the file name.

Additionally, in your command, you don't have any organism name specified after the -SO flag.  You'll probably want to do so, e.g. 

$ python DIAMOND_specific_organism_retriever.py -I view_RMTA1a.m8 -D nr1_1653.faa -SO Prevotella

Best,
Sam

--
You received this message because you are subscribed to the Google Groups "SAMSA bioinformatics group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to samsa-bioinformatic...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/samsa-bioinformatics-group/0c8a48aa-5701-4ff3-890a-112d8bee38aan%40googlegroups.com.


--
Sam Westreich, PMP, PhD
Microbiome Scientist, DNAnexus, 
Reply all
Reply to author
Forward
0 new messages