MUMmer plot visualization
0
0
Entering edit mode
9 months ago
rthapa ▴ 50

Hi, I am using MUMmer plot to compare between a de novo assembly with reference genome. The percent identity between two genomes is more than 99% but when I plot the two genomes with MUMmer plot, the plot (https://ibb.co/q08xsZw) doesn't look like that. Does anyone have any idea? Thanks

enter image description here

MUMmer plot • 831 views
ADD COMMENT
0
Entering edit mode

Have you tried to reverse complement one of the sequences before trying this again? I don't recollect if MUMmer tries that.

ADD REPLY
0
Entering edit mode

No, I did't try the reverse complement. Do you have any suggestion on tools that we can use for getting reverse complement? The alignment looks like https://ibb.co/cL6ydYr with mauve. enter image description here

ADD REPLY
0
Entering edit mode

You can use reformat.sh from BBMap suite to do the reverse complement.

reformat.sh in=seq.fa out=rc.fa rcomp=t
ADD REPLY
0
Entering edit mode

It looks like the reverse complement is better to visualize. It seems like there are two inversions in the assembled genome https://ibb.co/k5MqYqZ.

enter image description here Is the red bar present in the right corner duplicated region in the reference genome? Do you have any idea if I need to get reverse complement to call the structural variants with Mummer or I can just use the de novo assembly?

ADD REPLY
0
Entering edit mode

Those blue segments represent inversions and translocations (which are apparent in your mauve plots too).

ADD REPLY
0
Entering edit mode

Thank you. But when I check the mummer results for structural variation, there are four inversion in the assembly compared to reference genome. I wonder why mummer plot is showing only two inversions.

tig00000001 INV 228534  226344  -2189
tig00000001 GAP 745271  742507  -2763   -232    -2531
tig00000001 GAP 765321  764991  -329    -151    -178
tig00000001 INV 872150  866517  -5632
tig00000001 GAP 1199745 1199745 1   -47 48
tig00000001 GAP 1900829 1900827 -1  -13978  13977
tig00000001 GAP 2499469 2498978 -490    -336    -154
tig00000001 DUP 2561758 2562745 988
tig00000001 BRK 2562746 2563397 652
tig00000001 DUP 2563398 2563566 169
tig00000001 BRK 2563567 2563570 4
tig00000001 INV 2769150 2763519 -5630
tig00000001 INV 3375024 3372831 -2192
tig00000001 JMP 3688215 3688214 0

``

ADD REPLY

Login before adding your answer.

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