Differences

This shows you the differences between two versions of the page.

--- whole-genome_alignment [2018/07/03 16:43] – hyjeong
+++ whole-genome_alignment [2021/03/17 13:09] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
 ====== Whole-genome alignment (MUMmer) ======
-[[http://mummer.sourceforge.net/|MUMmer]]는 suffix tree algorithm을 이용하여 genome sequence 수준의 alignment를 생성하는 도구이다. 발표된지는 매우 오래되었으나 매우 완성도가 높아서 여전히 많은 프로그램에서 활용되고 있다. MUMmer는 alignment를 시각화하는 도구는 탑재하고 있지 않으므로 [[http://www.gnuplot.info/|gnuplot]]에 의존한다. 핵심 프로그램인 mummer는 짧은 exact match(20 bp)의 단위를 빠르게 탐색하는 역할을 하며, 이들이 어떤 기준을 충족하여 모여있으면 이를 확장해 나가는 alignment generator script(예: nucmer)와 기타 utility가 패키지에 포함된다.
+[[http://mummer.sourceforge.net/|MUMmer]]는 suffix tree algorithm을 이용하여 genome sequence 수준의 alignment를 생성하는 도구이다. 발표된지는 매우 오래되었으나 매우 완성도가 높아서 여전히 많은 프로그램에서 활용되고 있다. 기본 컨셉트는 20 nt의 짧은 exact match(-l 옵션으로 조정 가능)을 reference와 query genome 모두에서 찾아내는 것이다. MUMmer는 alignment를 시각화하는 도구는 탑재하고 있지 않으므로 [[http://www.gnuplot.info/|gnuplot]]에 의존한다. 핵심 프로그램인 mummer는 짧은 exact match(20 bp)의 단위를 빠르게 탐색하는 역할을 하며, 이들이 어떤 기준을 충족하여 모여있으면 이를 확장해 나가는 alignment generator script(예: nucmer)와 기타 utility가 패키지에 포함된다.
+bp의 match는 genome에 반복하여 나타날 수 있다. 이론상으로는 하나의 20 bp는 4^20 = 1.1x10^12, 즉 약 1.1 테라염기쌍마다 한번씩 나타날 것이겠지만 말이다. 반복하여 나타나는 match의 처리는 다음의 옵션으로 조정한다.
+  \-mum           compute maximal matches that are unique in both sequences
+  -mumcand       same as -mumreference
+  -mumreference  compute maximal matches that are unique in the reference-
+                 sequence but not necessarily in the query-sequence (default)
+  -maxmatch      compute all maximal matches regardless of their uniqueness
 ===== MUMmer로 할 수 있는 작업 =====
   - Aligning two finished sequences (highly similar with/without rearrangements, fairly similar/dissimilar sequences)
@@ Line 36: / Line 45: @@
 ==== [2] Nucmer 사용 ====
+MUMmer를 가장 편하게 사용하는 것은 alignment generator인 nucmer Perl script를 쓰는 것이다. 일정 기준에 맞게 인접한 maximal exact match는 하나로 병합할 수 있는데, 이는 -g(maxgap; default 90)과 -l(min length; default 20)으로 조정 가능하다. show-coords는 alignment block의 위치를 표시한다. -c는 결과에 percent coverage 포함, -l은 결과에 서열 길이 포함, -r은 reference coordinate에 맞춤을 의미한다.
   $ nucmer -p nucmer_test E681.fa M1.fa
   $ mummerplot nucmer_test.delta