Question: Mummer alignment of two chromosomes
0
gravatar for wjidea
5.0 years ago by
wjidea50
United States
wjidea50 wrote:

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!

 

alignment mummer • 2.4k views
ADD COMMENTlink modified 5.0 years ago • written 5.0 years ago by wjidea50
1
gravatar for wjidea
5.0 years ago by
wjidea50
United States
wjidea50 wrote:

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 COMMENTlink modified 11 months ago by _r_am31k • written 5.0 years ago by wjidea50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1998 users visited in the last hour