Nucleotide alignments from mmseqs "tblastn"
0
0
Entering edit mode
10 months ago
saladi ▴ 30

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 • 473 views
ADD COMMENT

Login before adding your answer.

Traffic: 1531 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6