This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to align a genome fasta file from NCBI to a reference genome?

Hi,

I want to align a genome published on NCBI (scaffold level) to a reference genome (chromosome level). I though bwa mem could work, but it is stuck, running for too long without any error message. Do you know any software that could do this? Also, I would like to call SNPs after the alignment.

Thanks!

genome genetics alignment
• 2,195 views • link
updated by Antonio R. Franco Antonio R. Franco • written by ritafonsa • ritafonsa

mummer is probably a better option for large pairwise alignments.

• 0 views • link
written by Joe Joe

Please use a long read aligner like minimap2 or better still lastz (LINK) which is designed for chromosome sized pieces.

• 0 views • link
written by GenoMax GenoMax

I agree either mummer or minimap2 are good choices (don't know lastz). But as an aside, if the ploidy of the organism is anything but haploid you will lose a lot of variants that are lost when collapsing into a haploid reference sequence. It might be worth trying to find the original short or long read data they used to assembly the genome to call SNPs. They should be available in most bioproject accessions.

• 0 views • link
written by dthorbur dthorbur

hii can u pls tell me the file formats of the scaffold and chromosome level genomes please I'm also working on genome assemblies and I want to align a contig to a chromosome reference thanks

• 0 views • link
written by vasudhapai • vasudhapai

What's the purpose of aligning the contig to the chromosome reference? May change the tools that are suggested.

• 0 views • link
written by dthorbur dthorbur

1 answer

There is a version of a CUDA graphic accelerated version of lastz named SegAlign

• 0 views • link
written by Antonio R. Franco Antonio R. Franco

Log in to answer this question.

More posts like this