Mummer alignment of two chromosomes
1
0
Entering edit mode
5.6 years ago
wjidea ▴ 50

I was trying to align two full chromosome sequences (one contig for each chromosome) with Mummer v3.1, and I get mummer dot plot to visualize the alignment. These two chromosomes were actually very similar, however, the aligned region (red dots) were not shown on the diagonal. I am wondering is there a parameter I should change to let mummer put all the aligned region on the diagonal or I need to do some other manipulatons on my sequence file (fasta files).

The commands I used for this alignment was:

nucmer --maxgap=500 --mincluster=100 --prefix=ref_qry ref.fa qry.fa
show-coords -r ref_qry.delta > ref_qry.coords
show-aligns ref_qry.delta "ref" "qry" > ref_qry.aligns
delta-filter -q -r ref_qry.delta > ref_qry.filter
mummerplot --postscript ref_qry.filter -R ref.fa -Q qry.fa

 

Thank you for your help!

 

Mummer alignment • 2.6k views
ADD COMMENT
1
Entering edit mode
5.6 years ago
wjidea ▴ 50

I will answer to my own question then. I was trying to look for a flag I can modify to solve this problem. However, I didnt find a solution from Mummer itself. Therefore, I decided to use a stupid method combining blastn and python::faidx. I posted my solution to github gist, see below:

ADD COMMENT

Login before adding your answer.

Traffic: 1795 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