This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Aligning fasta sequences against a reference genome with scaffolds

Hello,

I have a reference genome that consists only of scaffolds without the assignment of scaffolds to any chromosome. My reference genome is in fasta format.

I would like to extract those sequences that possibly belong to chromosome number 5 of my species. For that, I have a transcriptome de novo assembly from which I have obtained a list of fasta sequences that belong to chromosome 5.

What I would like to do now is to align my fasta sequences belonging to chromosome 5 obtained from transcriptome assembly against the reference genome scaffolds to obtain those scaffolds that belong to chromosome number 5.

My question is that what the best software is for this purpose?

Thank you very much for your help in n advance.

alignment

1 answer

There are many good choices among them I'd recommend lastal or lastz...

Note: despite the similar names the tools are completely unrelated

Log in to answer this question.