Question: local alignment with a sets of sequences in a fasta file to a reference sequence.
gravatar for kanshenglong
4.4 years ago by
kanshenglong0 wrote:

I am in need of some help again.

I have a set of sequences (in a multifasta file), and I need to align them against a reference with local alinment. I have tried EMBL-EBI pairwise alignment tools, but it just calculate the local alignment of two sequences. And 'bwa mem' does a local alignment, it only aligns the best matching parts of my contigs while the remainder is clipped off, but I need the intact sequences. Other software that can be useful?

Thanks in advance.

alignment • 1.5k views
ADD COMMENTlink modified 4.4 years ago by Biostar ♦♦ 20 • written 4.4 years ago by kanshenglong0

Use bwa mem -Y....

ADD REPLYlink written 4.4 years ago by lh332k

Can I retrieve such alignment in a fasta format? Thanks a lot.

ADD REPLYlink written 4.4 years ago by kanshenglong0

I have not tried it before but bwasw option of bwa might work.

ADD REPLYlink written 4.4 years ago by Sej Modha4.7k

Look at the pair-wise alignment output formats from EMBOSS on this page (markxN series). I assume you are referring to EMBOSS tools on the EBI site. It may be that you just need to select -aformat markx0 as output for your pairwise alignment.

ADD REPLYlink modified 4.4 years ago • written 4.4 years ago by GenoMax94k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1608 users visited in the last hour