katholt / srst2

Short Read Sequence Typing for Bacterial Pathogens
Other
123 stars 65 forks source link

grep gene db syntax error #92

Open IshwaryaMullaiSelvapandian opened 6 years ago

IshwaryaMullaiSelvapandian commented 6 years ago

Hi,

I have this issue at the end of running srst2 using resistance gene database. Can anyone please let me know why this error occurs and how to resolve it?

sh: -c: line 0: syntax error near unexpected token (' sh: -c: line 0:grep 69mef(A)mef(A)_31542 /testmount/efs/SharedData/Spn_Scripts_Reference-master/Spn_Reference_DB/ResFinder.fasta' sh: -c: line 0: syntax error near unexpected token (' sh: -c: line 0:grep 44msr(D)msr(D)_21944 /testmount/efs/SharedData/Spn_Scripts_Reference-master/Spn_Reference_DB/ResFinder.fasta' sh: -c: line 0: syntax error near unexpected token (' sh: -c: line 0:grep 70mef(A)__mef(A)_101225 /testmount/efs/SharedData/Spn_Scripts_Reference-master/Spn_Reference_DB/ResFinder.fasta'

Thanks, ishwarya

oliverdrechsel commented 4 years ago

Dear ishwarya,

this issue is pretty old, but maybe this is helpful for future users.

Your Fasta file probably contains '|' symbols in the fasta headers, e.g. ">MK268740.1|Bacillus cereus strain|"

As srst2 uses some bash logics the pipe symbol is interpreted as such and causes issues in a grep command.

Replace all non-alphanumeric letters in the fasta headers and srst2 should run better.