Close genome comparison
Entering edit mode
21 months ago
A_heath ▴ 160

Hi all,

I have two fasta files corresponding to two close genomes (same species, and very close strains). However, the issue is that one of these genome is very fragmented and contains a lot of contigs (346 contigs) which makes the comparison very difficult. Ideally, I would like to find unique regions within the reference genome when compared to the fragmented one.

Most of the time, I use Mugsy to align genomes but in this case, the fragmented genome prevent this tool from performing a complete alignment. I also tried with Mauve or GSAlign but I can not identify unique genomic region on the reference genome.

Do you have any other solutions I could try to solve this, please ?


Mauve Mugsy genome-comparison GSAlign • 428 views
Entering edit mode

Have you tried reference guided scaffolding of the contigs? You might be able to place them and pad with Ns to get around some of the fragmentation, but you'll still be missing a fair bit of info from those regions I imagine.


Login before adding your answer.

Traffic: 2262 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6