Short sequences against short sequences alignment
1
0
Entering edit mode
8.8 years ago

Hi,

I have thousand short sequences (~30 bp) in fasta format and I want to compare all of them against a database of short sequences (also around 30 bp), which are also in a fasta format.

I know blastn and USEARCH can probably handle such comparisons, but have you some advice?

Thank you

sequence alignment • 1.9k views
ADD COMMENT
0
Entering edit mode

A little up!

Any advice?

ADD REPLY
0
Entering edit mode

Can you provide some additional information?

Are these sequences homologous or not so (if not what % identity on average)? Are you ok with allowing gaps in alignment and do you want to the alignments to be global (end-to-end)?

Depending on your answer, blat may work. You could also give BBMap a try.

ADD REPLY
2
Entering edit mode
8.8 years ago

Why not just try blastn? It sounds like it shouldn't take a long time to build a BLAST database and query it.

ADD COMMENT
1
Entering edit mode

Yeah, blastn -task blastn-short

ADD REPLY

Login before adding your answer.

Traffic: 1414 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