Biostar Beta. Not for public use.
Mummer alignment of two chromosomes
0
Entering edit mode
24 months ago
wjidea • 50
United States

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 • 1.8k views
ADD COMMENTlink
1
Entering edit mode
24 months ago
wjidea • 50
United States

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:

\gist bc53e3065f2962f19e82

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1