Assessing Repeat Content Among Assembled Contigs
1
1
Entering edit mode
12.7 years ago
Abhi ★ 1.6k

Hi All

I am trying to asses and compare genomes assemblies for the same data set and suspect there can be some contigs assembled which are overlapping in the genome region they cover.

So wondering if I can detect and visualize these differences. In essence what I would like to do a given a multi-fasta file, do a pariwise comparison between all the contigs and see how many contigs have common regions from the genome.

PS: I have asked somewhat similar question but this time this is more specific.

http://biostar.stackexchange.com/questions/13370/comparing-multiple-genome-assemblies-with-reference

Thanks! -Abhi

fasta assembly • 2.3k views
ADD COMMENT
2
Entering edit mode
12.7 years ago
Torst ▴ 980

You should be able to use nucmer from the MuMMer package to do this, using your contigs multifasta file as the Query and as the Subject and using --nosimplify option. You can then visualize the dotplot using mummerplot to get a feel for the repetitiveness, then formally using the mummer output files.

ADD COMMENT

Login before adding your answer.

Traffic: 2892 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6