Skip to main content
Figure 4 | BMC Bioinformatics

Figure 4

From: Genome-wide DNA polymorphism analyses using VariScan

Figure 4

Application of the MRA analysis to the coalescent-simulated data set. The data contains 10 sequences of 2,000,000 bp each, and it was generated applying a per-site value of θ = 0.01. Upon this raw data set, we made two different levels of changes: i) two wide reductions in nucleotide diversity levels (g1: α = 1/3, β = 500,000; g2: α = 1/2, β = 500,000); and ii) 11 local valleys of reduced variability (v1: α = 1/4, β = 20,000; v2: α = 1/4, β = 15,000; v3: α = 1/4, β = 10,000; v4: α = 1/4, β = 5,000; v5: α = 1/4, β = 2,000; v6: α = 1/3, β = 20,000; v7: α = 1/3, β = 10,000; v8: α = 1/3, β = 5,000; v9: α = 1/2, β = 10,000; v10: α = 1/2, β = 5,000; v11: α = 1/2, β = 2,000). (a) nucleotide diversity profile obtained by SW using non-overlapping windows of 50 bp; (b) Signal reconstruction of low-frequency bands with information from 7 to 8 MRA levels, showing the location (in arrows) of 5 depleted-variation regions (v4–5, v8, v10–11; β ≤ 5,000). c) Signal reconstruction from 9 to 12 MRA levels, showing the location (in arrows) of 9 depleted-variation regions (v1–4, v6–10; 5,000 ≤ β ≤ 20,000). d) Signal reconstruction from 13 to 15 MRA levels, showing the location (in arrows) of the two broad areas with reduced levels of variation (g1–2; β = 500,000). The nucleotide sequence positions (X axis) are given in kb.

.

Back to article page