diff --git a/AAI.py b/AAI.py new file mode 100755 index 0000000..bec794d --- /dev/null +++ b/AAI.py @@ -0,0 +1,114 @@ +#!/usr/bin/env python + + +import os +from argparse import ArgumentParser +from glob import glob +from numpy import mean, std +from subprocess import Popen +from sys import exit +from time import strftime + +def parseArgs(): + parser = ArgumentParser(description='Computes the average amino acid ' + 'identity (AAI) between two protein sets', add_help=False) + req = parser.add_argument_group('Required') + req.add_argument('-1', '--set1', required=True, metavar='FILE', + help='first input FastA sequence file') + req.add_argument('-2', '--set2', required=True, metavar='FILE', + help='second input FastA sequence file') + opt = parser.add_argument_group('Optional') + opt.add_argument('-c', '--cpus', type=str, metavar='INT', + default='1', help='number of CPUs [1]') + opt.add_argument('-f', '--fraction', type=float, metavar='FLOAT', + default=70.0, help='minimum alignment length percentage [70.0]') + opt.add_argument('-h', '--help', action='help', + help='show this help message and exit') + opt.add_argument('-i', '--identity', type=float, metavar='FLOAT', + default=30.0, help='minimum percent identity [30.0]') + opt.add_argument('-l', '--length', type=int, metavar='INT', + default=0, help='minimum alignment character length (sum of all ' + 'aligned segments and all gaps) [0]') + opt.add_argument('-o', '--outpath', metavar='PATH', + default=None, help='output directory [./AAI--_