This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Nucleotide alignments from mmseqs "tblastn"

I am searching a protein sequences against a nucleotide database ("tblastn") using mmseqs and need nucleotide sequence alignments as an output.

I start by the following (pg 23 of https://mmseqs.com/latest/userguide.pdf)

mmseqs search ecoli_proteins ecoli_genome alnDB tmp

and then get alignments as follows:

mmseqs convertalis ecoli_proteins ecoli_genome alnDB alnDB.html --format-mode 3

However, this gives me protein sequence alignments, when I really want the nucleotide base alignments (i.e. would be aligned codon-by-codon).

What is the proper way to get this output from mmseqs?

mmseqs2 tblastn mmseqs

0 answers

No answers yet.

Log in to answer this question.