Alignment of genomes and get statistics about alignment

Hi,

I have 3 genomes: 1 reference and 2 query. I would like to know which of the 2 query is more similar to the reference.

I use to use nucmer for genome alignment but I can't find a way to get statistics about global percentage of identity ... or maybe I missed something ?

Do you have any recommendation for doing that ?

Thanks

genome

• 1,484 views • link

updated 4.8 years ago by manaswwm

• written 4.8 years ago by Picasa

whole genome alignment using nucmer

• 0 views • link

written 4.8 years ago by manaswwm

I know how to align genome with nucmer and got the plot etc. But I am looking for alignment statistics

• 0 views • link

written 4.8 years ago by Picasa

Maybe the answer from one of the user is useful for you, hence the comment

• 0 views • link

written 4.8 years ago by manaswwm

0 answers

No answers yet.

More posts like this

Detecting Inversions from Genome Alignments
written 3.0 years ago by saamhasan55 •

Hello, I have whole genome pairwise alignments generated using Nucmer. The alignment is a table file with chromosome number, starting and end positions of all …
BlastP results based query
written 4.6 years ago by Habibatariq68 •

which is better, if 1. query coverage is 100% and percentage identity of 69%, Max score: 675 2. query coverage is 100% and perc. identity …
Comparing two highly similar genomes
written 5.1 years ago by seta

Dear all, I have two highly similar mammalian genomes at the chromosome level, one is a wild species and another is a domestic species. For …
How to run nucmer without reference? (viruses metagenomics)
written 5.5 years ago by Arsenal •

Hi! I'm testing detection and abundance of viruses in metagenomics (viromics). In ([Roux et al, 2017][1]), The authors state the use of [nucmer][2] (from [mummer][3]) …
ggsearch: not returning all alignments for short sequences
written 5.7 years ago by felix.teufel •

Hi. I'm using ggsearch36 from the FASTA package to create a similarity matrix of protein sequences, using global identity as the similarity. I realized that …
How can I obtain column statistics from a multiple sequence alignment?
written 5.8 years ago by ej6474 •

I have aligned about ~2k viral genomes (highly similar in identity about 99%) using MAFFT to a reference genome to build a full MSA and …
Mummer for running multiple query files
written 6.7 years ago by AP •

Hello everyone, I am trying to use Mummer to differentiate core and accessory chromosomes in some genome assembly files I have. I used the following …
MUMmer overall identity & output
written 6.9 years ago by yp19

Hi guys! I want to use nucmer to assess how similar two species' assemblies are to one another. I used the command like this: nucmer …
Inversion detection between 2 genomes
written 8.2 years ago by Picasa

Hi all, I am comparing two very similar genomes and looking to see whether there are some part that are rearranged such as inversion. Do …
Why alignment in "blat" and "nucmer" are different?
written 8.5 years ago by Yuyayuya •

Hello everyone, I'm doing genome assembly using SOAPdenovo2-MUMmer3 pipeline. After scaffolds were built, I used `nucmer` to align each scaffold to chromosome based on the …