This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Find chromosomal location of a scaffold

For two closely related species, is there an easy way to align genome assemblies?

Species A is a chromosomal assembly (although it has gaps) and Species B is assembled into smaller scaffolds. I would like to figure out where Species B scaffolds lie on the Species A assembly in order to leverage annotations that are available for A but not B.

Thank you very much!

genome mapping
• 2,321 views • link
updated by h.mon h.mon • written by Jautis Jautis

1 answer

See this answer. You can also use bwa mem, minimap2, Blat, among others.

For annotation transfer, see RATT.

• 0 views • link
written by h.mon h.mon

Thanks for the response. However, the answer and methods you refer me to have to do with mapping reads to a genome. What I'm interested in doing is aligning fasta genomes to each other (or rather, finding the chromosome location scaffold X corresponds with).

• 0 views • link
written by Jautis Jautis

The question indeed is about mapping reads to a second reference genome, but no, the answer I linked is not about mapping reads to a genome, it is about aligning scaffolds to a complete genome:

Don't map the reads back to the Zebra Finch genome.

MUMmer is one of the programs commonly used for aligning contigs / scaffolds to reference genomes, and is a requirement for RATT, which is a software to transfer annotation from a reference (annotated) genome to an unannotated query genome.

edit: just to clarify, bwa mem, minimap2 and Blast are all also used for contig / scaffold mapping to a reference genome, with the requirement the reference genome should be close and reasonably similar to the query contig / scaffolds - I believe at most 10-15% difference, someone will correct me if I got the numbers wrong.

• 0 views • link
written by h.mon h.mon

Thanks for the clarification!

• 0 views • link
written by Jautis Jautis

I had indeed provided the wrong link, I linked to a comment to the answer I wanted to link. I've corrected it now.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this