Skip to content

Commit

Permalink
Add option to perform translated DIAMOND alignment
Browse files Browse the repository at this point in the history
  • Loading branch information
lauraluebbert authored Nov 12, 2024
1 parent d221f2e commit 3152c01
Showing 1 changed file with 10 additions and 1 deletion.
11 changes: 10 additions & 1 deletion gget/main.py
Original file line number Diff line number Diff line change
Expand Up @@ -456,7 +456,15 @@ def main():
type=str,
nargs="+",
required=True,
help="Reference sequences (str or list) or path to FASTA file containing reference sequences.",
help="Reference sequences (str or list) or path to FASTA file containing reference sequences. Add `-x` flag if reference sequences are amino acid sequences and query sequences are nucleotide sequences.",
)
parser_diamond.add_argument(
"-x",
"--translated",
default=False,
action="store_true",
required=False,
help="Perform translated alignment of nucleotide sequences to amino acid reference sequences.",
)
parser_diamond.add_argument(
"-db",
Expand Down Expand Up @@ -2692,6 +2700,7 @@ def main():
diamond_results = diamond(
query=args.query,
reference=args.reference,
translated=args.translated,
diamond_db=args.diamond_db,
sensitivity=args.sensitivity,
threads=args.threads,
Expand Down

0 comments on commit 3152c01

Please sign in to comment.