Bash not finding file
1
0
Entering edit mode
5.5 years ago
mac03pat ▴ 30

I'm running Debian 9 as VM from serial console on Azure.

Link to files I'm attempting to input: https://drive.google.com/open?id=1aIIjWO70clU8u4_gV2X_8c17HmEnqjr_

I was given some code from the help people at antiSMASH. It should iterate through genbank files in a folder and process them with the antiSMASH package:

for infile in inputs/*.gbk; do
    antismash $infile --taxon fungi --input-type nucl --knownclusterblast
done

I tried running it while at the directory containing the files, and I tried running it while at the directory containing the folder "inputs" which contains the .gb files:

(antismash) macpat@Debian9:~/inputs$ for infile in inputs/*.gb; do
> antismash $infile --taxon fungi --input-type nucl --knownclusterblast
> done
ERROR   25/05 21:53:00   No sequence file found at 'inputs/*.gb'

This is the email I was sent by the antiSMASH people:

Dear a ntiSMASH user,

For running antiSMASH on many input files, I usually write a loop in bash, like so:

for infile in inputs/*.gbk; do
    antismash $infile --your --other-options --here
done

Assuming your input files are in GenBank format and located in a subdirectory of your current directory called "input", antiSMASH will run on all your input files sequentially. I'm aware you asked for "all at once", but as antiSMASH has pretty decent CPU and memory requirements, especially when running ClusterBlast, I would not recommend that.

Best regards, Kai

bash linux debian9 • 1.9k views
ADD COMMENT
1
Entering edit mode

Can you show us what you get when you do: for infile in ~/inputs/*.gb; do echo $infile; done?

ADD REPLY
0
Entering edit mode
(antismash) macpat@Debian9:~$ for infile in ~/inputs/*.gb; do echo $infile; done
/home/macpat/inputs/DQ660910.gbor infile in ~/inputs/*.gb; do echo $infile; done
/home/macpat/inputs/EU872212.gb
/home/macpat/inputs/GU930713.gb
/home/macpat/inputs/GU930714.gb
/home/macpat/inputs/HM180407.gb
/home/macpat/inputs/HM180409.gb
/home/macpat/inputs/HQ823618.gb
/home/macpat/inputs/HQ823619.gb
/home/macpat/inputs/HQ823620.gb
/home/macpat/inputs/HQ823621.gb
/home/macpat/inputs/JN408682.gb
/home/macpat/inputs/JQ340775.gb
/home/macpat/inputs/JX067626.gb
/home/macpat/inputs/JX067627.gb
/home/macpat/inputs/JX232185.gb
/home/macpat/inputs/JX232186.gb
/home/macpat/inputs/JX232187.gb
/home/macpat/inputs/JX232188.gb
/home/macpat/inputs/KJ501919.gb
/home/macpat/inputs/MG777489.gb
/home/macpat/inputs/MG777490.gb
/home/macpat/inputs/MG777491.gb
/home/macpat/inputs/MG777492.gb
/home/macpat/inputs/MG777493.gb
/home/macpat/inputs/MG777494.gb
/home/macpat/inputs/MG777495.gb
/home/macpat/inputs/MG777496.gb
/home/macpat/inputs/MG777497.gb
/home/macpat/inputs/MG777498.gb
/home/macpat/inputs/MG777499.gb
/home/macpat/inputs/MG777500.gb
/home/macpat/inputs/MG777501.gb
/home/macpat/inputs/MG777502.gb
ADD REPLY
4
Entering edit mode
5.5 years ago
AK ★ 2.2k

I see. Try again with:

for infile in ~/inputs/*.gb; do antismash $infile --taxon fungi --input-type nucl --knownclusterblast; done
ADD COMMENT
0
Entering edit mode

Seems to have worked, thank you!. I beleive now I'm running into problems with the package itself. The online version returns results for the same file that I'm trying to run on the downloaded package version of the program. However, the downloaded version is returning me this error:

WARNING 25/05 23:06:49   Error parsing best knowncluster hit; knownclusters array = []. Possibly no significant hits to known biosynthetic gene clusters.
ADD REPLY
0
Entering edit mode

You're welcome. This is already a new question, please consider posting to a new one.

ADD REPLY
0
Entering edit mode

Not sure why the error appeared but files were output! Good to go, thank you.

ADD REPLY

Login before adding your answer.

Traffic: 2021 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6