Skip to content
Snippets Groups Projects
Commit c8fafcb6 authored by Lukas Jelonek's avatar Lukas Jelonek
Browse files

Add max hits argument to run_ghostx.py

parent fb9ff632
No related branches found
No related tags found
Loading
......@@ -8,6 +8,7 @@ ghostx_tool = 'ghostx'
parser = argparse.ArgumentParser(description='Identify homologues in the swissprot database')
parser.add_argument('--fasta', '-f', required=True, help='A fasta file with aminoacid sequences')
parser.add_argument('--database', '-d', required=True, help='Database to search in')
parser.add_argument('--max_hits', '-m', type=int, default=20, help='Maximal number of reported hits')
parser.add_argument('--output', '-o', required=True, help='The result directory. Will contain info.json, results.tsv and enum_headers.tsv.')
args = parser.parse_args()
......@@ -22,4 +23,4 @@ print(" VERSION=\$(ghostx 2>&1 | grep -Eo '[[:digit:]]\\\\.[[:digit:]]\\\\.[[
toolconfig = '{\\\\"name\\\\": \\\\"ghostx\\\\", \\\\"database\\\\":' \
+ '\\\\"' + args.database + '\\\\"' + ', \\\\"version\\\\":\\\\"\$VERSION\\\\" }'
print(' echo ' + toolconfig + ' > ' + args.output + '/info.json')
print(' ' + ghostx_tool + " aln -d " + args.database + " -o " + args.output + "/results.tsv -i " + args.fasta)
print(' ' + ghostx_tool + " aln -v 5 -b " + str(args.max_hits) + " -d " + args.database + " -o " + args.output + "/results.tsv -i " + args.fasta)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment