diff --git a/lab_modules.py b/lab_modules.py index a56d13d..352190b 100755 --- a/lab_modules.py +++ b/lab_modules.py @@ -11,12 +11,12 @@ def main(): pass -def call_iqtree_with_SNP_aln(args, d_count): +def call_iqtree_with_SNP_aln(args, d_count, query): constants=','.join([str(d_count.get('A')), str(d_count.get('C')), str(d_count.get('G')), str(d_count.get('T'))]) - print ('constants',constants) - call(['nice', 'iqtree', '-s', args.out + '_trimmed_SNP.aln', + print ('constants', query, constants) + call(['nice', 'iqtree', '-s', query + '_trimmed_SNP.aln', '-nt', str(args.threads), '-m', 'GTR+G', '-bb', '1000', '-czb', '-fconst', constants])