antismash error
2
0
Entering edit mode
5 days ago
szb6246 • 0

Hello, I have ran antismah in web and now trying to use through conda. I see that input file required is gbk. However even if we annotate fasta file with Rast or Prokka, the output is gbgff file and providing that always run into error such as

 Record c00014_NODE_14.. contains no genes and no genefinding tool specified


I have tried it with assemblies downloaded from NCBI but its same error. Can some one please guide me Command I am using is

antismash BP1199.gbff

antismash gbkfile • 251 views
0
Entering edit mode

That doesn't look like an error. Probably the NODE_14 is a short contig with no genes

2
Entering edit mode
4 days ago
kblin ▴ 30

You can run antiSMASH on FASTA inputs, antiSMASH will then use prodigal for gene finding. On GenBank inputs, antiSMASH will assume that you have already run gene finding yourself. As Andres mentioned, you're likely stumbling over a contig that's too small to have any genes on it. antiSMASH errors out on these when run on the command line to force you to pick a strategy for dealing with these (--genefinding-tool error). So if you just want to ignore contigs without gene calls, you could just run with --genefinding-tool none instead.

0
Entering edit mode
4 days ago
szb6246 • 0

I tried it with fasta file as well. I have assembled genomes using SPAdes. When I load them on web, I do get results but I need the output file gbk from the anti smash in command line to run the BIG-SCAPE.

How do I run antismash in command line with fasta file. Can you share the direct code please.