An Alternative To Blastx With Speed Similar To Bowtie
5
4
Entering edit mode
10.3 years ago
Tianyang Li ▴ 500

Hi,

Is there any fast aligner that aligns nucleotide sequences to protein databases?

I've used blastx but blastx is too slow and takes too much memory, I'm wondering if anything like bowtie exists.

Thanks!

alignment protein nucleotide blast bowtie • 3.7k views
3
Entering edit mode
10.3 years ago

Maybe Usearch might be useful.

2
Entering edit mode
10.3 years ago
Neilfws 49k

-t=prot -q=dnax


I've only used it for nucleotide-nucleotide searches, but it should be faster than BLAST. You just need sufficient RAM for the in-memory database indices.

0
Entering edit mode

Apparently blat can't do this. I got this when I ran it: d and q must both be either protein or dna

0
Entering edit mode

Really? Interesting, I've never tried and just assumed it had a blastx-like function.

0
Entering edit mode

The error message "d and q must both be either protein or dna" is not correct; it should begin "t and q". Furthermore, this error occurs when -q=dna, not -q=dnax as I suggested. The error in the latter case is "Unrecognized combination of target and query types."

An answer in this thread suggests simply switching the -t and -q options (i.e. -q=prot -t=dnax).

0
Entering edit mode

Blat with options -t=prot -q=dnax gives results similar to tBLASTn but not BLASTx. Be careful during the comparisons.

2
Entering edit mode
8.2 years ago
Dgg32 ▴ 90

RAPsearch 2 can do DNA vs protein or protein vs protein, very fast.

1
Entering edit mode
7.1 years ago
Len Trigg ★ 1.6k

The mapx command from Real Time Genomics is for this specific scenario. As well as being inherently very fast, it also supports multiple cores. It is included in our commercial RTG Metagenomics and RTG Core products, as well as the non-commercial version of RTG Core.

1
Entering edit mode
3.6 years ago
h.mon 34k

I know I am reviving a really old thread, but it is surprising no one has cited DIAMOND so far. It is not particularly light on resources, but it is incredibly fast.

0
Entering edit mode

All true, but unfortunately it stinks at doing blastX-like searches

0
Entering edit mode

really? besides protein-protein searches, I understood this was one of the principal objectives of the algorithm. What do you mean?