library(dartRverse)
library(tidyverse)
library(RColorBrewer)
library(patchwork)
library(scales)
library(ggplot2)
library(gridExtra)
library(grid)
library(raster)
library(otuSummary)
13 Combining Genomic Resources
Session Presenters
Introduction
Chromosomal rearrangements are often associated with local adaptation and speciation because they suppress recombination, and as a result, rearrangements have been implicated in disrupting geneflow. Although there is strong evidence to suggest that chromosome rearrangements are a factor in genetic isolation of divergent populations, the underlying mechanism remains elusive. Here we applied an integrative cytogenetics and genomics approach testing whether chromosomal rearrangements are the initial process, or a consequence, of population divergence in the dwarf goanna, Varanus acanthurus. Specifically, we tested whether chromosome rearrangements are indicators of genetic barriers that can be used to identify divergent populations by looking at geneflow within and between populations with rearrangements. We found that geneflow was present between individuals with chromosome rearrangements within populations, but there was no geneflow between populations that had similar chromosome rearrangements. Moreover, we identified a correlation between reduced genetic variation in populations with a higher frequency of homozygous submetacentric individuals. These findings suggest that chromosomal rearrangements were widespread prior to divergence and because we found populations with higher frequencies of submetacentric chromosomes were associated with lower genetic diversity, this could indicate that polymorphisms within populations are early indicators of genetic drift.
Required packages
make sure you have the packages installed, see Install dartRverse
Geneflow between Karyotypes
# setup -----------------------
## load ---
<- readRDS("./data/spiny_gl1.rds")
vac.glch2
# missing karyotype MM in East population
@other$ind.metrics$karyo[grep("AA095619", vac.glch2@ind.names)] <- "MM"
vac.glch2
## overview --------------------------------------------------------------------
vac.glch2
/// GENLIGHT OBJECT /////////
// 233 genotypes, 301,738 binary SNPs, size: 302 Mb
19405083 (27.6 %) missing data
// Basic content
@gen: list of 233 SNPbin
@ploidy: ploidy of each individual (range: 2-2)
// Optional content
@ind.names: 233 individual labels
@loc.names: 301738 locus labels
@loc.all: 301738 alleles
@pop: population of each individual (group size range: 1-36)
@other: a list containing: loc.metrics ind.metrics latlon loc.metrics.flags history verbose
pop(vac.glch2) %>% table
.
b baritji bb c cbb dr dw ea East ekim
4 15 1 1 15 1 4 5 10 4
gd kim mi nb North q rg sd South storri
3 30 4 1 9 1 9 8 9 5
vd wa wam West wm
3 36 32 21 2
table(vac.glch2@other$ind.metrics$PCR1)
f m
172 23 38
@other$ind.metrics$sex = NULL
vac.glch2
= toupper(as.character(vac.glch2@other$ind.metrics$PCR1))
sex # sex[sex==""] = "M" #TODO FIX na
@other$ind.metrics$sex <- sex
vac.glch2
#gl.smearplot(vac.glch2, group.pop = FALSE)
## filtering -------------------------------------------------------------------
<-gl.keep.pop(vac.glch2, pop.list = c("North","South","East","West")) vac.pops
Starting gl.keep.pop
Processing genlight object with SNP data
Checking for presence of nominated populations
Retaining only populations North, South, East, West
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.pop
#gl.smearplot(vac.pops, group.pop = TRUE)
<- gl.filter.sexlinked(vac.glch2, system = "zw")
sex table(sex$w.linked@other$loc.metrics$Chrom_Varanus_JD_infer_w)
scaf_10323 scaf_1064 scaf_1092 scaf_1121 scaf_11303 scaf_1133
299862 2 14 3 19 1 1
scaf_13939 scaf_14539 scaf_15833 scaf_159 scaf_164 scaf_16462 scaf_16941
5 5 1 12 28 1 6
scaf_17264 scaf_17378 scaf_175 scaf_17533 scaf_176 scaf_1768 scaf_177
2 3 63 2 28 4 17
scaf_17745 scaf_181 scaf_18976 scaf_19496 scaf_19546 scaf_19676 scaf_19678
1 23 7 5 12 1 15
scaf_19686 scaf_19710 scaf_20029 scaf_20122 scaf_20137 scaf_20324 scaf_2063
1 6 20 3 1 9 1
scaf_2296 scaf_24228 scaf_2674 scaf_2686 scaf_28831 scaf_29153 scaf_29340
1 3 4 2 9 1 7
scaf_30879 scaf_31396 scaf_32040 scaf_33358 scaf_33490 scaf_33621 scaf_33754
3 4 3 4 4 3 9
scaf_33755 scaf_33815 scaf_33834 scaf_33937 scaf_33995 scaf_34187 scaf_34257
10 6 3 5 3 1 5
scaf_34264 scaf_34293 scaf_34296 scaf_34780 scaf_35217 scaf_37485 scaf_37914
3 59 13 1 2 34 4
scaf_38564 scaf_39169 scaf_39480 scaf_39564 scaf_39585 scaf_39615 scaf_40621
10 29 1 4 4 1 1
scaf_40908 scaf_40937 scaf_41277 scaf_41356 scaf_41363 scaf_41402 scaf_4141
2 8 9 1 11 17 6
scaf_41492 scaf_41662 scaf_42222 scaf_42472 scaf_42597 scaf_42603 scaf_42605
25 6 1 16 9 18 8
scaf_42608 scaf_42630 scaf_43039 scaf_43122 scaf_43126 scaf_43132 scaf_43207
13 47 15 3 32 5 39
scaf_43209 scaf_43227 scaf_43364 scaf_43405 scaf_43511 scaf_43848 scaf_43901
13 7 20 16 3 2 14
scaf_43906 scaf_43915 scaf_43977 scaf_44080 scaf_44104 scaf_44191 scaf_44251
19 53 89 70 91 2 1
scaf_44262 scaf_44332 scaf_44349 scaf_44350 scaf_44435 scaf_44535 scaf_44569
75 8 8 30 13 52 8
scaf_591 scaf_619 scaf_6829 scaf_7848 scaf_8162 scaf_8185 scaf_9006
33 25 1 21 3 5 6
scaf_9015 scaf_9025 scaf_9028 scaf_9029 scaf_9033 scaf_9043 scaf_939
6 5 22 1 1 2 45
scaf_10851 scaf_11034 scaf_18166 scaf_20103 scaf_22094 scaf_2413 scaf_30511
2 4 1 2 2 6 2
scaf_33868 scaf_33956 scaf_34282 scaf_39568 scaf_40828 scaf_41352 scaf_41577
2 2 3 2 4 5 6
scaf_42194 scaf_42206 scaf_42407 scaf_42590 scaf_43134 scaf_43255 scaf_43347
8 2 1 2 11 1 1
scaf_835 scaf_9036 scaf_9133 scaf_944 scaf_1049 scaf_10750 scaf_10952
2 6 8 4 3 1 2
scaf_14926 scaf_16087 scaf_16565 scaf_17565 scaf_18106 scaf_18615 scaf_19081
2 1 5 2 4 1 1
scaf_19305 scaf_19657 scaf_19658 scaf_19659 scaf_19794 scaf_19885 scaf_20048
1 1 1 1 1 3 3
scaf_20059 scaf_20116 scaf_20213 scaf_20222 scaf_20326 scaf_23061 scaf_2445
2 1 1 1 1 2 1
scaf_2667 scaf_28826 scaf_29810 scaf_30302 scaf_32488 scaf_32653 scaf_32815
1 1 3 2 2 1 1
scaf_33422 scaf_33451 scaf_33595 scaf_33665 scaf_33841 scaf_33845 scaf_34167
2 3 1 2 1 10 1
scaf_34261 scaf_34265 scaf_34360 scaf_35141 scaf_37830 scaf_37837 scaf_37982
1 2 2 1 1 1 1
scaf_38014 scaf_38163 scaf_38481 scaf_39075 scaf_39298 scaf_39668 scaf_39916
1 1 1 1 1 5 4
scaf_40887 scaf_41264 scaf_41274 scaf_41390 scaf_4168 scaf_41789 scaf_42039
4 4 5 1 1 4 3
scaf_42444 scaf_42617 scaf_43226 scaf_43673 scaf_44264 scaf_44620 scaf_44621
2 1 1 1 2 1 2
scaf_4750 scaf_4773 scaf_526 scaf_574 scaf_8262 scaf_8885 scaf_8931
3 1 1 1 3 3 3
scaf_8974 scaf_901 scaf_9081 scaf_950
1 1 1 2
## start ----------
<- gl.filter.callrate(vac.pops, method = "loc", threshold = 0.80) vac.gl1
Starting gl.filter.callrate
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Warning: Data may include monomorphic loci in call rate
calculations for filtering
Recalculating Call Rate
Removing loci based on Call Rate, threshold = 0.8
Completed: gl.filter.callrate
<- gl.filter.callrate(vac.gl1, method = "ind", threshold = 0.75) vac.gl1
Starting gl.filter.callrate
Processing genlight object with SNP data
Warning: Data may include monomorphic loci in call rate
calculations for filtering
Recalculating Call Rate
Removing individuals based on Call Rate, threshold = 0.75
Note: Locus metrics not recalculated
Note: Resultant monomorphic loci not deleted
Completed: gl.filter.callrate
<- gl.filter.reproducibility(vac.gl1, t = 0.99) vac.gl1
Starting gl.filter.reproducibility
Processing genlight object with SNP data
Removing loci with repeatability less than 0.99
Completed: gl.filter.reproducibility
<- gl.filter.monomorphs(vac.gl1, v = 0)
vac.gl1
## colours -----------------------------------
<- c(East = "#0571b0", North = "#404040",
popcolors South = "#5e3c99", West = "#e66101")
<- c("MM" = "red", "MA" = "black", "AA" = "blue") karyocolors
## pcoa ------------------------------------------------------------------------
<- gl.pcoa(vac.gl1, nfactors = 5) pc
Starting gl.pcoa
Processing genlight object with SNP data
Performing a PCA, individuals as entities, loci as attributes, SNP genotype as state
Starting gl.colors
Selected color type 2
Completed: gl.colors
Completed: gl.pcoa
## ibd data --------------------------------------------------------------------
# pop
<- complete.cases(vac.gl1@other$latlon) # missing locations
indexXY <- gl.keep.ind(vac.gl1,ind.list = indNames(vac.gl1)[indexXY]) vac.gl17
Starting gl.keep.ind
Processing genlight object with SNP data
Deleting all but the listed individuals ABTC11930_aca, ABTC11932_aca, ABTC72881_aca, ABTC8929_aca, ABTC8938_aca, AA095616, AA093614, AA095550, AA095620, AA095553, AA095548, AA095617, AA093615, AA095626, AA095636, AA095546, AA096028, AA095562, AA095623, AA095635, AA095547, AA096030, AA095629, AA095624, AA095637, AA095545, AA095568, AA095625, AA095542, AA095570, AA095549, AA095551, AA095543, AA096029, AA095627, AA095540, AA095544, AA095634, AA095569, AA095539a
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.ind
<- runif(nInd(vac.gl17))/10000 # unique individual locations
latAdjust @other$latlon$lat <- vac.gl17@other$latlon$lat+ latAdjust
vac.gl17
# karyo
<- gl.drop.pop(vac.gl17, as.pop = "karyo", pop.list = "") vac.gl18
Starting gl.drop.pop
Processing genlight object with SNP data
Temporarily assigning karyo as population
Checking for presence of nominated populations, deleting them
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Restoring population assignments to initial state
Completed: gl.drop.pop
pop(vac.gl18) <- vac.gl18@other$ind.metrics$karyo
# figure 3 ---------------------------------------------------------------------
<- levels(vac.gl1@pop)
popi <- vac.gl1
vac.edit # all East population assumed to be MM
@other$ind.metrics$karyo[pop(vac.gl1)=="East"] <- "MM" vac.edit
Within population PCoA and IBD
for(i in 1:4){
#pcoa
<- gl.keep.pop(vac.edit, pop.list = popi[i], as.pop = "pop")
vac.popi <- gl.filter.monomorphs(vac.popi, verbose = 0)
vacp pop(vacp) <- vacp@other$ind.metrics$karyo
<- gl.keep.pop(vacp, as.pop = "karyo", pop.list = names(karyocolors), verbose = 0)
vacpK <- gl.pcoa(vacpK, nfactors = 5)
pca gl.pcoa.plot(pca, vacpK)#plot.pcoa(pca, vacpK, xpc = 1, ypc = 2, popColors = karyocolors)
<- list.files(tempdir())[grep("Plot", list.files(tempdir()))]
plotfiles unlink(paste0(tempdir(), "/", plotfiles))
#ibd
gl.ibd(vacpK, distance = "propShared", paircols = "pop",
Dgeo_trans = "Dgeo")
}
Starting gl.keep.pop
Processing genlight object with SNP data
Temporarily assigning pop as population
Checking for presence of nominated populations
Retaining only populations East
Locus metrics not recalculated
Resetting population assignments to initial state
Completed: gl.keep.pop
Warning: Listed population MA not present in the dataset -- ignored
Warning: Listed population AA not present in the dataset -- ignored
Starting gl.pcoa
Processing genlight object with SNP data
Performing a PCA, individuals as entities, loci as attributes, SNP genotype as state
Starting gl.colors
Selected color type 2
Completed: gl.colors
Completed: gl.pcoa
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Warning: Removed 9 rows containing non-finite outside the scale range
(`stat_smooth()`).
Warning: Removed 9 rows containing missing values or values outside the scale range
(`geom_point()`).
Removed 9 rows containing missing values or values outside the scale range
(`geom_point()`).
Removed 9 rows containing missing values or values outside the scale range
(`geom_point()`).
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.6174
Significance: 0.003
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.290 0.372 0.431 0.507
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Starting gl.keep.pop
Processing genlight object with SNP data
Temporarily assigning pop as population
Checking for presence of nominated populations
Retaining only populations North
Locus metrics not recalculated
Resetting population assignments to initial state
Completed: gl.keep.pop
Warning: Listed population MM not present in the dataset -- ignored
Starting gl.pcoa
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Performing a PCA, individuals as entities, loci as attributes, SNP genotype as state
Starting gl.colors
Selected color type 2
Completed: gl.colors
Completed: gl.pcoa
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.3604
Significance: 0.025
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.214 0.286 0.353 0.411
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Starting gl.keep.pop
Processing genlight object with SNP data
Temporarily assigning pop as population
Checking for presence of nominated populations
Retaining only populations South
Locus metrics not recalculated
Resetting population assignments to initial state
Completed: gl.keep.pop
Warning: Listed population AA not present in the dataset -- ignored
Starting gl.pcoa
Processing genlight object with SNP data
Performing a PCA, individuals as entities, loci as attributes, SNP genotype as state
Starting gl.colors
Selected color type 2
Completed: gl.colors
Completed: gl.pcoa
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.643
Significance: 0.003
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.307 0.376 0.637 0.643
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Starting gl.keep.pop
Processing genlight object with SNP data
Temporarily assigning pop as population
Checking for presence of nominated populations
Retaining only populations West
Locus metrics not recalculated
Resetting population assignments to initial state
Completed: gl.keep.pop
Starting gl.pcoa
Processing genlight object with SNP data
Performing a PCA, individuals as entities, loci as attributes, SNP genotype as state
Starting gl.colors
Selected color type 2
Completed: gl.colors
Completed: gl.pcoa
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Warning: Removed 36 rows containing non-finite outside the scale range
(`stat_smooth()`).
Warning: Removed 36 rows containing missing values or values outside the scale range
(`geom_point()`).
Removed 36 rows containing missing values or values outside the scale range
(`geom_point()`).
Removed 36 rows containing missing values or values outside the scale range
(`geom_point()`).
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5041
Significance: 0.014
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.282 0.370 0.431 0.520
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Between population PCoA and IBD
# figure 4 ---------------------------------------------------------------------
## pcoa figures ----------------------------------------------------------------
gl.pcoa.plot(pc, vac.gl1, xaxis = 1, yaxis = 2,pt.colors = popcolors)
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
gl.pcoa.plot(pc, vac.gl1, xaxis = 3, yaxis = 4, pt.colors = popcolors)
Starting gl.pcoa.plot
Processing an ordination file (glPca)
Processing genlight object with SNP data
Plotting populations in a space defined by the SNPs
Preparing plot .... please wait
Completed: gl.pcoa.plot
## ibd figures -----------------------------------------------------------------
# pop
gl.ibd(vac.gl17, distance = "propShared",
paircols = "pop", Dgeo_trans = "Dgeo")
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5168
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.0795 0.1093 0.1340 0.1654
Permutation: free
Number of permutations: 999
Completed: gl.ibd
$Dgen
ABTC11930_aca ABTC11932_aca ABTC72881_aca ABTC8929_aca
ABTC11932_aca 0.08454080
ABTC72881_aca 0.09019160 0.09276463
ABTC8929_aca 0.09041666 0.09189282 0.09039321
ABTC8938_aca 0.08500393 0.08827010 0.08764128 0.08227897
AA095616 0.14224510 0.14444736 0.13983353 0.14064384
AA093614 0.13119301 0.13448232 0.13044816 0.13082531
AA095550 0.12975810 0.13164186 0.12685137 0.12859441
AA095620 0.09165329 0.09182624 0.09206501 0.08888659
AA095553 0.31563647 0.31520435 0.31885137 0.32395215
AA095548 0.31082402 0.31180427 0.31435175 0.32007863
AA095617 0.14237363 0.14541258 0.14028226 0.14107857
AA093615 0.13113039 0.13273985 0.13040335 0.13007040
AA095626 0.13033422 0.13312409 0.13003856 0.13026280
AA095636 0.09487763 0.09645876 0.09775272 0.09232150
AA095546 0.30893456 0.30981928 0.31402059 0.31772570
AA096028 0.14419338 0.14585513 0.14119674 0.14266633
AA095562 0.12966304 0.13319552 0.12914342 0.13094168
AA095623 0.13126143 0.13356033 0.12987764 0.13087248
AA095635 0.09249621 0.09352207 0.09418239 0.08801618
AA095547 0.30606513 0.30645016 0.31142805 0.31474297
AA096030 0.14322856 0.14503143 0.14134988 0.14173929
AA095629 0.12968866 0.13258012 0.12829711 0.12909769
AA095624 0.12983838 0.13278626 0.12911250 0.13000746
AA095637 0.09294763 0.09432971 0.09540032 0.08963028
AA095545 0.31091314 0.31031033 0.31445160 0.32056864
AA095568 0.12862718 0.13168002 0.12706821 0.12927880
AA095625 0.13053987 0.13229771 0.12933121 0.13076769
AA095542 0.31102069 0.31275665 0.31571472 0.31905866
AA095570 0.14390257 0.14625238 0.14182635 0.14276606
AA095549 0.13034483 0.13305645 0.12854411 0.12953764
AA095551 0.14404059 0.14530848 0.14083547 0.14210848
AA095543 0.31407635 0.31277976 0.31812137 0.32117153
AA096029 0.14515080 0.14674154 0.14195874 0.14332123
AA095627 0.14235573 0.14473684 0.14041481 0.14077953
AA095540 0.13115636 0.13349027 0.13001218 0.13025591
AA095544 0.31095217 0.31076378 0.31572589 0.31948844
AA095634 0.14348749 0.14623684 0.14175393 0.14237073
AA095569 0.13041644 0.13422845 0.12986408 0.13103968
AA095539a 0.31163857 0.31154875 0.31704579 0.31969287
ABTC8938_aca AA095616 AA093614 AA095550 AA095620
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616 0.14080744
AA093614 0.13065543 0.08678606
AA095550 0.12737336 0.08626468 0.06292768
AA095620 0.08508519 0.14121521 0.12965296 0.12677222
AA095553 0.32049149 0.32998325 0.31788160 0.31575779 0.31824280
AA095548 0.31601644 0.32376396 0.31340504 0.31138502 0.31375667
AA095617 0.14068990 0.04621888 0.08742573 0.08650916 0.14140623
AA093615 0.12977855 0.08841884 0.06305928 0.06167350 0.12855024
AA095626 0.12955431 0.08667520 0.06343646 0.06060468 0.12962919
AA095636 0.08637781 0.14572803 0.13587765 0.13223899 0.08789149
AA095546 0.31362790 0.32429397 0.31249170 0.31110395 0.31095495
AA096028 0.14251904 0.05059806 0.08862362 0.08521521 0.14207077
AA095562 0.12937856 0.08782714 0.05961851 0.06241020 0.12924603
AA095623 0.12932203 0.08756431 0.06422775 0.06061574 0.12932390
AA095635 0.08542981 0.14254812 0.13117772 0.12801478 0.08521251
AA095547 0.31254715 0.32151812 0.30977227 0.30868648 0.31000947
AA096030 0.14189850 0.05053547 0.08806557 0.08551761 0.14125691
AA095629 0.12882123 0.08619119 0.06437394 0.05970354 0.12853692
AA095624 0.12833686 0.08627574 0.06229852 0.05984714 0.12816978
AA095637 0.08614294 0.14745168 0.13509229 0.13170790 0.08641127
AA095545 0.31602387 0.32495526 0.31336418 0.31270592 0.31432967
AA095568 0.12759969 0.08611324 0.06378782 0.06009052 0.12735961
AA095625 0.12938245 0.08696716 0.06325949 0.05972007 0.12873366
AA095542 0.31661065 0.32540439 0.31401151 0.31246329 0.31431394
AA095570 0.14327493 0.05170547 0.08973352 0.08717393 0.14286584
AA095549 0.12892716 0.08774660 0.06116545 0.06216928 0.12857554
AA095551 0.14228480 0.05154941 0.09015148 0.08739275 0.14229508
AA095543 0.31918316 0.32754955 0.31557041 0.31386374 0.31621192
AA096029 0.14301233 0.05329794 0.09180275 0.08879880 0.14311669
AA095627 0.14095951 0.05025390 0.08603755 0.08504056 0.14005179
AA095540 0.12939861 0.08779664 0.06551237 0.06225029 0.12987653
AA095544 0.31579404 0.32453682 0.31230384 0.31102992 0.31394411
AA095634 0.14254177 0.05200826 0.09017185 0.08735654 0.14213807
AA095569 0.12921947 0.08781870 0.06429757 0.06114564 0.12912740
AA095539a 0.31581439 0.32352226 0.31305452 0.31195574 0.31291427
AA095553 AA095548 AA095617 AA093615 AA095626 AA095636
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548 0.13256956
AA095617 0.33089508 0.32464320
AA093615 0.31724608 0.31193351 0.08839889
AA095626 0.31527201 0.31075000 0.08694741 0.03364938
AA095636 0.32506875 0.31822304 0.14628306 0.13415417 0.13457442
AA095546 0.13643491 0.13229662 0.32558840 0.31157869 0.30984602 0.31666481
AA096028 0.32914573 0.32416224 0.04992833 0.08827145 0.08693631 0.14716467
AA095562 0.31608899 0.31171344 0.08830582 0.06097394 0.06085567 0.13443139
AA095623 0.31637475 0.31208733 0.08823392 0.03539683 0.02798696 0.13472416
AA095635 0.32174888 0.31548405 0.14306562 0.12967495 0.12932462 0.08565490
AA095547 0.13787410 0.13184008 0.32255030 0.31029298 0.30892445 0.31549152
AA096030 0.32942168 0.32292085 0.05045667 0.08874684 0.08738178 0.14642960
AA095629 0.31512170 0.31085470 0.08667726 0.05506653 0.03393092 0.13368724
AA095624 0.31654082 0.31211803 0.08679105 0.03530623 0.03363794 0.13384089
AA095637 0.32080777 0.31806665 0.14716351 0.13395858 0.13460763 0.08786402
AA095545 0.13651325 0.13323251 0.32681564 0.31353957 0.31223785 0.31887713
AA095568 0.31649923 0.31119569 0.08663958 0.05769056 0.05673038 0.13268005
AA095625 0.31612087 0.31132771 0.08703557 0.03487633 0.02959312 0.13447549
AA095542 0.13577795 0.13287985 0.32634633 0.31327768 0.31228976 0.31829619
AA095570 0.32902223 0.32287158 0.05127300 0.08920974 0.08714534 0.14842900
AA095549 0.31587783 0.31065472 0.08760149 0.06302146 0.06219534 0.13417333
AA095551 0.32990947 0.32503217 0.05106563 0.09047166 0.08976934 0.14725102
AA095543 0.13425109 0.12989502 0.32823199 0.31459298 0.31365191 0.32087072
AA096029 0.33015877 0.32378540 0.05391313 0.09131175 0.08974954 0.14816340
AA095627 0.32834276 0.32255565 0.04947125 0.08750058 0.08565806 0.14576422
AA095540 0.31593782 0.31145889 0.08761641 0.06193108 0.06149244 0.13461171
AA095544 0.13366468 0.13062818 0.32531973 0.31199330 0.31119519 0.31770948
AA095634 0.32956274 0.32407482 0.05115303 0.08920753 0.08764361 0.14659825
AA095569 0.31697344 0.31241060 0.08835102 0.06159852 0.06075713 0.13437068
AA095539a 0.13431550 0.13100307 0.32520523 0.31298616 0.31202731 0.31858067
AA095546 AA096028 AA095562 AA095623 AA095635 AA095547
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028 0.32455743
AA095562 0.31221677 0.08794511
AA095623 0.31131322 0.08768589 0.06094550
AA095635 0.31559237 0.14293665 0.13059436 0.12974316
AA095547 0.12539362 0.32230058 0.30966075 0.31016163 0.31229070
AA096030 0.32430196 0.04735877 0.08836485 0.08805626 0.14282961 0.32162986
AA095629 0.31084479 0.08632537 0.06115443 0.03355674 0.12888741 0.30909260
AA095624 0.31132376 0.08571299 0.06056771 0.02962980 0.12919027 0.30952003
AA095637 0.31642781 0.14721800 0.13503086 0.13427310 0.08341461 0.31333001
AA095545 0.13144768 0.32619475 0.31298818 0.31334090 0.31617039 0.13084569
AA095568 0.31064090 0.08705635 0.06192790 0.05737427 0.12791055 0.30863806
AA095625 0.31109569 0.08699698 0.06102429 0.02986604 0.12892985 0.30944548
AA095542 0.13114873 0.32483646 0.31323019 0.31356067 0.31613885 0.13276991
AA095570 0.32395479 0.03821203 0.08823189 0.08878645 0.14385914 0.32119241
AA095549 0.31028711 0.08847863 0.06390318 0.06296549 0.12945600 0.30884350
AA095551 0.32557097 0.04906167 0.08891378 0.09039065 0.14326554 0.32177308
AA095543 0.13007471 0.32706687 0.31516148 0.31501365 0.31947099 0.13114116
AA096029 0.32422984 0.04767337 0.09102516 0.09082262 0.14430205 0.32133449
AA095627 0.32341082 0.04464449 0.08719926 0.08698304 0.14114280 0.32044824
AA095540 0.31200280 0.08709663 0.06370568 0.06143676 0.13069307 0.30879758
AA095544 0.13178003 0.32475974 0.31232949 0.31233624 0.31651224 0.13071585
AA095634 0.32484949 0.04598886 0.08874634 0.08904473 0.14308616 0.32237819
AA095569 0.31258888 0.08840510 0.06335043 0.06078681 0.12983439 0.30915980
AA095539a 0.13152827 0.32472480 0.31243759 0.31309883 0.31780010 0.12930911
AA096030 AA095629 AA095624 AA095637 AA095545 AA095568
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629 0.08696658
AA095624 0.08626202 0.03330168
AA095637 0.14685931 0.13334048 0.13359050
AA095545 0.32514975 0.31263063 0.31367701 0.31855084
AA095568 0.08733609 0.05725295 0.05578466 0.13267915 0.31208080
AA095625 0.08687419 0.03357051 0.02841547 0.13435548 0.31244552 0.05396822
AA095542 0.32334546 0.31266621 0.31335517 0.31890594 0.13078863 0.31226555
AA095570 0.04853950 0.08681556 0.08627419 0.14777161 0.32504723 0.08869328
AA095549 0.08738697 0.06172205 0.06112322 0.13324189 0.31205608 0.06039514
AA095551 0.04918395 0.08889146 0.08839442 0.14807163 0.32583802 0.08890859
AA095543 0.32652951 0.31408346 0.31404784 0.31944097 0.12723080 0.31335272
AA096029 0.05121461 0.08876394 0.08902335 0.14922636 0.32657999 0.09038502
AA095627 0.04773444 0.08604834 0.08472028 0.14601258 0.32471395 0.08624569
AA095540 0.08826607 0.06161968 0.06096895 0.13499545 0.31264790 0.06200474
AA095544 0.32449216 0.31129795 0.31193242 0.31684114 0.12879307 0.31012859
AA095634 0.04878956 0.08773148 0.08725959 0.14808463 0.32628736 0.08800230
AA095569 0.08756439 0.06062960 0.06035364 0.13409335 0.31376765 0.06288358
AA095539a 0.32405304 0.31162055 0.31274707 0.31735737 0.13268788 0.31195003
AA095625 AA095542 AA095570 AA095549 AA095551 AA095543
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542 0.31266575
AA095570 0.08740318 0.32539532
AA095549 0.06197565 0.31191635 0.08966095
AA095551 0.08911189 0.32564075 0.04956806 0.09000905
AA095543 0.31400658 0.12441947 0.32805086 0.31362910 0.32733366
AA096029 0.08937366 0.32582974 0.04907738 0.09083063 0.04942303 0.32839328
AA095627 0.08559181 0.32409451 0.04110066 0.08718791 0.04936365 0.32652236
AA095540 0.06140551 0.31208592 0.08758136 0.06438560 0.08926535 0.31476137
AA095544 0.31165319 0.12388355 0.32464949 0.31040625 0.32560650 0.12329120
AA095634 0.08772173 0.32616434 0.04782081 0.08963205 0.05024179 0.32772451
AA095569 0.05967139 0.31391008 0.08853196 0.06395083 0.08950890 0.31518086
AA095539a 0.31255259 0.13192559 0.32418748 0.31151328 0.32534540 0.13071839
AA096029 AA095627 AA095540 AA095544 AA095634 AA095569
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542
AA095570
AA095549
AA095551
AA095543
AA096029
AA095627 0.04966153
AA095540 0.08978857 0.08732479
AA095544 0.32497694 0.32375252 0.31138196
AA095634 0.05006785 0.04752622 0.08864646 0.32527701
AA095569 0.08998039 0.08654804 0.06300766 0.31173884 0.08800186
AA095539a 0.32503093 0.32269768 0.31156178 0.13118967 0.32574257 0.31221582
$Dgeo
ABTC11930_aca ABTC11932_aca ABTC72881_aca ABTC8929_aca
ABTC11932_aca 1.479506e-01
ABTC72881_aca 2.559400e+05 2.559401e+05
ABTC8929_aca 2.647148e+05 2.647146e+05 3.105123e+05
ABTC8938_aca 1.806023e+05 1.806022e+05 1.979688e+05 1.206174e+05
AA095616 7.732494e+05 7.732492e+05 1.008649e+06 7.735525e+05
AA093614 4.289002e+05 4.289002e+05 6.680925e+05 6.395803e+05
AA095550 4.315584e+05 4.315584e+05 6.759576e+05 6.314426e+05
AA095620 1.745281e+05 1.745280e+05 1.842651e+05 1.353572e+05
AA095553 4.648706e+05 4.648707e+05 5.853559e+05 7.294889e+05
AA095548 4.647510e+05 4.647511e+05 5.844383e+05 7.293445e+05
AA095617 7.466274e+05 7.466273e+05 9.820426e+05 7.480617e+05
AA093615 4.288560e+05 4.288560e+05 6.680680e+05 6.395041e+05
AA095626 4.288551e+05 4.288551e+05 6.680683e+05 6.395011e+05
AA095636 1.815843e+05 1.815842e+05 2.155105e+05 1.059797e+05
AA095546 4.648706e+05 4.648708e+05 5.853559e+05 7.294890e+05
AA096028 6.002876e+05 6.002875e+05 8.345819e+05 6.060720e+05
AA095562 4.288563e+05 4.288563e+05 6.680679e+05 6.395051e+05
AA095623 4.288549e+05 4.288549e+05 6.680683e+05 6.395004e+05
AA095635 1.815874e+05 1.815873e+05 2.155151e+05 1.059746e+05
AA095547 4.648673e+05 4.648675e+05 5.853538e+05 7.294857e+05
AA096030 6.002854e+05 6.002853e+05 8.345800e+05 6.060709e+05
AA095629 4.288561e+05 4.288561e+05 6.680680e+05 6.395043e+05
AA095624 4.288558e+05 4.288558e+05 6.680681e+05 6.395033e+05
AA095637 1.745293e+05 1.745293e+05 1.842675e+05 1.353545e+05
AA095545 4.917123e+05 4.917124e+05 5.912690e+05 7.550987e+05
AA095568 4.447169e+05 4.447169e+05 6.839464e+05 6.542060e+05
AA095625 4.288562e+05 4.288563e+05 6.680679e+05 6.395049e+05
AA095542 4.759495e+05 4.759496e+05 5.794111e+05 7.395803e+05
AA095570 6.002926e+05 6.002925e+05 8.345862e+05 6.060745e+05
AA095549 4.447172e+05 4.447172e+05 6.839463e+05 6.542070e+05
AA095551 6.009385e+05 6.009384e+05 8.349173e+05 6.057533e+05
AA095543 4.759425e+05 4.759426e+05 5.793617e+05 7.395693e+05
AA096029 6.003361e+05 6.003360e+05 8.346206e+05 6.060886e+05
AA095627 6.002880e+05 6.002879e+05 8.345822e+05 6.060722e+05
AA095540 4.409964e+05 4.409964e+05 6.824284e+05 6.465696e+05
AA095544 4.757600e+05 4.757601e+05 5.792618e+05 7.393929e+05
AA095634 6.003280e+05 6.003279e+05 8.346130e+05 6.060825e+05
AA095569 4.324983e+05 4.324983e+05 6.766541e+05 6.328738e+05
AA095539a 4.647550e+05 4.647551e+05 5.852496e+05 7.293733e+05
ABTC8938_aca AA095616 AA093614 AA095550 AA095620
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616 8.327245e+05
AA093614 5.982879e+05 6.174851e+05
AA095550 5.963781e+05 5.835782e+05 3.391609e+04
AA095620 1.479670e+04 8.419934e+05 5.960629e+05 5.949874e+05
AA095553 6.334209e+05 9.798642e+05 3.712245e+05 4.037582e+05 6.231797e+05
AA095548 6.330729e+05 9.812641e+05 3.728293e+05 4.053365e+05 6.227998e+05
AA095617 8.064215e+05 2.665915e+04 5.962914e+05 5.623756e+05 8.156321e+05
AA093615 5.982301e+05 6.173752e+05 1.211815e+02 3.380511e+04 5.960078e+05
AA095626 5.982283e+05 6.173684e+05 1.279293e+02 3.379828e+04 5.960062e+05
AA095636 1.759474e+04 8.168964e+05 5.934916e+05 5.906269e+05 3.156877e+04
AA095546 6.334210e+05 9.798643e+05 3.712245e+05 4.037583e+05 6.231798e+05
AA096028 6.600335e+05 1.740869e+05 4.953333e+05 4.621578e+05 6.689795e+05
AA095562 5.982307e+05 6.173776e+05 1.188520e+02 3.380748e+04 5.960083e+05
AA095623 5.982279e+05 6.173668e+05 1.294963e+02 3.379670e+04 5.960058e+05
AA095635 1.759910e+04 8.168936e+05 5.934930e+05 5.906280e+05 3.157358e+04
AA095547 6.334180e+05 9.798605e+05 3.712210e+05 4.037547e+05 6.231768e+05
AA096030 6.600319e+05 1.740889e+05 4.953303e+05 4.621548e+05 6.689779e+05
AA095629 5.982302e+05 6.173757e+05 1.206880e+02 3.380561e+04 5.960079e+05
AA095624 5.982296e+05 6.173735e+05 1.228971e+02 3.380337e+04 5.960074e+05
AA095637 1.479403e+04 8.419919e+05 5.960636e+05 5.949879e+05 2.682530e+00
AA095545 6.544681e+05 1.029991e+06 4.221159e+05 4.546300e+05 6.434768e+05
AA095568 6.137958e+05 6.151320e+05 1.589271e+04 3.540922e+04 6.116504e+05
AA095625 5.982305e+05 6.173771e+05 1.193680e+02 3.380695e+04 5.960082e+05
AA095542 6.396142e+05 1.015144e+06 4.088057e+05 4.410892e+05 6.287432e+05
AA095570 6.600370e+05 1.740826e+05 4.953402e+05 4.621647e+05 6.689832e+05
AA095549 6.137964e+05 6.151343e+05 1.589278e+04 3.541152e+04 6.116509e+05
AA095551 6.600919e+05 1.737321e+05 4.972879e+05 4.641276e+05 6.690772e+05
AA095543 6.395936e+05 1.015212e+06 4.088855e+05 4.411675e+05 6.287210e+05
AA096029 6.600631e+05 1.740471e+05 4.954125e+05 4.622372e+05 6.690104e+05
AA095627 6.600337e+05 1.740866e+05 4.953339e+05 4.621584e+05 6.689798e+05
AA095540 6.084028e+05 6.018546e+05 1.913132e+04 2.159549e+04 6.065700e+05
AA095544 6.394333e+05 1.014979e+06 4.086623e+05 4.409426e+05 6.285635e+05
AA095634 6.600561e+05 1.740548e+05 4.954047e+05 4.622294e+05 6.690034e+05
AA095569 5.975481e+05 5.850930e+05 3.239465e+04 1.908246e+03 5.961199e+05
AA095539a 6.333057e+05 9.797935e+05 3.711706e+05 4.037018e+05 6.230647e+05
AA095553 AA095548 AA095617 AA093615 AA095626
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548 1.916028e+03
AA095617 9.570266e+05 9.584078e+05
AA093615 3.713162e+05 3.729208e+05 5.961802e+05
AA095626 3.713224e+05 3.729270e+05 5.961734e+05 6.937866e+00
AA095636 6.391648e+05 6.388662e+05 7.906359e+05 5.934308e+05 5.934289e+05
AA095546 6.285262e-02 1.915997e+03 9.570267e+05 3.713163e+05 3.713225e+05
AA096028 8.400489e+05 8.412906e+05 1.474645e+05 4.952148e+05 4.952079e+05
AA095562 3.713141e+05 3.729187e+05 5.961826e+05 2.400323e+00 9.338189e+00
AA095623 3.713238e+05 3.729284e+05 5.961718e+05 8.546075e+00 1.608209e+00
AA095635 6.391688e+05 6.388701e+05 7.906332e+05 5.934322e+05 5.934303e+05
AA095547 3.875288e+00 1.917960e+03 9.570229e+05 3.713128e+05 3.713190e+05
AA096030 8.400460e+05 8.412876e+05 1.474664e+05 4.952118e+05 4.952049e+05
AA095629 3.713158e+05 3.729204e+05 5.961807e+05 5.082852e-01 7.446151e+00
AA095624 3.713178e+05 3.729224e+05 5.961785e+05 1.765922e+00 5.171943e+00
AA095637 6.231817e+05 6.228018e+05 8.156306e+05 5.960084e+05 5.960068e+05
AA095545 5.089205e+04 4.929430e+04 1.007028e+06 4.222076e+05 4.222138e+05
AA095568 3.782367e+05 3.798801e+05 5.943659e+05 1.592592e+04 1.592619e+04
AA095625 3.713146e+05 3.729191e+05 5.961821e+05 1.868420e+00 8.806286e+00
AA095542 3.881126e+04 3.699378e+04 9.920808e+05 4.088955e+05 4.089016e+05
AA095570 8.400557e+05 8.412973e+05 1.474601e+05 4.952217e+05 4.952148e+05
AA095549 3.782347e+05 3.798780e+05 5.943683e+05 1.592601e+04 1.592629e+04
AA095551 8.417367e+05 8.429759e+05 1.471312e+05 4.971693e+05 4.971624e+05
AA095543 3.890289e+04 3.708478e+04 9.921476e+05 4.089752e+05 4.089813e+05
AA096029 8.401240e+05 8.413656e+05 1.474252e+05 4.952940e+05 4.952871e+05
AA095627 8.400495e+05 8.412911e+05 1.474641e+05 4.952154e+05 4.952084e+05
AA095540 3.893781e+05 3.909992e+05 5.809683e+05 1.907035e+04 1.906554e+04
AA095544 3.870400e+04 3.688393e+04 9.919143e+05 4.087520e+05 4.087582e+05
AA095634 8.401155e+05 8.413571e+05 1.474328e+05 4.952861e+05 4.952792e+05
AA095569 4.025826e+05 4.041650e+05 5.639225e+05 3.228544e+04 3.227865e+04
AA095539a 1.156729e+02 1.908561e+03 9.569547e+05 3.712623e+05 3.712685e+05
AA095636 AA095546 AA096028 AA095562 AA095623
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546 6.391649e+05
AA096028 6.444533e+05 8.400490e+05
AA095562 5.934315e+05 3.713141e+05 4.952172e+05
AA095623 5.934284e+05 3.713239e+05 4.952063e+05 1.094640e+01
AA095635 5.104533e+00 6.391688e+05 6.444507e+05 5.934329e+05 5.934298e+05
AA095547 6.391618e+05 3.938140e+00 8.400451e+05 3.713106e+05 3.713204e+05
AA096030 6.444518e+05 8.400460e+05 3.014050e+00 4.952142e+05 4.952033e+05
AA095629 5.934309e+05 3.713158e+05 4.952153e+05 1.892038e+00 9.054360e+00
AA095624 5.934303e+05 3.713179e+05 4.952131e+05 4.166246e+00 6.780152e+00
AA095637 3.156624e+04 6.231818e+05 6.689781e+05 5.960089e+05 5.960064e+05
AA095545 6.613808e+05 5.089200e+04 8.889012e+05 4.222055e+05 4.222152e+05
AA095568 6.088960e+05 3.782368e+05 4.965936e+05 1.592582e+04 1.592626e+04
AA095625 5.934313e+05 3.713146e+05 4.952167e+05 5.319035e-01 1.041449e+01
AA095542 6.463415e+05 3.881121e+04 8.733243e+05 4.088934e+05 4.089030e+05
AA095570 6.444568e+05 8.400557e+05 6.861102e+00 4.952241e+05 4.952132e+05
AA095549 6.088967e+05 3.782347e+05 4.965960e+05 1.592592e+04 1.592636e+04
AA095551 6.444840e+05 8.417368e+05 2.216430e+03 4.971717e+05 4.971608e+05
AA095543 6.463233e+05 3.890284e+04 8.733838e+05 4.089731e+05 4.089828e+05
AA096029 6.444820e+05 8.401240e+05 7.992089e+01 4.952964e+05 4.952855e+05
AA095627 6.444536e+05 8.400495e+05 5.364397e-01 4.952178e+05 4.952068e+05
AA095540 6.031469e+05 3.893782e+05 4.826337e+05 1.907201e+04 1.906443e+04
AA095544 6.461587e+05 3.870395e+04 8.731489e+05 4.087499e+05 4.087596e+05
AA095634 6.444750e+05 8.401156e+05 7.277465e+01 4.952885e+05 4.952776e+05
AA095569 5.918387e+05 4.025827e+05 4.639055e+05 3.228779e+04 3.227708e+04
AA095539a 6.390494e+05 1.157242e+02 8.399682e+05 3.712602e+05 3.712699e+05
AA095635 AA095547 AA096030 AA095629 AA095624
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547 6.391658e+05
AA096030 6.444491e+05 8.400422e+05
AA095629 5.934324e+05 3.713123e+05 4.952123e+05
AA095624 5.934317e+05 3.713143e+05 4.952101e+05 2.274208e+00
AA095637 3.157105e+04 6.231788e+05 6.689765e+05 5.960085e+05 5.960080e+05
AA095545 6.613850e+05 5.089550e+04 8.888983e+05 4.222072e+05 4.222092e+05
AA095568 6.088974e+05 3.782334e+05 4.965906e+05 1.592589e+04 1.592599e+04
AA095625 5.934328e+05 3.713111e+05 4.952137e+05 1.360135e+00 3.634342e+00
AA095542 6.463456e+05 3.881413e+04 8.733213e+05 4.088950e+05 4.088970e+05
AA095570 6.444542e+05 8.400519e+05 9.875152e+00 4.952222e+05 4.952199e+05
AA095549 6.088981e+05 3.782313e+05 4.965930e+05 1.592599e+04 1.592608e+04
AA095551 6.444814e+05 8.417329e+05 2.219042e+03 4.971699e+05 4.971676e+05
AA095543 6.463274e+05 3.890576e+04 8.733809e+05 4.089748e+05 4.089768e+05
AA096029 6.444794e+05 8.401202e+05 8.289407e+01 4.952945e+05 4.952922e+05
AA095627 6.444510e+05 8.400456e+05 3.550489e+00 4.952159e+05 4.952136e+05
AA095540 6.031482e+05 3.893747e+05 4.826307e+05 1.907070e+04 1.906913e+04
AA095544 6.461629e+05 3.870686e+04 8.731459e+05 4.087516e+05 4.087536e+05
AA095634 6.444724e+05 8.401117e+05 7.571145e+01 4.952867e+05 4.952844e+05
AA095569 5.918399e+05 4.025791e+05 4.639025e+05 3.228594e+04 3.228371e+04
AA095539a 6.390533e+05 1.125312e+02 8.399653e+05 3.712619e+05 3.712639e+05
AA095637 AA095545 AA095568 AA095625 AA095542
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545 6.434790e+05
AA095568 6.116510e+05 4.291045e+05
AA095625 5.960088e+05 4.222060e+05 1.592584e+04
AA095542 6.287453e+05 1.669474e+04 4.161431e+05 4.088938e+05
AA095570 6.689818e+05 8.889080e+05 4.966005e+05 4.952236e+05 8.733310e+05
AA095549 6.116516e+05 4.291024e+05 2.401226e+00 1.592594e+04 4.161411e+05
AA095551 6.690757e+05 8.905650e+05 4.985805e+05 4.971712e+05 8.749781e+05
AA095543 6.287231e+05 1.666843e+04 4.162249e+05 4.089736e+05 9.944637e+01
AA096029 6.690090e+05 8.889758e+05 4.966731e+05 4.952959e+05 8.733987e+05
AA095627 6.689784e+05 8.889018e+05 4.965942e+05 4.952172e+05 8.733248e+05
AA095540 6.065705e+05 4.402676e+05 1.402144e+04 1.907165e+04 4.270970e+05
AA095544 6.285656e+05 1.688971e+04 4.160044e+05 4.087504e+05 1.957239e+02
AA095634 6.690019e+05 8.889673e+05 4.966654e+05 4.952880e+05 8.733902e+05
AA095569 5.961204e+05 4.534591e+05 3.353010e+04 3.228727e+04 4.399516e+05
AA095539a 6.230667e+05 5.094551e+04 3.781868e+05 3.712607e+05 3.883660e+04
AA095570 AA095549 AA095551 AA095543 AA096029
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542
AA095570
AA095549 4.966029e+05
AA095551 2.210488e+03 4.985829e+05
AA095543 8.733906e+05 4.162228e+05 8.750375e+05
AA096029 7.316530e+01 4.966755e+05 2.141666e+03 8.734582e+05
AA095627 6.324663e+00 4.965966e+05 2.215965e+03 8.733844e+05 7.939205e+01
AA095540 4.826406e+05 1.402381e+04 4.846180e+05 4.271775e+05 4.827132e+05
AA095544 8.731556e+05 4.160023e+05 8.748025e+05 2.350750e+02 8.732232e+05
AA095634 6.611534e+01 4.966678e+05 2.146821e+03 8.734497e+05 8.635395e+00
AA095569 4.639124e+05 3.353241e+04 4.658770e+05 4.400301e+05 4.639849e+05
AA095539a 8.399750e+05 3.781847e+05 8.416559e+05 3.892812e+04 8.400433e+05
AA095627 AA095540 AA095544 AA095634 AA095569
ABTC11932_aca
ABTC72881_aca
ABTC8929_aca
ABTC8938_aca
AA095616
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542
AA095570
AA095549
AA095551
AA095543
AA096029
AA095627
AA095540 4.826343e+05
AA095544 8.731494e+05 4.269554e+05
AA095634 7.225263e+01 4.827054e+05 8.732148e+05
AA095569 4.639061e+05 1.969448e+04 4.398054e+05 4.639771e+05
AA095539a 8.399688e+05 3.893259e+05 3.872887e+04 8.400348e+05 4.025266e+05
$mantel
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5168
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.0795 0.1093 0.1340 0.1654
Permutation: free
Number of permutations: 999
# karyo
gl.ibd(vac.gl18, distance = "propShared",
paircols = "pop", Dgeo_trans = "Dgeo")
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5604
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.090 0.119 0.154 0.183
Permutation: free
Number of permutations: 999
Completed: gl.ibd
$Dgen
AA095616 AA093614 AA095550 AA095620 AA095553 AA095548
AA093614 0.08678606
AA095550 0.08626468 0.06292768
AA095620 0.14121521 0.12965296 0.12677222
AA095553 0.32998325 0.31788160 0.31575779 0.31824280
AA095548 0.32376396 0.31340504 0.31138502 0.31375667 0.13256956
AA095617 0.04621888 0.08742573 0.08650916 0.14140623 0.33089508 0.32464320
AA093615 0.08841884 0.06305928 0.06167350 0.12855024 0.31724608 0.31193351
AA095626 0.08667520 0.06343646 0.06060468 0.12962919 0.31527201 0.31075000
AA095636 0.14572803 0.13587765 0.13223899 0.08789149 0.32506875 0.31822304
AA095546 0.32429397 0.31249170 0.31110395 0.31095495 0.13643491 0.13229662
AA096028 0.05059806 0.08862362 0.08521521 0.14207077 0.32914573 0.32416224
AA095562 0.08782714 0.05961851 0.06241020 0.12924603 0.31608899 0.31171344
AA095623 0.08756431 0.06422775 0.06061574 0.12932390 0.31637475 0.31208733
AA095635 0.14254812 0.13117772 0.12801478 0.08521251 0.32174888 0.31548405
AA095547 0.32151812 0.30977227 0.30868648 0.31000947 0.13787410 0.13184008
AA096030 0.05053547 0.08806557 0.08551761 0.14125691 0.32942168 0.32292085
AA095629 0.08619119 0.06437394 0.05970354 0.12853692 0.31512170 0.31085470
AA095624 0.08627574 0.06229852 0.05984714 0.12816978 0.31654082 0.31211803
AA095637 0.14745168 0.13509229 0.13170790 0.08641127 0.32080777 0.31806665
AA095545 0.32495526 0.31336418 0.31270592 0.31432967 0.13651325 0.13323251
AA095568 0.08611324 0.06378782 0.06009052 0.12735961 0.31649923 0.31119569
AA095625 0.08696716 0.06325949 0.05972007 0.12873366 0.31612087 0.31132771
AA095542 0.32540439 0.31401151 0.31246329 0.31431394 0.13577795 0.13287985
AA095570 0.05170547 0.08973352 0.08717393 0.14286584 0.32902223 0.32287158
AA095543 0.32754955 0.31557041 0.31386374 0.31621192 0.13425109 0.12989502
AA096029 0.05329794 0.09180275 0.08879880 0.14311669 0.33015877 0.32378540
AA095627 0.05025390 0.08603755 0.08504056 0.14005179 0.32834276 0.32255565
AA095544 0.32453682 0.31230384 0.31102992 0.31394411 0.13366468 0.13062818
AA095634 0.05200826 0.09017185 0.08735654 0.14213807 0.32956274 0.32407482
AA095569 0.08781870 0.06429757 0.06114564 0.12912740 0.31697344 0.31241060
AA095617 AA093615 AA095626 AA095636 AA095546 AA096028
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615 0.08839889
AA095626 0.08694741 0.03364938
AA095636 0.14628306 0.13415417 0.13457442
AA095546 0.32558840 0.31157869 0.30984602 0.31666481
AA096028 0.04992833 0.08827145 0.08693631 0.14716467 0.32455743
AA095562 0.08830582 0.06097394 0.06085567 0.13443139 0.31221677 0.08794511
AA095623 0.08823392 0.03539683 0.02798696 0.13472416 0.31131322 0.08768589
AA095635 0.14306562 0.12967495 0.12932462 0.08565490 0.31559237 0.14293665
AA095547 0.32255030 0.31029298 0.30892445 0.31549152 0.12539362 0.32230058
AA096030 0.05045667 0.08874684 0.08738178 0.14642960 0.32430196 0.04735877
AA095629 0.08667726 0.05506653 0.03393092 0.13368724 0.31084479 0.08632537
AA095624 0.08679105 0.03530623 0.03363794 0.13384089 0.31132376 0.08571299
AA095637 0.14716351 0.13395858 0.13460763 0.08786402 0.31642781 0.14721800
AA095545 0.32681564 0.31353957 0.31223785 0.31887713 0.13144768 0.32619475
AA095568 0.08663958 0.05769056 0.05673038 0.13268005 0.31064090 0.08705635
AA095625 0.08703557 0.03487633 0.02959312 0.13447549 0.31109569 0.08699698
AA095542 0.32634633 0.31327768 0.31228976 0.31829619 0.13114873 0.32483646
AA095570 0.05127300 0.08920974 0.08714534 0.14842900 0.32395479 0.03821203
AA095543 0.32823199 0.31459298 0.31365191 0.32087072 0.13007471 0.32706687
AA096029 0.05391313 0.09131175 0.08974954 0.14816340 0.32422984 0.04767337
AA095627 0.04947125 0.08750058 0.08565806 0.14576422 0.32341082 0.04464449
AA095544 0.32531973 0.31199330 0.31119519 0.31770948 0.13178003 0.32475974
AA095634 0.05115303 0.08920753 0.08764361 0.14659825 0.32484949 0.04598886
AA095569 0.08835102 0.06159852 0.06075713 0.13437068 0.31258888 0.08840510
AA095562 AA095623 AA095635 AA095547 AA096030 AA095629
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623 0.06094550
AA095635 0.13059436 0.12974316
AA095547 0.30966075 0.31016163 0.31229070
AA096030 0.08836485 0.08805626 0.14282961 0.32162986
AA095629 0.06115443 0.03355674 0.12888741 0.30909260 0.08696658
AA095624 0.06056771 0.02962980 0.12919027 0.30952003 0.08626202 0.03330168
AA095637 0.13503086 0.13427310 0.08341461 0.31333001 0.14685931 0.13334048
AA095545 0.31298818 0.31334090 0.31617039 0.13084569 0.32514975 0.31263063
AA095568 0.06192790 0.05737427 0.12791055 0.30863806 0.08733609 0.05725295
AA095625 0.06102429 0.02986604 0.12892985 0.30944548 0.08687419 0.03357051
AA095542 0.31323019 0.31356067 0.31613885 0.13276991 0.32334546 0.31266621
AA095570 0.08823189 0.08878645 0.14385914 0.32119241 0.04853950 0.08681556
AA095543 0.31516148 0.31501365 0.31947099 0.13114116 0.32652951 0.31408346
AA096029 0.09102516 0.09082262 0.14430205 0.32133449 0.05121461 0.08876394
AA095627 0.08719926 0.08698304 0.14114280 0.32044824 0.04773444 0.08604834
AA095544 0.31232949 0.31233624 0.31651224 0.13071585 0.32449216 0.31129795
AA095634 0.08874634 0.08904473 0.14308616 0.32237819 0.04878956 0.08773148
AA095569 0.06335043 0.06078681 0.12983439 0.30915980 0.08756439 0.06062960
AA095624 AA095637 AA095545 AA095568 AA095625 AA095542
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637 0.13359050
AA095545 0.31367701 0.31855084
AA095568 0.05578466 0.13267915 0.31208080
AA095625 0.02841547 0.13435548 0.31244552 0.05396822
AA095542 0.31335517 0.31890594 0.13078863 0.31226555 0.31266575
AA095570 0.08627419 0.14777161 0.32504723 0.08869328 0.08740318 0.32539532
AA095543 0.31404784 0.31944097 0.12723080 0.31335272 0.31400658 0.12441947
AA096029 0.08902335 0.14922636 0.32657999 0.09038502 0.08937366 0.32582974
AA095627 0.08472028 0.14601258 0.32471395 0.08624569 0.08559181 0.32409451
AA095544 0.31193242 0.31684114 0.12879307 0.31012859 0.31165319 0.12388355
AA095634 0.08725959 0.14808463 0.32628736 0.08800230 0.08772173 0.32616434
AA095569 0.06035364 0.13409335 0.31376765 0.06288358 0.05967139 0.31391008
AA095570 AA095543 AA096029 AA095627 AA095544 AA095634
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542
AA095570
AA095543 0.32805086
AA096029 0.04907738 0.32839328
AA095627 0.04110066 0.32652236 0.04966153
AA095544 0.32464949 0.12329120 0.32497694 0.32375252
AA095634 0.04782081 0.32772451 0.05006785 0.04752622 0.32527701
AA095569 0.08853196 0.31518086 0.08998039 0.08654804 0.31173884 0.08800186
$Dgeo
AA095616 AA093614 AA095550 AA095620 AA095553
AA093614 6.174851e+05
AA095550 5.835782e+05 3.391609e+04
AA095620 8.419934e+05 5.960629e+05 5.949874e+05
AA095553 9.798642e+05 3.712245e+05 4.037582e+05 6.231797e+05
AA095548 9.812641e+05 3.728293e+05 4.053365e+05 6.227998e+05 1.916028e+03
AA095617 2.665915e+04 5.962914e+05 5.623756e+05 8.156321e+05 9.570266e+05
AA093615 6.173752e+05 1.211815e+02 3.380511e+04 5.960078e+05 3.713162e+05
AA095626 6.173684e+05 1.279293e+02 3.379828e+04 5.960062e+05 3.713224e+05
AA095636 8.168964e+05 5.934916e+05 5.906269e+05 3.156877e+04 6.391648e+05
AA095546 9.798643e+05 3.712245e+05 4.037583e+05 6.231798e+05 6.285262e-02
AA096028 1.740869e+05 4.953333e+05 4.621578e+05 6.689795e+05 8.400489e+05
AA095562 6.173776e+05 1.188520e+02 3.380748e+04 5.960083e+05 3.713141e+05
AA095623 6.173668e+05 1.294963e+02 3.379670e+04 5.960058e+05 3.713238e+05
AA095635 8.168936e+05 5.934930e+05 5.906280e+05 3.157358e+04 6.391688e+05
AA095547 9.798605e+05 3.712210e+05 4.037547e+05 6.231768e+05 3.875288e+00
AA096030 1.740889e+05 4.953303e+05 4.621548e+05 6.689779e+05 8.400460e+05
AA095629 6.173757e+05 1.206880e+02 3.380561e+04 5.960079e+05 3.713158e+05
AA095624 6.173735e+05 1.228971e+02 3.380337e+04 5.960074e+05 3.713178e+05
AA095637 8.419919e+05 5.960636e+05 5.949879e+05 2.682530e+00 6.231817e+05
AA095545 1.029991e+06 4.221159e+05 4.546300e+05 6.434768e+05 5.089205e+04
AA095568 6.151320e+05 1.589271e+04 3.540922e+04 6.116504e+05 3.782367e+05
AA095625 6.173771e+05 1.193680e+02 3.380695e+04 5.960082e+05 3.713146e+05
AA095542 1.015144e+06 4.088057e+05 4.410892e+05 6.287432e+05 3.881126e+04
AA095570 1.740826e+05 4.953402e+05 4.621647e+05 6.689832e+05 8.400557e+05
AA095543 1.015212e+06 4.088855e+05 4.411675e+05 6.287210e+05 3.890289e+04
AA096029 1.740471e+05 4.954125e+05 4.622372e+05 6.690104e+05 8.401240e+05
AA095627 1.740866e+05 4.953339e+05 4.621584e+05 6.689798e+05 8.400495e+05
AA095544 1.014979e+06 4.086623e+05 4.409426e+05 6.285635e+05 3.870400e+04
AA095634 1.740548e+05 4.954047e+05 4.622294e+05 6.690034e+05 8.401155e+05
AA095569 5.850930e+05 3.239465e+04 1.908246e+03 5.961199e+05 4.025826e+05
AA095548 AA095617 AA093615 AA095626 AA095636
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617 9.584078e+05
AA093615 3.729208e+05 5.961802e+05
AA095626 3.729270e+05 5.961734e+05 6.937866e+00
AA095636 6.388662e+05 7.906359e+05 5.934308e+05 5.934289e+05
AA095546 1.915997e+03 9.570267e+05 3.713163e+05 3.713225e+05 6.391649e+05
AA096028 8.412906e+05 1.474645e+05 4.952148e+05 4.952079e+05 6.444533e+05
AA095562 3.729187e+05 5.961826e+05 2.400323e+00 9.338189e+00 5.934315e+05
AA095623 3.729284e+05 5.961718e+05 8.546075e+00 1.608209e+00 5.934284e+05
AA095635 6.388701e+05 7.906332e+05 5.934322e+05 5.934303e+05 5.104533e+00
AA095547 1.917960e+03 9.570229e+05 3.713128e+05 3.713190e+05 6.391618e+05
AA096030 8.412876e+05 1.474664e+05 4.952118e+05 4.952049e+05 6.444518e+05
AA095629 3.729204e+05 5.961807e+05 5.082852e-01 7.446151e+00 5.934309e+05
AA095624 3.729224e+05 5.961785e+05 1.765922e+00 5.171943e+00 5.934303e+05
AA095637 6.228018e+05 8.156306e+05 5.960084e+05 5.960068e+05 3.156624e+04
AA095545 4.929430e+04 1.007028e+06 4.222076e+05 4.222138e+05 6.613808e+05
AA095568 3.798801e+05 5.943659e+05 1.592592e+04 1.592619e+04 6.088960e+05
AA095625 3.729191e+05 5.961821e+05 1.868420e+00 8.806286e+00 5.934313e+05
AA095542 3.699378e+04 9.920808e+05 4.088955e+05 4.089016e+05 6.463415e+05
AA095570 8.412973e+05 1.474601e+05 4.952217e+05 4.952148e+05 6.444568e+05
AA095543 3.708478e+04 9.921476e+05 4.089752e+05 4.089813e+05 6.463233e+05
AA096029 8.413656e+05 1.474252e+05 4.952940e+05 4.952871e+05 6.444820e+05
AA095627 8.412911e+05 1.474641e+05 4.952154e+05 4.952084e+05 6.444536e+05
AA095544 3.688393e+04 9.919143e+05 4.087520e+05 4.087582e+05 6.461587e+05
AA095634 8.413571e+05 1.474328e+05 4.952861e+05 4.952792e+05 6.444750e+05
AA095569 4.041650e+05 5.639225e+05 3.228544e+04 3.227865e+04 5.918387e+05
AA095546 AA096028 AA095562 AA095623 AA095635
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028 8.400490e+05
AA095562 3.713141e+05 4.952172e+05
AA095623 3.713239e+05 4.952063e+05 1.094640e+01
AA095635 6.391688e+05 6.444507e+05 5.934329e+05 5.934298e+05
AA095547 3.938140e+00 8.400451e+05 3.713106e+05 3.713204e+05 6.391658e+05
AA096030 8.400460e+05 3.014050e+00 4.952142e+05 4.952033e+05 6.444491e+05
AA095629 3.713158e+05 4.952153e+05 1.892038e+00 9.054360e+00 5.934324e+05
AA095624 3.713179e+05 4.952131e+05 4.166246e+00 6.780152e+00 5.934317e+05
AA095637 6.231818e+05 6.689781e+05 5.960089e+05 5.960064e+05 3.157105e+04
AA095545 5.089200e+04 8.889012e+05 4.222055e+05 4.222152e+05 6.613850e+05
AA095568 3.782368e+05 4.965936e+05 1.592582e+04 1.592626e+04 6.088974e+05
AA095625 3.713146e+05 4.952167e+05 5.319035e-01 1.041449e+01 5.934328e+05
AA095542 3.881121e+04 8.733243e+05 4.088934e+05 4.089030e+05 6.463456e+05
AA095570 8.400557e+05 6.861102e+00 4.952241e+05 4.952132e+05 6.444542e+05
AA095543 3.890284e+04 8.733838e+05 4.089731e+05 4.089828e+05 6.463274e+05
AA096029 8.401240e+05 7.992089e+01 4.952964e+05 4.952855e+05 6.444794e+05
AA095627 8.400495e+05 5.364397e-01 4.952178e+05 4.952068e+05 6.444510e+05
AA095544 3.870395e+04 8.731489e+05 4.087499e+05 4.087596e+05 6.461629e+05
AA095634 8.401156e+05 7.277465e+01 4.952885e+05 4.952776e+05 6.444724e+05
AA095569 4.025827e+05 4.639055e+05 3.228779e+04 3.227708e+04 5.918399e+05
AA095547 AA096030 AA095629 AA095624 AA095637
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030 8.400422e+05
AA095629 3.713123e+05 4.952123e+05
AA095624 3.713143e+05 4.952101e+05 2.274208e+00
AA095637 6.231788e+05 6.689765e+05 5.960085e+05 5.960080e+05
AA095545 5.089550e+04 8.888983e+05 4.222072e+05 4.222092e+05 6.434790e+05
AA095568 3.782334e+05 4.965906e+05 1.592589e+04 1.592599e+04 6.116510e+05
AA095625 3.713111e+05 4.952137e+05 1.360135e+00 3.634342e+00 5.960088e+05
AA095542 3.881413e+04 8.733213e+05 4.088950e+05 4.088970e+05 6.287453e+05
AA095570 8.400519e+05 9.875152e+00 4.952222e+05 4.952199e+05 6.689818e+05
AA095543 3.890576e+04 8.733809e+05 4.089748e+05 4.089768e+05 6.287231e+05
AA096029 8.401202e+05 8.289407e+01 4.952945e+05 4.952922e+05 6.690090e+05
AA095627 8.400456e+05 3.550489e+00 4.952159e+05 4.952136e+05 6.689784e+05
AA095544 3.870686e+04 8.731459e+05 4.087516e+05 4.087536e+05 6.285656e+05
AA095634 8.401117e+05 7.571145e+01 4.952867e+05 4.952844e+05 6.690019e+05
AA095569 4.025791e+05 4.639025e+05 3.228594e+04 3.228371e+04 5.961204e+05
AA095545 AA095568 AA095625 AA095542 AA095570
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568 4.291045e+05
AA095625 4.222060e+05 1.592584e+04
AA095542 1.669474e+04 4.161431e+05 4.088938e+05
AA095570 8.889080e+05 4.966005e+05 4.952236e+05 8.733310e+05
AA095543 1.666843e+04 4.162249e+05 4.089736e+05 9.944637e+01 8.733906e+05
AA096029 8.889758e+05 4.966731e+05 4.952959e+05 8.733987e+05 7.316530e+01
AA095627 8.889018e+05 4.965942e+05 4.952172e+05 8.733248e+05 6.324663e+00
AA095544 1.688971e+04 4.160044e+05 4.087504e+05 1.957239e+02 8.731556e+05
AA095634 8.889673e+05 4.966654e+05 4.952880e+05 8.733902e+05 6.611534e+01
AA095569 4.534591e+05 3.353010e+04 3.228727e+04 4.399516e+05 4.639124e+05
AA095543 AA096029 AA095627 AA095544 AA095634
AA093614
AA095550
AA095620
AA095553
AA095548
AA095617
AA093615
AA095626
AA095636
AA095546
AA096028
AA095562
AA095623
AA095635
AA095547
AA096030
AA095629
AA095624
AA095637
AA095545
AA095568
AA095625
AA095542
AA095570
AA095543
AA096029 8.734582e+05
AA095627 8.733844e+05 7.939205e+01
AA095544 2.350750e+02 8.732232e+05 8.731494e+05
AA095634 8.734497e+05 8.635395e+00 7.225263e+01 8.732148e+05
AA095569 4.400301e+05 4.639849e+05 4.639061e+05 4.398054e+05 4.639771e+05
$mantel
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5604
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.090 0.119 0.154 0.183
Permutation: free
Number of permutations: 999
Multiple linear regression model for within and between populations testing same vs different karyotypes
## reviewer figures -------------------------------------------------------------
# reviewer ibd karyo
<- function(gl, analysis = "between"){
review.analysis
if(analysis == "within"){ ndist <- 51; diffCol = "#FFC33B"}
if(analysis == "between"){ ndist <- 1050; diffCol = "#FF5AAF"}
<- gl.ibd(gl, distance = "propShared",
karyoibd paircols = "pop", Dgeo_trans = "Dgeo")
<- karyoibd$Dgeo/1000
mat <- karyoibd$Dgen
prop
<- matrixConvert(mat, c("id1", "id2", "geo"))
dfDist <- matrixConvert(prop, c("id1", "id2", "gen"))
dfGen <- gl
jason2 <- jason2@other$ind.metrics[,c("id", "karyo")]
dfkaryo <- jason2@other$ind.metrics[,c("id", "pop")]
dfpop
<- left_join(dfGen, dfDist) %>%
df left_join(dfkaryo, by = c("id1" = "id")) %>%
left_join(dfkaryo, by = c("id2" = "id")) %>%
left_join(dfpop, by = c("id1" = "id")) %>%
left_join(dfpop, by = c("id2" = "id")) %>%
mutate(pair = paste(karyo.x,"-", karyo.y)) %>%
filter(geo <= ndist) %>% # change based on question between or within
rowwise() %>%
mutate(karyo = ifelse(karyo.x == karyo.y, "same", "different"),
diffpop = pop.x == pop.y)
<- c(df$id1, df$id2)
indexPost
<- lm(gen ~ geo, data = df)
m summary(m)
<- lm(gen ~ geo + karyo, data = df)
m2 summary(m2)
# tab_model(m, m2, show.aic = TRUE, file = "temp.html")
# webshot("temp.html",file = paste0("model_karyo_ibd_", ndist, ".png") ,
# vwidth = 441, vheight = 351)
<- AIC(m, m2)
aic $delta <- aic$AIC - min(aic$AIC)
aic
<- 1:ndist
x
<- predict(m, newdata = data.frame(geo = x), interval = 'confidence')
y1 <- predict(m2, newdata = data.frame(geo = x, karyo = "same"))
y2 <- predict(m2, newdata = data.frame(geo = x, karyo = "different"))
y3
<- data.frame(geo = rep(x,2), gen = c(y2, y3),
predmodel karyo = rep(c("same", "different"), each = ndist))
<- summary(m2)
sm
<- round(m2$coefficients[1],3)
b0 <- signif(m2$coefficients[2], 1)
b1 <- round(m2$coefficients[3], 3)
b2
<- round(sm$r.squared,2)
r2 <- round(sm$coefficients[3,3], 2)
tb2
<- paste0(b0 ,
mExpress " + ", b1, " + ", b2)
<- round(sm$coefficients[3,4],2)
pval
<-ggplot(df, aes(geo, gen, fill = karyo)) +
reviewerPlot geom_point(alpha = 0.7, size = 4, pch = 21) +
theme_classic() +
geom_line(data = filter(predmodel, geo < ndist),
aes(x = geo, y = gen, colour = karyo), linewidth = 1, lty = 5)+
xlab("Distance (km)")+
ylab("1-(Proportion of Shared Alleles)")+
scale_fill_manual(values = c(different = diffCol, same = "black"))+
scale_colour_manual(values = c(different = diffCol, same = "black"))+
# theme(legend.position = "none") +
geom_text(aes(x = 39, y = 0.05, label = c("")), cex = 5)+
theme(legend.position = "bottom")+
labs(fill = "Karyotype", colour = "Karyotype")
return(reviewerPlot)
}
review.analysis(gl = vac.gl18, analysis = "between")
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5604
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.0922 0.1273 0.1537 0.1886
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Joining with `by = join_by(id1, id2)`
review.analysis(gl = vac.gl18, analysis = "within")
Starting gl.ibd
Analysis performed on the genlight object.
Coordinates transformed to Mercator (google) projection to calculate distances in meters.
Coordinates used from: x@other$latlon (Mercator transformed)
Transformation of Dgeo: Dgeo
Genetic distance: propShared
Tranformation of Dgen: Dgen
Mantel statistic based on Pearson's product-moment correlation
Call:
vegan::mantel(xdis = Dgen, ydis = Dgeo, permutations = permutations, na.rm = TRUE)
Mantel statistic r: 0.5604
Significance: 0.001
Upper quantiles of permutations (null model):
90% 95% 97.5% 99%
0.0939 0.1295 0.1625 0.2131
Permutation: free
Number of permutations: 999
Completed: gl.ibd
Joining with `by = join_by(id1, id2)`
Characterising a large chromosomal inversion
<- gl.load("./data/varanus_MBE_Karyotypes.rds") vac.gl
Starting gl.load
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Loaded object of type SNP from ./data/varanus_MBE_Karyotypes.rds
Completed: gl.load
class(vac.gl) <- "dartR"
####Add position for loci metadata
@position <- as.integer(runif(nLoc(vac.gl), 2, 49))
vac.gl
# subset SNPs that align to inferred chromosome pool scaffolds -----------------
# how many reads align to flow sorted chromosome pool?
table(vac.gl@other$loc.metrics$Chrom_Varanus_JD_vac_aligned_chr6_7 != "")
FALSE TRUE
283593 17974
# Index those loci
<-
index6 @other$loc.metrics$Chrom_Varanus_JD_vac_aligned_chr6_7 != ""
vac.gl
# assign alleles aligned to chromosome 6_7 scaffolds to a new genlight object --
<- vac.gl[, index6]
vac.gl6
gl.smearplot(vac.gl6, group.pop = FALSE)
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Starting gl.smearplot
Completed: gl.smearplot
# add DAPI Karyotype individual metric
# report private alleles and fixed alleles per pair of populations for vac.gl6
gl.report.pa(vac.gl6)
Starting gl.report.pa
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
p1 p2 pop1 pop2 N1 N2 fixed priv1 priv2 Chao1 Chao2 totalpriv AFD asym
1 1 2 AA MA 5 7 36 298 1213 72 121 1511 0.037 NA
2 1 3 AA MM 5 22 97 274 3973 72 697 4247 0.078 NA
3 2 3 MA MM 7 22 35 310 3346 121 697 3656 0.063 NA
asym.sig
1 NA
2 NA
3 NA
Table of private alleles and fixed differences returned
Completed: gl.report.pa
p1 p2 pop1 pop2 N1 N2 fixed priv1 priv2 Chao1 Chao2 totalpriv AFD asym
1 1 2 AA MA 5 7 36 298 1213 72 121 1511 0.037 NA
2 1 3 AA MM 5 22 97 274 3973 72 697 4247 0.078 NA
3 2 3 MA MM 7 22 35 310 3346 121 697 3656 0.063 NA
asym.sig
1 NA
2 NA
3 NA
#pop(vac.gl6) <- vac.gl6@other$ind.metrics$DAPI
table(pop(vac.gl6))
AA MA MM
5 7 22
#s <- gl.smearplot(glkaryo, group.pop = TRUE)
pop(vac.gl6) <- vac.gl6@other$ind.metrics$pop
Investigate the West only - all three karyotypes with geneflow
### Check the genlight object for dartR requirements
<- gl.compliance.check(vac.gl6) vac.gl6
Starting gl.compliance.check
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking coding of SNPs
SNP data scored NA, 0, 1 or 2 confirmed
Checking for population assignments
Population assignments confirmed
Checking locus metrics and flags
Recalculating locus metrics
Checking for monomorphic loci
Dataset contains monomorphic loci
Checking for loci with all missing data
Dataset contains loci with all missing dat
Checking whether individual names are unique.
Checking for individual metrics
Individual metrics confirmed
Spelling of coordinates checked and changed if necessary to
lat/lon
Completed: gl.compliance.check
@loc.names %>% duplicated %>% table vac.gl6
.
FALSE
17974
# Subset the western population ------------------------------------------------
@other$ind.metrics$ind vac.gl6
[1] SW1 bh1 BH9 MI1 CC1 CC10 SW2 bh2 BHA MI2a CC2 BHF RR1 BH3 BHB
[16] MI5 CC3 RR2a bh4 BHC MI4a BH5 BHD CC6 RR3 MI3 CC7 BHJ RR4 RR5
[31] CC8 BHK RR6 BH8
65 Levels: bh1 bh2 BH3 bh4 BH5 BH6 BH7 BH8 BH9 BHA BHB BHC BHD BHE BHF ... zoo2
<- gl.keep.pop(vac.gl6, c("bh")) west
Starting gl.keep.pop
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking for presence of nominated populations
Retaining only populations bh
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.pop
# Group individuals of the western population by karyotype
pop(west) <- west@other$ind.metrics$karyo
# Reassign karyotype to a new genlight object ----------------------------------
<- gl.keep.pop(west, c("AA", "MA", "MM")) west
Starting gl.keep.pop
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking for presence of nominated populations
Retaining only populations AA, MA, MM
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.pop
# Report private and fixed alleles for each karyotype
gl.report.pa(west)
Starting gl.report.pa
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
p1 p2 pop1 pop2 N1 N2 fixed priv1 priv2 Chao1 Chao2 totalpriv AFD asym
1 1 2 AA MA 5 4 36 416 649 1 2 1065 0.027 NA
2 1 3 AA MM 5 5 146 560 916 1 3 1476 0.043 NA
3 2 3 MA MM 4 5 28 429 599 2 3 1028 0.030 NA
asym.sig
1 NA
2 NA
3 NA
Table of private alleles and fixed differences returned
Completed: gl.report.pa
p1 p2 pop1 pop2 N1 N2 fixed priv1 priv2 Chao1 Chao2 totalpriv AFD asym
1 1 2 AA MA 5 4 36 416 649 1 2 1065 0.027 NA
2 1 3 AA MM 5 5 146 560 916 1 3 1476 0.043 NA
3 2 3 MA MM 4 5 28 429 599 2 3 1028 0.030 NA
asym.sig
1 NA
2 NA
3 NA
# Add individual metadata to specimen IDs --------------------------------------
pop(west) <-
paste(
@other$ind.metrics$karyo,
west@other$ind.metrics$id,
west@other$ind.metrics$pop
west
)
# Fixed difference analysis ----------------------------------------------------
<- gl.fixed.diff(west) fp6
Starting gl.fixed.diff
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Comparing populations for absolute fixed differences
Warning: Monomorphic loci retained, used in calculations
Comparing populations pairwise -- this may take time. Please be patient
Completed: gl.fixed.diff
$df fp6
NULL
table(pop(west))
AA AA095569 bh AA AA095623 bh AA AA095625 bh AA AA095626 bh AA AA095629 bh
1 1 1 1 1
MA AA093615 bh MA AA095559 bh MA AA095562 bh MA AA095624 bh MM AA093614 bh
1 1 1 1 1
MM AA095550 bh MM AA095554 bh MM AA095558 bh MM AA095568 bh
1 1 1 1
<- (gl.fixed.diff(west)) fp6
Starting gl.fixed.diff
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Comparing populations for absolute fixed differences
Warning: Monomorphic loci retained, used in calculations
Comparing populations pairwise -- this may take time. Please be patient
Completed: gl.fixed.diff
# Figure 2 ---------------------------------------------------------------------
<- as.matrix(gl.fixed.diff(west)$fd) res
Starting gl.fixed.diff
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Comparing populations for absolute fixed differences
Warning: Monomorphic loci retained, used in calculations
Comparing populations pairwise -- this may take time. Please be patient
Completed: gl.fixed.diff
gl.plot.heatmap(res, margins = c(10, 10))
Starting gl.plot.heatmap
Processing a data matrix
Starting gl.colors
Selected color type div
Completed: gl.colors
Completed: gl.plot.heatmap
# Table of SNPs per scaffold that align to chromosome 6 ------------------------
table(vac.gl6@other$loc.metrics["Chrom_Varanus_JD_vac_aligned_chr6_7"])
Chrom_Varanus_JD_vac_aligned_chr6_7
scaf_162 scaf_174 scaf_178 scaf_18 scaf_185 scaf_44936
0 10084 448 4056 584 1183 849
scaf_44937
770
# Fixed allele analysis for with predicted Karyotypes --------------
pop(vac.gl6) <- vac.gl6@other$ind.metrics$karyo
<- gl.keep.pop(vac.gl6, c("AA", "MA", "MM")) vac.karyo
Starting gl.keep.pop
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking for presence of nominated populations
Retaining only populations AA, MA, MM
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.pop
pop(vac.karyo) <- vac.karyo@other$ind.metrics$karyo
gl.fixed.diff(vac.karyo)
Starting gl.fixed.diff
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Comparing populations for absolute fixed differences
Warning: Monomorphic loci retained, used in calculations
Comparing populations pairwise -- this may take time. Please be patient
Completed: gl.fixed.diff
$gl
********************
*** DARTR OBJECT ***
********************
** 34 genotypes, 17,974 SNPs , size: 97.4 Mb
missing data: 122942 (=20.12 %) scored as NA
** Genetic data
@gen: list of 34 SNPbin
@ploidy: ploidy of each individual (range: 2-2)
** Additional data
@ind.names: 34 individual labels
@loc.names: 17974 locus labels
@loc.all: 17974 allele labels
@position: integer storing positions of the SNPs [within 69 base sequence]
@pop: population of each individual (group size range: 5-22)
@other: a list containing: loc.metrics, ind.metrics, latlon, loc.metrics.flags, history, verbose
@other$ind.metrics: id, species, ind, pop, sex, lat, lon, e, n, karyo, PCR1, DAPI, FISH
@other$loc.metrics: AlleleID, CloneID, AlleleSequence, TrimmedSequence, Chrom_Komodo_dragon_NCBIv1, ChromPos_Komodo_dragon_NCBIv1, AlnCnt_Komodo_dragon_NCBIv1, AlnEvalue_Komodo_dragon_NCBIv1, Chrom_Varanus_JD_infer_z, ChromPos_Varanus_JD_infer_z, AlnCnt_Varanus_JD_infer_z, AlnEvalue_Varanus_JD_infer_z, Chrom_Varanus_JD_vac_aligned_chr6_7, ChromPos_Varanus_JD_vac_aligned_chr6_7, AlnCnt_Varanus_JD_vac_aligned_chr6_7, AlnEvalue_Varanus_JD_vac_aligned_chr6_7, Chrom_Varanus_JD_infer_w, ChromPos_Varanus_JD_infer_w, AlnCnt_Varanus_JD_infer_w, AlnEvalue_Varanus_JD_infer_w, SNP, SnpPosition, CallRate, OneRatioRef, OneRatioSnp, FreqHomRef, FreqHomSnp, FreqHets, PICRef, PICSnp, AvgPIC, AvgCountRef, AvgCountSnp, RepAvg, clone, uid, rdepth, maf, monomorphs, OneRatio, PIC
@other$latlon[g]: coordinates for all individuals are attached
$fd
AA MA
MA 36
MM 97 35
$pcfd
AA MA
MA 0
MM 1 0
$nobs
AA MA MM
AA NA 10.8 21.7
MA 10.8 NA 23.3
MM 21.7 23.3 NA
$nloc
AA MA MM
AA NA 16319 16368
MA 16319 NA 17402
MM 16368 17402 NA
$expfpos
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA
$sdfpos
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA
$pval
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA
attr(,"class")
[1] "fd"
.1 <- gl.fixed.diff(vac.karyo) fpk6
Starting gl.fixed.diff
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Comparing populations for absolute fixed differences
Warning: Monomorphic loci retained, used in calculations
Comparing populations pairwise -- this may take time. Please be patient
Completed: gl.fixed.diff
.1$fd fpk6
AA MA
MA 36
MM 97 35
table(pop(vac.karyo))
AA MA MM
5 7 22
# group individuals by karyotype
<- seppop(vac.karyo)
pops
# Identify those loci that are fixed for one allele in AA and opposite allele in
# MM and MA are heterozygotes
<- which(((abs(
themany colMeans(as.matrix(pops$AA), na.rm = TRUE) - colMeans(as.matrix(pops$MM),
na.rm = TRUE)
== 2) & colMeans(as.matrix(pops$MA) == 1, na.rm = T))
))
themany
47530186-62-T/G 47510325-33-C/T 47506537-50-C/T 58564194-10-C/T 47530853-13-C/A
2475 2697 2943 3166 3178
58568357-53-C/T 47555226-11-T/A 47526162-36-G/T 58574714-5-C/G 47570599-32-G/C
4530 4876 5019 5086 5407
58564685-25-G/A 58568357-43-C/G 47571057-27-A/G 58568730-8-C/T 47534249-44-C/T
5800 5980 6197 6330 6429
58564990-43-G/C 58568170-47-G/A 47572514-30-T/C 47567360-61-C/G 47562299-15-T/A
7015 7392 7575 8894 8990
47490392-36-G/T 47557237-34-A/G 47539268-39-G/C 47584615-15-C/T 58581892-46-T/C
9300 9569 9614 9880 10097
58581176-56-T/C 58582604-32-G/A 47495246-42-A/G 58555628-11-A/C 58534605-26-C/T
10554 10746 13050 13395 14239
47572698-68-C/T 58548074-61-A/G 58548705-34-G/C 47553573-52-C/T 58552886-50-G/A
14309 14506 14826 14837 15026
58544938-37-C/T 47555650-59-C/T 47553573-53-G/T 47537859-68-G/T 58552335-7-G/A
15518 15526 15564 15699 16020
47553754-35-C/T 58554983-38-A/G 58526779-49-T/C 47495750-16-T/G 58529930-35-G/C
16725 16778 17129 17471 17788
# Supplemenarty Table 5 of SNPs sorted for fixed differences--------------------
@other$loc.metrics$TrimmedSequence[themany] vac.karyo
[1] TGCAGCTTTAAGTGTGCAATTGCCACACACTCCAATGCTTGTGTGAACTGACCTTCTCACACTCCACTG
[2] TGCAGATGCTGGGGGTGAATCTGGGCCCTCTTGCCTGCTCCCCAACTTTGCAAGGATGTGGTAGGGATT
[3] TGCAGAACCTCAGACCCACCAGGGAGCCCCAGATGGACACTCACCTCCCCCAAAGCTTCATGCCTTTCT
[4] TGCAGAAAATCGTCATTACCATAGAGGGGGAAAGGCAGCTAAAGGTGAGACAGGAAGGGCCCCGAAAGC
[5] TGCAGGCAAAAATCAAAAATTTCAAGCAAAACAATGAACTGTTTCACACATAATGGCAGGGAAGTGCAG
[6] TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
[7] TGCAGCGGAAATCTTTCAGCAAGGGCCACCCCATTCAATGGCAAAAGCCATGTGACAGACCATTAAACT
[8] TGCAGAAGTGGCGAGGACCAAAATACGGCGATCTTCGTGCCCATACACGCATTTTTAATACCTCAGCTG
[9] TGCAGCGAATGGTCCCTCAGATAAGAGCCTTAAGATGTTCCTGCAAAAGAACATGTGAGCATGAGCGGT
[10] TGCAGAATGAGAAGACCTGGTCTAGTAAATGTGTCATCAAATTCCAGAAAAGTGAAATGAACCATAGTG
[11] TGCAGCAGTGTCAGTAGATTCCACCGGTGCTTCCAACGCATAAAAAAGACAGAAATCACAGGCTATCCT
[12] TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
[13] TGCAGATTTTAAGAAGGAGGGGTTATGAAAGTGGAATGTATCATACAAGTTCTTCTAAGTAAGAGAAAT
[14] TGCAGCGTCGATCCTCCCGGGCATCCCACTCGCCCCGCCAGGCCATTTTGCCACTTCACTCGCCAGCAA
[15] TGCAGTTGCAGTACTAGCGGTAGTTTTTCAATGGGAGATTAAAGCAAGCTTGTTCATTCACATACCGTA
[16] TGCAGCAAACATCATTAACATAGTGCCTGCGGGCTGCTTTAGTGTCTGCTATCGATAGCAGAGACTAAA
[17] TGCAGAAATATGTAGAAATAAAGAAAACCGTATAAAATGATTCTAAAGGTACGATCCTGTACAAGATTA
[18] TGCAGGTGATTACAGCAACCCACCGTAAAATCAGTCTCGGGTTCAGGGAAGGCAACTACTTGGCATGGG
[19] TGCAGGTCCTAACTTCCTGCACTGACACTCTTGCAAAGACCAAGAAAATCCTGATATTCCCCATCCTCT
[20] TGCAGATGAGTTCTGTGATTATGCTCTTTATAGAATGCTCGTTTTCACCATATCAATGCTTCTAGGAGT
[21] TGCAGAAGAAGGTTGCTGGTTTGCACAACAGTAATTGCCTGAGAGTTGGTGGGCTTCCTCTCCGCTCAG
[22] TGCAGTACCTCCACCTGGCATGTATCTTTCTGTCACAATGACTGACAATAGATTGCGGACTGAACGAAG
[23] TGCAGACTTGCTACCCTGTTGAAGTGAATACCTCCCACAGACTCATTGCAAAGACATGGCGACATGCCA
[24] TGCAGACAGAAACTTCCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
[25] TGCAGGCGTCTTACACTCCAACAAAGTACTGGGTCCGACGCGATGCTTGATGTCTCGCACGGGCAACTC
[26] TGCAGTTACCATTTCCCCATCCTTATTAACCTCTAGATCAGTTACTGCTTAGGCACTACCCCCCAGTTC
[27] TGCAGACAGAAACTTTCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
[28] TGCAGCTGCCTCTATAACCTGTCTAGGCAGCCCTCACAAGACAGACTGGATTCTGCGGAAGAGGCCAAG
[29] TGCAGTCCTGGATTTTAAGTGCCTTTCTTTCTGCTTGAACAGGGAAGGCTTCTGAACCTGGTGAGATCC
[30] TGCAGAGGAGGGTAGGACAGGTTCCACTGGAACCATGCTCAGAGAACCAAGGTTCAATCTGGATGCCAG
[31] TGCAGTATGACACAGTAGAGCAGGAGTCCTTAAACTTTTTCTAAAGGGGGCCAGTTGACAGGTGCCAGC
[32] TGCAGCAGAAAATGCACTAAGAAATGCCACCAGCCAGCCATGCCTCCTCCCAGAACACTCCATTCCACG
[33] TGCAGGTGAGGAACAGGTGAAGAATGGTGGGCAAGTGCTCTTCCAGGGTAGCTTGTCCACATTTACTGA
[34] TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
[35] TGCAGGAGAGTTGCTTGCTTGCTTGTTTTCTGTACCACTCAATAGCCAAGGCTCTCTCGGACTGCTTTA
[36] TGCAGTTACAAGCAAACTGAACAACCAGAAAAGAAAACGTTGACCTAGATGTTTGCTACTCTTACTCAT
[37] TGCAGCTTGCAACCTCTAAAGGGCCGGGTGGTGTGCTCAGGCAGAAGTCCATTCCTTAACGTGAGAATA
[38] TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
[39] TGCAGAACCCTGGGCACTGTCAGCCCAGCCATTGATCTGAACTGTGGTGGCAATGGCCAATTATTCAGG
[40] TGCAGGTGTATTCAGAGTACAGTAGGCTCTCTGTTAACCGGAATTCAAACAACTGGCAAAAATAAAGAA
[41] TGCAGAGCCTGGCCTTAAGATAAACACCTCTGCTTCTCTCATCTCAACCTCTGGAATTGGGCCACATTC
[42] TGCAGCGTCTGCCTGCTGGGTCTCTTGGATTCCTCTGGATAGCTCCAGAGAGCCTCCATCTCTCCCTTC
[43] TGCAGCGCTGTTGTTAGCAAGCAAAACTAGGACAGCTTATGAGTATACGTTAGGCATGAGCGGTTCAGC
[44] TGCAGGAAAGACGGGCTCAATCTCAGCCATTCCCCTGCTAACCAGCATGATCACAGCCATCTCCCCACC
[45] TGCAGAGGGGCTCTAATCTTGATCTGCGCATATCAGGGTTTTTAAAAACAAATATACACTTCGTTGCAT
183938 Levels: TGCAGAAAAAAATCATTAAATTATTCCTGGACTGTACCATTTCCGGGCACAAGTAGAGGGATAAAGGTC ...
<- vac.karyo@other$loc.metrics[themany, ]
hits hits
AlleleID CloneID
39780 47530186|F|0-62:T>G-62:T>G 47530186
43463 47510325|F|0-33:C>T-33:C>T 47510325
46899 47506537|F|0-50:C>T-50:C>T 47506537
50141 58564194|F|0-10:C>T-10:C>T 58564194
50274 47530853|F|0-13:C>A-13:C>A 47530853
70972 58568357|F|0-53:C>T-53:C>T 58568357
76139 47555226|F|0-11:T>A-11:T>A 47555226
78208 47526162|F|0-36:G>T-36:G>T 47526162
79343 58574714|F|0-5:C>G-5:C>G 58574714
84162 47570599|F|0-32:G>C-32:G>C 47570599
90579 58564685|F|0-25:G>A-25:G>A 58564685
92968 58568357|F|0-43:C>G-43:C>G 58568357
95906 47571057|F|0-27:A>G-27:A>G 47571057
97972 58568730|F|0-8:C>T-8:C>T 58568730
99146 47534249|F|0-44:C>T-44:C>T 47534249
109333 58564990|F|0-43:G>C-43:G>C 58564990
115084 58568170|F|0-47:G>A-47:G>A 58568170
117389 47572514|F|0-30:T>C-30:T>C 47572514
139763 47567360|F|0-61:C>G-61:C>G 47567360
140882 47562299|F|0-15:T>A-15:T>A 47562299
146039 47490392|F|0-36:G>T-36:G>T 47490392
149807 47557237|F|0-34:A>G-34:A>G 47557237
150336 47539268|F|0-39:G>C-39:G>C 47539268
153551 47584615|F|0-15:C>T-15:C>T 47584615
156726 58581892|F|0-46:T>C-46:T>C 58581892
163050 58581176|F|0-56:T>C-56:T>C 58581176
166096 58582604|F|0-32:G>A-32:G>A 58582604
205464 47495246|F|0-42:A>G-42:A>G 47495246
211014 58555628|F|0-11:A>C-11:A>C 58555628
224380 58534605|F|0-26:C>T-26:C>T 58534605
225655 47572698|F|0-68:C>T-68:C>T 47572698
229644 58548074|F|0-61:A>G-61:A>G 58548074
235514 58548705|F|0-34:G>C-34:G>C 58548705
235706 47553573|F|0-52:C>T-52:C>T 47553573
238914 58552886|F|0-50:G>A-50:G>A 58552886
247070 58544938|F|0-37:C>T-37:C>T 58544938
247268 47555650|F|0-59:C>T-59:C>T 47555650
248073 47553573|F|0-53:G>T-53:G>T 47553573
250541 47537859|F|0-68:G>T-68:G>T 47537859
255750 58552335|F|0-7:G>A-7:G>A 58552335
268769 47553754|F|0-35:C>T-35:C>T 47553754
269844 58554983|F|0-38:A>G-38:A>G 58554983
277181 58526779|F|0-49:T>C-49:T>C 58526779
283756 47495750|F|0-16:T>G-16:T>G 47495750
295061 58529930|F|0-35:G>C-35:G>C 58529930
AlleleSequence
39780 TGCAGCTTTAAGTGTGCAATTGCCACACACTCCAATGCTTGTGTGAACTGACCTTCTCACACTCCACTG
43463 TGCAGATGCTGGGGGTGAATCTGGGCCCTCTTGCCTGCTCCCCAACTTTGCAAGGATGTGGTAGGGATT
46899 TGCAGAACCTCAGACCCACCAGGGAGCCCCAGATGGACACTCACCTCCCCCAAAGCTTCATGCCTTTCT
50141 TGCAGAAAATCGTCATTACCATAGAGGGGGAAAGGCAGCTAAAGGTGAGACAGGAAGGGCCCCGAAAGC
50274 TGCAGGCAAAAATCAAAAATTTCAAGCAAAACAATGAACTGTTTCACACATAATGGCAGGGAAGTGCAG
70972 TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
76139 TGCAGCGGAAATCTTTCAGCAAGGGCCACCCCATTCAATGGCAAAAGCCATGTGACAGACCATTAAACT
78208 TGCAGAAGTGGCGAGGACCAAAATACGGCGATCTTCGTGCCCATACACGCATTTTTAATACCTCAGCTG
79343 TGCAGCGAATGGTCCCTCAGATAAGAGCCTTAAGATGTTCCTGCAAAAGAACATGTGAGCATGAGCGGT
84162 TGCAGAATGAGAAGACCTGGTCTAGTAAATGTGTCATCAAATTCCAGAAAAGTGAAATGAACCATAGTG
90579 TGCAGCAGTGTCAGTAGATTCCACCGGTGCTTCCAACGCATAAAAAAGACAGAAATCACAGGCTATCCT
92968 TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
95906 TGCAGATTTTAAGAAGGAGGGGTTATGAAAGTGGAATGTATCATACAAGTTCTTCTAAGTAAGAGAAAT
97972 TGCAGCGTCGATCCTCCCGGGCATCCCACTCGCCCCGCCAGGCCATTTTGCCACTTCACTCGCCAGCAA
99146 TGCAGTTGCAGTACTAGCGGTAGTTTTTCAATGGGAGATTAAAGCAAGCTTGTTCATTCACATACCGTA
109333 TGCAGCAAACATCATTAACATAGTGCCTGCGGGCTGCTTTAGTGTCTGCTATCGATAGCAGAGACTAAA
115084 TGCAGAAATATGTAGAAATAAAGAAAACCGTATAAAATGATTCTAAAGGTACGATCCTGTACAAGATTA
117389 TGCAGGTGATTACAGCAACCCACCGTAAAATCAGTCTCGGGTTCAGGGAAGGCAACTACTTGGCATGGG
139763 TGCAGGTCCTAACTTCCTGCACTGACACTCTTGCAAAGACCAAGAAAATCCTGATATTCCCCATCCTCT
140882 TGCAGATGAGTTCTGTGATTATGCTCTTTATAGAATGCTCGTTTTCACCATATCAATGCTTCTAGGAGT
146039 TGCAGAAGAAGGTTGCTGGTTTGCACAACAGTAATTGCCTGAGAGTTGGTGGGCTTCCTCTCCGCTCAG
149807 TGCAGTACCTCCACCTGGCATGTATCTTTCTGTCACAATGACTGACAATAGATTGCGGACTGAACGAAG
150336 TGCAGACTTGCTACCCTGTTGAAGTGAATACCTCCCACAGACTCATTGCAAAGACATGGCGACATGCCA
153551 TGCAGACAGAAACTTCCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
156726 TGCAGGCGTCTTACACTCCAACAAAGTACTGGGTCCGACGCGATGCTTGATGTCTCGCACGGGCAACTC
163050 TGCAGTTACCATTTCCCCATCCTTATTAACCTCTAGATCAGTTACTGCTTAGGCACTACCCCCCAGTTC
166096 TGCAGACAGAAACTTTCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
205464 TGCAGCTGCCTCTATAACCTGTCTAGGCAGCCCTCACAAGACAGACTGGATTCTGCGGAAGAGGCCAAG
211014 TGCAGTCCTGGATTTTAAGTGCCTTTCTTTCTGCTTGAACAGGGAAGGCTTCTGAACCTGGTGAGATCC
224380 TGCAGAGGAGGGTAGGACAGGTTCCACTGGAACCATGCTCAGAGAACCAAGGTTCAATCTGGATGCCAG
225655 TGCAGTATGACACAGTAGAGCAGGAGTCCTTAAACTTTTTCTAAAGGGGGCCAGTTGACAGGTGCCAGC
229644 TGCAGCAGAAAATGCACTAAGAAATGCCACCAGCCAGCCATGCCTCCTCCCAGAACACTCCATTCCACG
235514 TGCAGGTGAGGAACAGGTGAAGAATGGTGGGCAAGTGCTCTTCCAGGGTAGCTTGTCCACATTTACTGA
235706 TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
238914 TGCAGGAGAGTTGCTTGCTTGCTTGTTTTCTGTACCACTCAATAGCCAAGGCTCTCTCGGACTGCTTTA
247070 TGCAGTTACAAGCAAACTGAACAACCAGAAAAGAAAACGTTGACCTAGATGTTTGCTACTCTTACTCAT
247268 TGCAGCTTGCAACCTCTAAAGGGCCGGGTGGTGTGCTCAGGCAGAAGTCCATTCCTTAACGTGAGAATA
248073 TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
250541 TGCAGAACCCTGGGCACTGTCAGCCCAGCCATTGATCTGAACTGTGGTGGCAATGGCCAATTATTCAGG
255750 TGCAGGTGTATTCAGAGTACAGTAGGCTCTCTGTTAACCGGAATTCAAACAACTGGCAAAAATAAAGAA
268769 TGCAGAGCCTGGCCTTAAGATAAACACCTCTGCTTCTCTCATCTCAACCTCTGGAATTGGGCCACATTC
269844 TGCAGCGTCTGCCTGCTGGGTCTCTTGGATTCCTCTGGATAGCTCCAGAGAGCCTCCATCTCTCCCTTC
277181 TGCAGCGCTGTTGTTAGCAAGCAAAACTAGGACAGCTTATGAGTATACGTTAGGCATGAGCGGTTCAGC
283756 TGCAGGAAAGACGGGCTCAATCTCAGCCATTCCCCTGCTAACCAGCATGATCACAGCCATCTCCCCACC
295061 TGCAGAGGGGCTCTAATCTTGATCTGCGCATATCAGGGTTTTTAAAAACAAATATACACTTCGTTGCAT
TrimmedSequence
39780 TGCAGCTTTAAGTGTGCAATTGCCACACACTCCAATGCTTGTGTGAACTGACCTTCTCACACTCCACTG
43463 TGCAGATGCTGGGGGTGAATCTGGGCCCTCTTGCCTGCTCCCCAACTTTGCAAGGATGTGGTAGGGATT
46899 TGCAGAACCTCAGACCCACCAGGGAGCCCCAGATGGACACTCACCTCCCCCAAAGCTTCATGCCTTTCT
50141 TGCAGAAAATCGTCATTACCATAGAGGGGGAAAGGCAGCTAAAGGTGAGACAGGAAGGGCCCCGAAAGC
50274 TGCAGGCAAAAATCAAAAATTTCAAGCAAAACAATGAACTGTTTCACACATAATGGCAGGGAAGTGCAG
70972 TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
76139 TGCAGCGGAAATCTTTCAGCAAGGGCCACCCCATTCAATGGCAAAAGCCATGTGACAGACCATTAAACT
78208 TGCAGAAGTGGCGAGGACCAAAATACGGCGATCTTCGTGCCCATACACGCATTTTTAATACCTCAGCTG
79343 TGCAGCGAATGGTCCCTCAGATAAGAGCCTTAAGATGTTCCTGCAAAAGAACATGTGAGCATGAGCGGT
84162 TGCAGAATGAGAAGACCTGGTCTAGTAAATGTGTCATCAAATTCCAGAAAAGTGAAATGAACCATAGTG
90579 TGCAGCAGTGTCAGTAGATTCCACCGGTGCTTCCAACGCATAAAAAAGACAGAAATCACAGGCTATCCT
92968 TGCAGAGAGGGACATGGGGGAAGAATTCCAGCAGCCAGCCAAGCCTGGTATGACTCCTGATCCCATGAT
95906 TGCAGATTTTAAGAAGGAGGGGTTATGAAAGTGGAATGTATCATACAAGTTCTTCTAAGTAAGAGAAAT
97972 TGCAGCGTCGATCCTCCCGGGCATCCCACTCGCCCCGCCAGGCCATTTTGCCACTTCACTCGCCAGCAA
99146 TGCAGTTGCAGTACTAGCGGTAGTTTTTCAATGGGAGATTAAAGCAAGCTTGTTCATTCACATACCGTA
109333 TGCAGCAAACATCATTAACATAGTGCCTGCGGGCTGCTTTAGTGTCTGCTATCGATAGCAGAGACTAAA
115084 TGCAGAAATATGTAGAAATAAAGAAAACCGTATAAAATGATTCTAAAGGTACGATCCTGTACAAGATTA
117389 TGCAGGTGATTACAGCAACCCACCGTAAAATCAGTCTCGGGTTCAGGGAAGGCAACTACTTGGCATGGG
139763 TGCAGGTCCTAACTTCCTGCACTGACACTCTTGCAAAGACCAAGAAAATCCTGATATTCCCCATCCTCT
140882 TGCAGATGAGTTCTGTGATTATGCTCTTTATAGAATGCTCGTTTTCACCATATCAATGCTTCTAGGAGT
146039 TGCAGAAGAAGGTTGCTGGTTTGCACAACAGTAATTGCCTGAGAGTTGGTGGGCTTCCTCTCCGCTCAG
149807 TGCAGTACCTCCACCTGGCATGTATCTTTCTGTCACAATGACTGACAATAGATTGCGGACTGAACGAAG
150336 TGCAGACTTGCTACCCTGTTGAAGTGAATACCTCCCACAGACTCATTGCAAAGACATGGCGACATGCCA
153551 TGCAGACAGAAACTTCCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
156726 TGCAGGCGTCTTACACTCCAACAAAGTACTGGGTCCGACGCGATGCTTGATGTCTCGCACGGGCAACTC
163050 TGCAGTTACCATTTCCCCATCCTTATTAACCTCTAGATCAGTTACTGCTTAGGCACTACCCCCCAGTTC
166096 TGCAGACAGAAACTTTCCCTTAGTCACTCTTCGGTCATCGTTTCCCAAGGAGGACCGTTTCCTACCCCT
205464 TGCAGCTGCCTCTATAACCTGTCTAGGCAGCCCTCACAAGACAGACTGGATTCTGCGGAAGAGGCCAAG
211014 TGCAGTCCTGGATTTTAAGTGCCTTTCTTTCTGCTTGAACAGGGAAGGCTTCTGAACCTGGTGAGATCC
224380 TGCAGAGGAGGGTAGGACAGGTTCCACTGGAACCATGCTCAGAGAACCAAGGTTCAATCTGGATGCCAG
225655 TGCAGTATGACACAGTAGAGCAGGAGTCCTTAAACTTTTTCTAAAGGGGGCCAGTTGACAGGTGCCAGC
229644 TGCAGCAGAAAATGCACTAAGAAATGCCACCAGCCAGCCATGCCTCCTCCCAGAACACTCCATTCCACG
235514 TGCAGGTGAGGAACAGGTGAAGAATGGTGGGCAAGTGCTCTTCCAGGGTAGCTTGTCCACATTTACTGA
235706 TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
238914 TGCAGGAGAGTTGCTTGCTTGCTTGTTTTCTGTACCACTCAATAGCCAAGGCTCTCTCGGACTGCTTTA
247070 TGCAGTTACAAGCAAACTGAACAACCAGAAAAGAAAACGTTGACCTAGATGTTTGCTACTCTTACTCAT
247268 TGCAGCTTGCAACCTCTAAAGGGCCGGGTGGTGTGCTCAGGCAGAAGTCCATTCCTTAACGTGAGAATA
248073 TGCAGACTGCCACAGAGGAGTTGCAGCTCAACACAAGTGGGATTAGTTTGCACGGCATGAGCGGTTCAG
250541 TGCAGAACCCTGGGCACTGTCAGCCCAGCCATTGATCTGAACTGTGGTGGCAATGGCCAATTATTCAGG
255750 TGCAGGTGTATTCAGAGTACAGTAGGCTCTCTGTTAACCGGAATTCAAACAACTGGCAAAAATAAAGAA
268769 TGCAGAGCCTGGCCTTAAGATAAACACCTCTGCTTCTCTCATCTCAACCTCTGGAATTGGGCCACATTC
269844 TGCAGCGTCTGCCTGCTGGGTCTCTTGGATTCCTCTGGATAGCTCCAGAGAGCCTCCATCTCTCCCTTC
277181 TGCAGCGCTGTTGTTAGCAAGCAAAACTAGGACAGCTTATGAGTATACGTTAGGCATGAGCGGTTCAGC
283756 TGCAGGAAAGACGGGCTCAATCTCAGCCATTCCCCTGCTAACCAGCATGATCACAGCCATCTCCCCACC
295061 TGCAGAGGGGCTCTAATCTTGATCTGCGCATATCAGGGTTTTTAAAAACAAATATACACTTCGTTGCAT
Chrom_Komodo_dragon_NCBIv1 ChromPos_Komodo_dragon_NCBIv1
39780 SJPD01000043.1_scaffold34_1 6324171
43463 SJPD01000031.1_scaffold24 16405804
46899 SJPD01000043.1_scaffold34_1 375374
50141 0
50274 SJPD01000031.1_scaffold24 11788603
70972 SJPD01000043.1_scaffold34_1 4580642
76139 SJPD01000043.1_scaffold34_1 7862835
78208 SJPD01000043.1_scaffold34_1 7938504
79343 SJPD01000031.1_scaffold24 10159958
84162 0
90579 0
92968 SJPD01000043.1_scaffold34_1 4580642
95906 0
97972 SJPD01000043.1_scaffold34_1 8699162
99146 SJPD01000031.1_scaffold24 7573219
109333 SJPD01000011.1_scaffold8 35951871
115084 0
117389 0
139763 SJPD01000031.1_scaffold24 10990793
140882 SJPD01000043.1_scaffold34_1 1337171
146039 SJPD01000031.1_scaffold24 10187800
149807 0
150336 SJPD01000031.1_scaffold24 16184650
153551 SJPD01000043.1_scaffold34_1 6564774
156726 0
163050 SJPD01000043.1_scaffold34_1 5237185
166096 0
205464 SJPD01000043.1_scaffold34_1 5621116
211014 0
224380 SJPD01000001.1_scaffold1 114327324
225655 SJPD01000031.1_scaffold24 12436780
229644 SJPD01000043.1_scaffold34_1 8251535
235514 SJPD01000031.1_scaffold24 17012576
235706 0
238914 SJPD01000031.1_scaffold24 8349203
247070 SJPD01000043.1_scaffold34_1 2154387
247268 SJPD01000001.1_scaffold1 113487419
248073 0
250541 SJPD01000031.1_scaffold24 17827486
255750 SJPD01000043.1_scaffold34_1 7881951
268769 SJPD01000031.1_scaffold24 16078793
269844 0
277181 SJPD01000031.1_scaffold24 12273432
283756 SJPD01000003.1_scaffold3 20807026
295061 0
AlnCnt_Komodo_dragon_NCBIv1 AlnEvalue_Komodo_dragon_NCBIv1
39780 1 2.21e-21
43463 1 1.32e-23
46899 1 2.19e-26
50141 0 9.99e+02
50274 1 2.84e-25
70972 2 1.32e-23
76139 1 2.84e-25
78208 1 2.84e-25
79343 1 6.14e-22
84162 0 9.99e+02
90579 0 9.99e+02
92968 2 1.32e-23
95906 0 9.99e+02
97972 1 3.67e-24
99146 1 1.71e-22
109333 10 1.33e-18
115084 0 9.99e+02
117389 0 9.99e+02
139763 1 6.10e-27
140882 1 2.84e-25
146039 1 3.70e-19
149807 0 9.99e+02
150336 1 1.71e-22
153551 1 6.14e-22
156726 0 9.99e+02
163050 1 6.14e-22
166096 0 9.99e+02
205464 1 1.33e-18
211014 0 9.99e+02
224380 1 4.78e-18
225655 1 6.10e-27
229644 1 3.72e-14
235514 1 7.89e-26
235706 0 9.99e+02
238914 1 6.10e-27
247070 1 6.10e-27
247268 1 6.10e-27
248073 0 9.99e+02
250541 1 2.84e-25
255750 1 6.10e-27
268769 1 6.10e-27
269844 0 9.99e+02
277181 1 8.00e-16
283756 1 6.14e-22
295061 0 9.99e+02
Chrom_Varanus_JD_infer_z ChromPos_Varanus_JD_infer_z
39780 0
43463 0
46899 0
50141 0
50274 0
70972 0
76139 0
78208 0
79343 0
84162 0
90579 0
92968 0
95906 0
97972 0
99146 0
109333 0
115084 0
117389 0
139763 0
140882 0
146039 0
149807 0
150336 0
153551 0
156726 0
163050 0
166096 0
205464 0
211014 0
224380 0
225655 0
229644 0
235514 0
235706 0
238914 0
247070 0
247268 0
248073 0
250541 0
255750 0
268769 0
269844 0
277181 0
283756 0
295061 0
AlnCnt_Varanus_JD_infer_z AlnEvalue_Varanus_JD_infer_z
39780 0 999
43463 0 999
46899 0 999
50141 0 999
50274 0 999
70972 0 999
76139 0 999
78208 0 999
79343 0 999
84162 0 999
90579 0 999
92968 0 999
95906 0 999
97972 0 999
99146 0 999
109333 0 999
115084 0 999
117389 0 999
139763 0 999
140882 0 999
146039 0 999
149807 0 999
150336 0 999
153551 0 999
156726 0 999
163050 0 999
166096 0 999
205464 0 999
211014 0 999
224380 0 999
225655 0 999
229644 0 999
235514 0 999
235706 0 999
238914 0 999
247070 0 999
247268 0 999
248073 0 999
250541 0 999
255750 0 999
268769 0 999
269844 0 999
277181 0 999
283756 0 999
295061 0 999
Chrom_Varanus_JD_vac_aligned_chr6_7
39780 scaf_185
43463 scaf_44937
46899 scaf_174
50141 scaf_44937
50274 scaf_178
70972 scaf_44936
76139 scaf_185
78208 scaf_185
79343 scaf_178
84162 scaf_178
90579 scaf_178
92968 scaf_44936
95906 scaf_44936
97972 scaf_185
99146 scaf_178
109333 scaf_178
115084 scaf_178
117389 scaf_185
139763 scaf_178
140882 scaf_174
146039 scaf_178
149807 scaf_178
150336 scaf_44937
153551 scaf_185
156726 scaf_185
163050 scaf_44936
166096 scaf_185
205464 scaf_185
211014 scaf_185
224380 scaf_178
225655 scaf_178
229644 scaf_185
235514 scaf_44937
235706 scaf_178
238914 scaf_178
247070 scaf_174
247268 scaf_178
248073 scaf_178
250541 scaf_44937
255750 scaf_185
268769 scaf_44937
269844 scaf_174
277181 scaf_178
283756 scaf_162
295061 scaf_44936
ChromPos_Varanus_JD_vac_aligned_chr6_7
39780 4535602
43463 1879053
46899 422662
50141 3057519
50274 17533506
70972 2239338
76139 3058006
78208 2982377
79343 15921589
84162 9531917
90579 8627028
92968 2239338
95906 626638
97972 2280435
99146 13361437
109333 16092795
115084 4799684
117389 3789392
139763 16744777
140882 1374442
146039 15947400
149807 12871115
150336 1656795
153551 4311211
156726 1611346
163050 2878581
166096 4311211
205464 5197021
211014 4008715
224380 959068
225655 18176452
229644 2688663
235514 2488100
235706 14281639
238914 14121110
247070 2162400
247268 1932940
248073 14281639
250541 3299754
255750 3038917
268769 1552585
269844 2162466
277181 18014400
283756 27088034
295061 2382323
AlnCnt_Varanus_JD_vac_aligned_chr6_7
39780 1
43463 1
46899 1
50141 1
50274 1
70972 1
76139 1
78208 1
79343 1
84162 1
90579 1
92968 1
95906 1
97972 1
99146 1
109333 6
115084 1
117389 1
139763 1
140882 1
146039 1
149807 1
150336 1
153551 1
156726 1
163050 1
166096 1
205464 1
211014 1
224380 1
225655 2
229644 1
235514 1
235706 1
238914 1
247070 1
247268 1
248073 1
250541 1
255750 1
268769 1
269844 1
277181 1
283756 2
295061 1
AlnEvalue_Varanus_JD_vac_aligned_chr6_7 Chrom_Varanus_JD_infer_w
39780 7.96e-30
43463 7.96e-30
46899 7.96e-30
50141 7.96e-30
50274 7.96e-30
70972 7.96e-30
76139 7.96e-30
78208 7.96e-30
79343 2.88e-24
84162 7.96e-30
90579 7.96e-30
92968 7.96e-30
95906 7.96e-30
97972 7.96e-30
99146 3.70e-28
109333 3.70e-28
115084 3.70e-28
117389 3.70e-28
139763 7.96e-30
140882 7.96e-30
146039 7.96e-30
149807 3.70e-28
150336 7.96e-30
153551 7.96e-30
156726 7.96e-30
163050 8.02e-25
166096 3.70e-28
205464 3.70e-28
211014 3.70e-28
224380 7.96e-30
225655 7.96e-30
229644 7.96e-30
235514 7.96e-30
235706 2.88e-24
238914 7.96e-30
247070 7.96e-30
247268 7.96e-30
248073 2.88e-24
250541 7.96e-30
255750 7.96e-30
268769 7.96e-30
269844 3.70e-28
277181 2.25e-20
283756 2.93e-14
295061 2.25e-20
ChromPos_Varanus_JD_infer_w AlnCnt_Varanus_JD_infer_w
39780 0 0
43463 0 0
46899 0 0
50141 0 0
50274 0 0
70972 0 0
76139 0 0
78208 0 0
79343 0 0
84162 0 0
90579 0 0
92968 0 0
95906 0 0
97972 0 0
99146 0 0
109333 0 0
115084 0 0
117389 0 0
139763 0 0
140882 0 0
146039 0 0
149807 0 0
150336 0 0
153551 0 0
156726 0 0
163050 0 0
166096 0 0
205464 0 0
211014 0 0
224380 0 0
225655 0 0
229644 0 0
235514 0 0
235706 0 0
238914 0 0
247070 0 0
247268 0 0
248073 0 0
250541 0 0
255750 0 0
268769 0 0
269844 0 0
277181 0 0
283756 0 0
295061 0 0
AlnEvalue_Varanus_JD_infer_w SNP SnpPosition CallRate OneRatioRef
39780 999 62:T>G 62 1.000000 0.8529412
43463 999 33:C>T 33 0.970588 0.8484848
46899 999 50:C>T 50 0.823529 0.7500000
50141 999 10:C>T 10 0.852941 0.8275862
50274 999 13:C>A 13 1.000000 0.8529412
70972 999 53:C>T 53 0.823529 0.8571429
76139 999 11:T>A 11 0.941176 0.8437500
78208 999 36:G>T 36 0.794118 0.8148148
79343 999 5:C>G 5 0.735294 0.8000000
84162 999 32:G>C 32 0.941176 0.8437500
90579 999 25:G>A 25 0.735294 0.8000000
92968 999 43:C>G 43 0.823529 0.8571429
95906 999 27:A>G 27 0.823529 0.8214286
97972 999 8:C>T 8 0.941176 0.8437500
99146 999 44:C>T 44 0.823529 0.4285714
109333 999 43:G>C 43 0.823529 0.7857143
115084 999 47:G>A 47 0.323529 0.3636364
117389 999 30:T>C 30 0.588235 0.6000000
139763 999 61:C>G 61 1.000000 0.8529412
140882 999 15:T>A 15 1.000000 0.8529412
146039 999 36:G>T 36 1.000000 0.8529412
149807 999 34:A>G 34 0.735294 0.6400000
150336 999 39:G>C 39 0.794118 0.8148148
153551 999 15:C>T 15 0.764706 0.6153846
156726 999 46:T>C 46 0.735294 0.6800000
163050 999 56:T>C 56 0.558824 0.4210526
166096 999 32:G>A 32 0.411765 0.8571429
205464 999 42:A>G 42 0.970588 0.3333333
211014 999 11:A>C 11 0.911765 0.2903226
224380 999 26:C>T 26 0.794118 0.8148148
225655 999 68:C>T 68 0.852941 0.7931034
229644 999 61:A>G 61 0.941176 0.8125000
235514 999 34:G>C 34 0.970588 0.7575758
235706 999 52:C>T 52 0.911765 0.8709677
238914 999 50:G>A 50 0.823529 0.8214286
247070 999 37:C>T 37 0.823529 0.7857143
247268 999 59:C>T 59 0.882353 0.8333333
248073 999 53:G>T 53 0.941176 0.8437500
250541 999 68:G>T 68 0.941176 0.8750000
255750 999 7:G>A 7 0.823529 0.6428571
268769 999 35:C>T 35 0.852941 0.9310345
269844 999 38:A>G 38 0.735294 0.6800000
277181 999 49:T>C 49 0.588235 0.4500000
283756 999 16:T>G 16 0.352941 0.5000000
295061 999 35:G>C 35 0.382353 0.8461538
OneRatioSnp FreqHomRef FreqHomSnp FreqHets PICRef PICSnp
39780 0.3529412 0.6470588 0.14705882 0.20588235 0.2508651 0.4567474
43463 0.3636364 0.6363636 0.15151515 0.21212121 0.2571166 0.4628099
46899 0.4285714 0.5714286 0.25000000 0.17857143 0.3750000 0.4897959
50141 0.3103448 0.6896552 0.17241379 0.13793103 0.2853746 0.4280618
50274 0.2647059 0.7352941 0.14705882 0.11764706 0.2508651 0.3892734
70972 0.2142857 0.7857143 0.14285714 0.07142857 0.2448980 0.3367347
76139 0.2500000 0.7500000 0.15625000 0.09375000 0.2636719 0.3750000
78208 0.2962963 0.7037037 0.18518519 0.11111111 0.3017833 0.4170096
79343 0.4000000 0.6000000 0.20000000 0.20000000 0.3200000 0.4800000
84162 0.3750000 0.6250000 0.15625000 0.21875000 0.2636719 0.4687500
90579 0.3600000 0.6400000 0.20000000 0.16000000 0.3200000 0.4608000
92968 0.2500000 0.7500000 0.14285714 0.10714286 0.2448980 0.3750000
95906 0.2857143 0.7142857 0.17857143 0.10714286 0.2933673 0.4081633
97972 0.3125000 0.6875000 0.15625000 0.15625000 0.2636719 0.4296875
99146 0.8214286 0.1785714 0.57142857 0.25000000 0.4897959 0.2933673
109333 0.3214286 0.6785714 0.21428571 0.10714286 0.3367347 0.4362245
115084 0.7272727 0.2727273 0.63636364 0.09090909 0.4628099 0.3966942
117389 0.5500000 0.4500000 0.40000000 0.15000000 0.4800000 0.4950000
139763 0.3235294 0.6764706 0.14705882 0.17647059 0.2508651 0.4377163
140882 0.2058824 0.7941176 0.14705882 0.05882353 0.2508651 0.3269896
146039 0.2941176 0.7058824 0.14705882 0.14705882 0.2508651 0.4152249
149807 0.4800000 0.5200000 0.36000000 0.12000000 0.4608000 0.4992000
150336 0.3333333 0.6666667 0.18518519 0.14814815 0.3017833 0.4444444
153551 0.4615385 0.5384615 0.38461538 0.07692308 0.4733728 0.4970414
156726 0.4400000 0.5600000 0.32000000 0.12000000 0.4352000 0.4928000
163050 0.6315789 0.3684211 0.57894737 0.05263158 0.4875346 0.4653740
166096 0.2142857 0.7857143 0.14285714 0.07142857 0.2448980 0.3367347
205464 0.7878788 0.2121212 0.66666667 0.12121212 0.4444444 0.3342516
211014 0.8064516 0.1935484 0.70967742 0.09677419 0.4120708 0.3121748
224380 0.4074074 0.5925926 0.18518519 0.22222222 0.3017833 0.4828532
225655 0.3793103 0.6206897 0.20689655 0.17241379 0.3281807 0.4708680
229644 0.2812500 0.7187500 0.18750000 0.09375000 0.3046875 0.4042969
235514 0.3636364 0.6363636 0.24242424 0.12121212 0.3673095 0.4628099
235706 0.1612903 0.8387097 0.12903226 0.03225806 0.2247659 0.2705515
238914 0.4285714 0.5714286 0.17857143 0.25000000 0.2933673 0.4897959
247070 0.3928571 0.6071429 0.21428571 0.17857143 0.3367347 0.4770408
247268 0.4000000 0.6000000 0.16666667 0.23333333 0.2777778 0.4800000
248073 0.2812500 0.7187500 0.15625000 0.12500000 0.2636719 0.4042969
250541 0.2500000 0.7500000 0.12500000 0.12500000 0.2187500 0.3750000
255750 0.3928571 0.6071429 0.35714286 0.03571429 0.4591837 0.4770408
268769 0.1724138 0.8275862 0.06896552 0.10344828 0.1284185 0.2853746
269844 0.4000000 0.6000000 0.32000000 0.08000000 0.4352000 0.4800000
277181 0.6500000 0.3500000 0.55000000 0.10000000 0.4950000 0.4550000
283756 0.5833333 0.4166667 0.50000000 0.08333333 0.5000000 0.4861111
295061 0.2307692 0.7692308 0.15384615 0.07692308 0.2603550 0.3550296
AvgPIC AvgCountRef AvgCountSnp RepAvg clone uid rdepth
39780 0.3538062 6.39256 5.26549 0.992308 47530186 47530186-62 7.3
43463 0.3599633 6.98502 5.53659 0.992188 47510325 47510325-33 8.0
46899 0.4323980 3.43028 3.11765 0.990741 47506537 47506537-50 3.6
50141 0.3567182 4.65323 3.38182 0.963636 58564194 58564194-10 4.9
50274 0.3200692 5.37970 3.55556 0.991379 47530853 47530853-13 5.5
70972 0.2908163 3.15179 2.61290 0.970588 58568357 58568357-53 3.4
76139 0.3193359 7.28685 6.76923 1.000000 47555226 47555226-11 7.6
78208 0.3593964 9.28947 8.66000 0.991071 47526162 47526162-36 9.8
79343 0.4000000 10.18595 9.39583 1.000000 58574714 58574714-5 11.0
84162 0.3662109 7.08808 5.27273 0.982143 47570599 47570599-32 7.5
90579 0.3904000 3.51101 3.31034 0.980392 58564685 58564685-25 3.7
92968 0.3099490 3.15179 2.68421 0.989583 58568357 58568357-43 3.2
95906 0.3507653 4.36697 3.50000 0.990385 47571057 47571057-27 4.5
97972 0.3466797 3.17225 3.04878 0.959184 58568730 58568730-8 3.5
99146 0.3915816 4.96226 6.76136 0.981818 47534249 47534249-44 7.4
109333 0.3864796 3.43085 3.18182 1.000000 58564990 58564990-43 3.6
115084 0.4297521 5.90370 4.81967 0.987500 58568170 58568170-47 5.9
117389 0.4875000 4.94949 4.31683 1.000000 47572514 47572514-30 5.0
139763 0.3442907 11.39085 9.09804 1.000000 47567360 47567360-61 12.3
140882 0.2889273 12.79012 7.41053 0.992424 47562299 47562299-15 13.1
146039 0.3330450 9.19712 7.79825 0.984127 47490392 47490392-36 10.2
149807 0.4800000 13.12022 9.26087 0.991071 47557237 47557237-34 13.9
150336 0.3731139 16.48000 15.55000 1.000000 47539268 47539268-39 18.1
153551 0.4852071 5.48864 5.31343 0.979167 47584615 47584615-15 5.8
156726 0.4640000 3.57059 3.45714 0.980769 58581892 58581892-46 4.0
163050 0.4764543 14.11538 11.69748 1.000000 58581176 58581176-56 12.4
166096 0.2908163 5.31343 5.20000 0.982143 58582604 58582604-32 5.5
205464 0.3893480 6.20000 6.56981 0.991228 47495246 47495246-42 7.0
211014 0.3621228 4.43333 4.13139 1.000000 58555628 58555628-11 4.3
224380 0.3923182 4.41758 3.63768 0.975410 58534605 58534605-26 4.7
225655 0.3995244 3.65238 2.97059 0.953704 47572698 47572698-68 4.0
229644 0.3544922 3.46185 2.94872 0.956897 58548074 58548074-61 3.6
235514 0.4150597 3.06186 2.65179 0.928571 58548705 58548705-34 3.5
235706 0.2476587 13.40948 9.10769 0.981818 47553573 47553573-52 14.0
238914 0.3915816 4.14583 3.35849 0.981481 58552886 58552886-50 4.5
247070 0.4068878 4.12121 3.25210 0.937500 58544938 58544938-37 4.4
247268 0.3788889 4.54237 3.67442 0.981132 47555650 47555650-59 4.9
248073 0.3339844 13.40948 13.02439 0.990566 47553573 47553573-53 14.7
250541 0.2968750 6.90000 5.63636 1.000000 47537859 47537859-68 7.1
255750 0.4681122 3.20879 2.91304 0.976744 58552335 58552335-7 3.5
268769 0.2068966 6.87374 6.16667 0.988372 47553754 47553754-35 7.0
269844 0.4576000 3.95302 3.07229 0.980392 58554983 58554983-38 4.1
277181 0.4750000 16.03704 9.48322 1.000000 58526779 58526779-49 10.6
283756 0.4930556 6.35000 5.91228 0.983333 47495750 47495750-16 6.3
295061 0.3076923 5.02151 4.66667 0.981481 58529930 58529930-35 5.1
maf monomorphs OneRatio PIC
39780 0.2500 NA NA NA
43463 0.2576 NA NA NA
46899 0.3393 NA NA NA
50141 0.2414 NA NA NA
50274 0.2059 NA NA NA
70972 0.1786 NA NA NA
76139 0.2031 NA NA NA
78208 0.2407 NA NA NA
79343 0.3000 NA NA NA
84162 0.2656 NA NA NA
90579 0.2800 NA NA NA
92968 0.1964 NA NA NA
95906 0.2321 NA NA NA
97972 0.2344 NA NA NA
99146 0.3036 NA NA NA
109333 0.2679 NA NA NA
115084 0.3182 NA NA NA
117389 0.4750 NA NA NA
139763 0.2353 NA NA NA
140882 0.1765 NA NA NA
146039 0.2206 NA NA NA
149807 0.4200 NA NA NA
150336 0.2593 NA NA NA
153551 0.4231 NA NA NA
156726 0.3800 NA NA NA
163050 0.3947 NA NA NA
166096 0.1786 NA NA NA
205464 0.2727 NA NA NA
211014 0.2419 NA NA NA
224380 0.2963 NA NA NA
225655 0.2931 NA NA NA
229644 0.2344 NA NA NA
235514 0.3030 NA NA NA
235706 0.1452 NA NA NA
238914 0.3036 NA NA NA
247070 0.3036 NA NA NA
247268 0.2833 NA NA NA
248073 0.2188 NA NA NA
250541 0.1875 NA NA NA
255750 0.3750 NA NA NA
268769 0.1207 NA NA NA
269844 0.3600 NA NA NA
277181 0.4000 NA NA NA
283756 0.4583 NA NA NA
295061 0.1923 NA NA NA
<- vac.karyo
vac.pops
# SNPs aligning to scaffolds of interest following fixed allele analysis
table(vac.pops@other$loc.metrics[themany, "Chrom_Varanus_JD_vac_aligned_chr6_7"])
scaf_162 scaf_174 scaf_178 scaf_18 scaf_185 scaf_44936
0 1 4 17 0 12 5
scaf_44937
6
# create an index for only those loci on scaf_178-------------------------------
<- which(vac.pops@other$loc.metrics$Chrom_Varanus_JD_vac_aligned_chr6_7 == "scaf_178")
index_178
# check the index to verify genlight object
index_178
[1] 4 14 19 24 42 45 46 47 54 55 62 66
[13] 67 74 78 82 97 104 119 120 122 125 129 132
[25] 141 143 144 155 157 161 174 175 193 194 199 217
[37] 221 228 229 234 239 240 241 246 247 248 262 266
[49] 269 275 276 281 293 294 299 306 307 312 316 320
[61] 330 333 334 335 349 352 363 367 374 381 385 394
[73] 396 406 407 417 430 439 440 441 442 444 445 446
[85] 447 457 484 495 496 498 502 503 507 524 548 551
[97] 562 571 576 578 579 586 610 612 614 617 618 620
[109] 627 628 645 653 661 663 664 665 666 670 671 672
[121] 677 682 684 685 686 687 693 695 696 698 699 701
[133] 702 704 707 708 709 710 711 720 721 722 723 728
[145] 731 732 734 738 743 746 748 750 751 754 757 760
[157] 765 766 768 770 771 776 778 782 784 787 794 799
[169] 827 828 835 839 850 863 864 865 867 868 869 870
[181] 872 879 881 882 883 886 892 893 897 898 900 915
[193] 920 924 929 930 936 937 938 951 952 953 955 959
[205] 963 969 970 971 973 974 975 977 981 990 991 992
[217] 993 1012 1014 1015 1016 1021 1022 1036 1040 1041 1044 1045
[229] 1046 1047 1049 1051 1053 1055 1057 1059 1061 1063 1075 1081
[241] 1082 1089 1092 1105 1108 1111 1121 1129 1139 1140 1141 1146
[253] 1155 1157 1162 1176 1177 1178 1180 1183 1184 1194 1201 1204
[265] 1205 1211 1212 1215 1230 1239 1240 1243 1244 1252 1254 1263
[277] 1265 1266 1268 1269 1270 1283 1290 1291 1292 1294 1295 1297
[289] 1301 1305 1306 1307 1310 1314 1315 1323 1327 1328 1329 1335
[301] 1336 1337 1349 1350 1354 1364 1365 1366 1367 1368 1369 1371
[313] 1375 1379 1383 1393 1395 1402 1404 1406 1411 1413 1414 1418
[325] 1428 1429 1441 1444 1445 1446 1448 1452 1453 1463 1469 1472
[337] 1473 1475 1477 1479 1488 1492 1493 1501 1505 1507 1508 1509
[349] 1511 1515 1531 1536 1537 1545 1546 1549 1561 1562 1579 1580
[361] 1586 1589 1594 1595 1597 1601 1604 1605 1606 1607 1608 1609
[373] 1611 1616 1645 1649 1652 1653 1675 1677 1681 1682 1683 1684
[385] 1691 1696 1698 1705 1717 1718 1719 1720 1743 1748 1750 1763
[397] 1766 1767 1774 1778 1782 1788 1790 1796 1799 1801 1803 1819
[409] 1820 1821 1823 1825 1826 1828 1831 1836 1837 1838 1839 1840
[421] 1848 1857 1858 1862 1863 1866 1874 1876 1882 1883 1884 1891
[433] 1892 1893 1909 1911 1912 1917 1922 1926 1927 1934 1935 1937
[445] 1938 1949 1950 1951 1952 1953 1955 1957 1958 1961 1967 1970
[457] 1974 1979 1980 1981 1982 1994 2003 2005 2017 2021 2022 2024
[469] 2025 2028 2037 2040 2041 2042 2075 2077 2081 2087 2089 2092
[481] 2094 2097 2099 2101 2103 2104 2107 2111 2117 2119 2120 2121
[493] 2125 2126 2130 2133 2147 2148 2165 2175 2182 2186 2195 2198
[505] 2199 2201 2202 2204 2214 2218 2219 2221 2227 2230 2231 2232
[517] 2233 2235 2238 2240 2243 2254 2255 2256 2264 2268 2275 2284
[529] 2292 2295 2312 2315 2316 2317 2326 2333 2334 2340 2341 2342
[541] 2343 2344 2346 2348 2358 2367 2375 2378 2380 2382 2394 2397
[553] 2406 2413 2414 2421 2424 2425 2434 2436 2446 2447 2450 2451
[565] 2452 2454 2457 2458 2459 2467 2468 2470 2474 2476 2478 2484
[577] 2490 2491 2494 2499 2500 2503 2505 2511 2525 2527 2528 2529
[589] 2530 2539 2540 2559 2569 2572 2573 2574 2575 2576 2582 2583
[601] 2586 2587 2589 2590 2593 2596 2597 2599 2602 2609 2618 2626
[613] 2628 2632 2633 2643 2647 2648 2649 2670 2675 2676 2685 2687
[625] 2702 2712 2716 2724 2726 2732 2733 2747 2752 2760 2761 2762
[637] 2763 2764 2766 2798 2802 2806 2808 2812 2813 2828 2832 2837
[649] 2841 2848 2850 2852 2858 2859 2866 2871 2872 2879 2885 2886
[661] 2887 2888 2917 2920 2922 2925 2926 2941 2942 2950 2956 2957
[673] 2963 2986 2996 2999 3002 3006 3016 3025 3028 3029 3030 3031
[685] 3035 3041 3049 3050 3052 3054 3055 3056 3058 3068 3076 3078
[697] 3082 3086 3100 3108 3109 3125 3126 3129 3130 3131 3141 3142
[709] 3143 3144 3147 3154 3159 3161 3167 3168 3171 3172 3175 3176
[721] 3178 3180 3184 3186 3188 3193 3196 3198 3203 3207 3208 3211
[733] 3212 3213 3221 3223 3229 3232 3236 3237 3241 3242 3244 3250
[745] 3257 3263 3267 3273 3275 3276 3293 3299 3301 3307 3315 3317
[757] 3339 3345 3348 3349 3350 3353 3354 3386 3387 3400 3402 3406
[769] 3412 3416 3427 3430 3432 3434 3436 3440 3447 3450 3464 3467
[781] 3473 3495 3498 3499 3512 3514 3515 3522 3524 3528 3533 3537
[793] 3538 3541 3543 3549 3557 3558 3561 3562 3566 3581 3585 3596
[805] 3597 3613 3617 3626 3628 3630 3637 3638 3640 3644 3648 3652
[817] 3653 3654 3655 3657 3659 3665 3670 3677 3686 3689 3691 3695
[829] 3724 3730 3731 3734 3742 3746 3761 3763 3766 3767 3781 3784
[841] 3791 3793 3794 3796 3797 3801 3810 3811 3816 3819 3820 3824
[853] 3828 3829 3833 3834 3837 3839 3843 3847 3851 3852 3854 3863
[865] 3865 3877 3885 3894 3898 3909 3910 3911 3918 3922 3925 3927
[877] 3932 3938 3941 3942 3947 3950 3954 3955 3956 3957 3958 3959
[889] 3967 3973 3980 3981 3983 3985 3986 3987 3992 3994 3999 4000
[901] 4016 4017 4018 4021 4026 4032 4035 4038 4044 4047 4048 4049
[913] 4051 4055 4057 4058 4059 4060 4071 4073 4076 4084 4089 4100
[925] 4101 4103 4107 4111 4113 4124 4125 4127 4128 4133 4136 4137
[937] 4144 4145 4157 4160 4161 4164 4165 4167 4171 4176 4181 4185
[949] 4192 4193 4194 4198 4201 4204 4214 4220 4223 4226 4227 4235
[961] 4239 4248 4249 4250 4255 4261 4267 4274 4281 4283 4284 4292
[973] 4298 4302 4305 4331 4336 4337 4349 4351 4356 4359 4364 4369
[985] 4372 4375 4386 4387 4389 4390 4394 4395 4403 4425 4428 4430
[997] 4439 4440 4442 4445 4450 4455 4457 4459 4460 4461 4466 4467
[1009] 4472 4476 4479 4481 4482 4485 4491 4495 4505 4507 4519 4520
[1021] 4523 4540 4543 4551 4560 4561 4569 4570 4572 4577 4579 4582
[1033] 4583 4589 4590 4591 4593 4595 4606 4611 4616 4626 4629 4639
[1045] 4642 4648 4650 4651 4656 4657 4666 4667 4673 4674 4678 4679
[1057] 4680 4693 4719 4722 4728 4738 4755 4758 4763 4767 4779 4784
[1069] 4790 4793 4795 4800 4808 4812 4813 4837 4838 4840 4841 4846
[1081] 4852 4853 4854 4862 4864 4874 4878 4885 4891 4906 4915 4930
[1093] 4936 4940 4943 4954 4955 4956 4962 4966 4972 4980 4982 5002
[1105] 5003 5004 5005 5007 5009 5015 5020 5027 5029 5036 5044 5046
[1117] 5056 5059 5068 5069 5086 5087 5088 5092 5101 5109 5117 5119
[1129] 5132 5136 5138 5141 5143 5151 5152 5159 5160 5181 5188 5194
[1141] 5204 5215 5218 5226 5228 5240 5242 5244 5255 5269 5270 5271
[1153] 5273 5274 5279 5281 5282 5284 5291 5314 5315 5316 5335 5337
[1165] 5343 5351 5353 5354 5358 5360 5367 5368 5388 5389 5390 5391
[1177] 5397 5405 5406 5407 5413 5414 5419 5437 5440 5441 5446 5450
[1189] 5456 5457 5458 5459 5460 5462 5464 5465 5468 5475 5480 5481
[1201] 5486 5488 5489 5493 5494 5495 5503 5511 5515 5516 5522 5529
[1213] 5534 5536 5538 5541 5551 5552 5553 5568 5574 5578 5579 5581
[1225] 5583 5585 5594 5596 5614 5623 5625 5627 5628 5631 5633 5642
[1237] 5647 5649 5651 5652 5653 5655 5657 5662 5664 5669 5670 5673
[1249] 5675 5679 5689 5691 5698 5703 5704 5717 5718 5724 5727 5728
[1261] 5730 5732 5743 5745 5747 5752 5758 5764 5765 5768 5770 5777
[1273] 5781 5782 5787 5793 5794 5796 5800 5804 5808 5810 5817 5818
[1285] 5826 5831 5833 5834 5835 5836 5843 5844 5847 5850 5854 5860
[1297] 5861 5865 5871 5872 5880 5881 5888 5890 5895 5903 5904 5906
[1309] 5920 5923 5924 5928 5934 5935 5936 5942 5946 5947 5949 5950
[1321] 5952 5953 5960 5961 5972 5979 5983 5986 5988 5989 5990 5992
[1333] 6009 6010 6020 6022 6025 6034 6035 6037 6044 6045 6046 6047
[1345] 6052 6053 6059 6060 6068 6069 6072 6081 6082 6083 6088 6094
[1357] 6105 6107 6108 6115 6119 6122 6130 6132 6136 6140 6159 6160
[1369] 6161 6170 6186 6188 6190 6195 6199 6214 6218 6223 6227 6228
[1381] 6231 6250 6252 6253 6259 6272 6274 6276 6285 6286 6287 6303
[1393] 6316 6319 6320 6323 6324 6329 6335 6340 6344 6347 6355 6359
[1405] 6363 6364 6370 6385 6397 6399 6401 6403 6406 6417 6423 6428
[1417] 6429 6437 6438 6446 6449 6451 6454 6460 6469 6474 6479 6482
[1429] 6490 6491 6496 6497 6498 6499 6506 6511 6515 6518 6526 6539
[1441] 6542 6544 6545 6550 6561 6577 6580 6582 6586 6589 6590 6600
[1453] 6604 6610 6611 6612 6621 6625 6626 6629 6634 6635 6644 6658
[1465] 6659 6662 6666 6670 6674 6677 6689 6690 6712 6718 6723 6726
[1477] 6727 6728 6732 6735 6738 6739 6742 6748 6749 6754 6756 6760
[1489] 6763 6765 6771 6787 6801 6804 6815 6816 6817 6820 6824 6829
[1501] 6830 6835 6847 6852 6856 6857 6858 6862 6864 6865 6871 6876
[1513] 6881 6886 6894 6895 6913 6918 6921 6925 6935 6941 6942 6944
[1525] 6947 6949 6953 6956 6960 6963 6964 6965 6974 6986 6992 6993
[1537] 6994 6996 6997 6999 7002 7004 7015 7018 7021 7025 7028 7029
[1549] 7030 7038 7056 7057 7058 7066 7071 7072 7077 7083 7101 7102
[1561] 7105 7107 7108 7109 7116 7133 7147 7153 7155 7160 7167 7173
[1573] 7176 7177 7181 7185 7200 7205 7209 7210 7218 7232 7235 7252
[1585] 7263 7266 7267 7272 7273 7276 7278 7279 7283 7292 7302 7307
[1597] 7309 7312 7315 7318 7319 7322 7326 7328 7331 7351 7354 7355
[1609] 7357 7361 7370 7371 7373 7374 7378 7382 7386 7392 7396 7401
[1621] 7402 7403 7416 7417 7424 7425 7428 7431 7439 7440 7441 7446
[1633] 7458 7459 7465 7470 7471 7482 7485 7491 7494 7498 7500 7504
[1645] 7510 7511 7514 7519 7520 7523 7532 7535 7536 7539 7546 7547
[1657] 7548 7550 7552 7565 7569 7577 7583 7588 7595 7599 7601 7602
[1669] 7607 7609 7611 7622 7623 7624 7626 7636 7637 7642 7645 7646
[1681] 7648 7662 7668 7669 7670 7671 7677 7679 7685 7686 7689 7692
[1693] 7693 7708 7709 7713 7721 7726 7728 7734 7736 7739 7748 7751
[1705] 7752 7755 7757 7767 7773 7785 7793 7796 7797 7800 7802 7805
[1717] 7806 7808 7812 7814 7818 7824 7828 7835 7837 7839 7851 7858
[1729] 7859 7862 7865 7866 7868 7869 7872 7875 7883 7886 7891 7894
[1741] 7898 7906 7911 7922 7928 7933 7942 7944 7948 7949 7950 7953
[1753] 7960 7961 7971 7972 7973 7990 7993 7994 7997 7999 8001 8005
[1765] 8008 8015 8016 8017 8021 8023 8025 8026 8027 8039 8044 8049
[1777] 8053 8058 8069 8070 8072 8080 8083 8088 8090 8093 8097 8108
[1789] 8123 8124 8137 8139 8140 8142 8144 8146 8151 8154 8157 8159
[1801] 8161 8170 8175 8178 8183 8192 8194 8208 8219 8228 8240 8241
[1813] 8245 8257 8262 8268 8272 8273 8274 8277 8293 8295 8296 8299
[1825] 8317 8319 8321 8325 8327 8329 8339 8342 8345 8347 8350 8351
[1837] 8352 8357 8358 8359 8363 8381 8387 8388 8394 8402 8406 8409
[1849] 8411 8415 8428 8445 8447 8451 8455 8461 8462 8464 8466 8467
[1861] 8469 8471 8479 8485 8490 8495 8499 8501 8502 8503 8505 8506
[1873] 8517 8523 8525 8527 8532 8533 8534 8540 8543 8546 8551 8552
[1885] 8561 8563 8568 8575 8589 8590 8595 8598 8607 8609 8612 8613
[1897] 8616 8623 8626 8631 8636 8637 8643 8645 8646 8649 8650 8654
[1909] 8656 8663 8664 8669 8674 8680 8686 8688 8690 8696 8697 8701
[1921] 8702 8706 8710 8713 8715 8719 8733 8736 8739 8740 8741 8744
[1933] 8745 8751 8753 8755 8778 8779 8784 8785 8788 8793 8801 8802
[1945] 8803 8806 8808 8811 8812 8815 8816 8817 8821 8827 8831 8836
[1957] 8842 8845 8846 8854 8861 8865 8867 8871 8872 8881 8883 8884
[1969] 8894 8897 8899 8900 8904 8905 8911 8912 8914 8915 8916 8927
[1981] 8931 8932 8933 8934 8943 8947 8948 8949 8950 8954 8955 8956
[1993] 8965 8968 8973 8974 8975 8977 8979 8980 8987 8988 8994 8999
[2005] 9007 9008 9009 9019 9021 9026 9036 9039 9046 9047 9048 9051
[2017] 9052 9053 9054 9056 9057 9068 9070 9072 9073 9074 9077 9078
[2029] 9079 9088 9089 9091 9094 9112 9113 9116 9117 9124 9125 9138
[2041] 9139 9143 9152 9154 9155 9157 9158 9159 9160 9166 9175 9177
[2053] 9178 9190 9195 9196 9197 9198 9199 9209 9210 9212 9213 9216
[2065] 9220 9221 9226 9228 9233 9240 9241 9244 9249 9251 9255 9258
[2077] 9261 9265 9273 9283 9286 9288 9298 9300 9309 9319 9330 9331
[2089] 9336 9342 9343 9344 9354 9361 9366 9368 9369 9376 9381 9387
[2101] 9388 9394 9395 9400 9405 9416 9420 9441 9449 9453 9468 9470
[2113] 9482 9486 9488 9491 9492 9495 9502 9507 9516 9526 9531 9533
[2125] 9541 9544 9545 9546 9551 9552 9553 9556 9557 9558 9561 9564
[2137] 9567 9569 9584 9586 9589 9590 9594 9608 9609 9610 9624 9632
[2149] 9637 9645 9648 9661 9666 9667 9671 9679 9680 9682 9684 9690
[2161] 9691 9695 9696 9698 9703 9704 9713 9715 9716 9717 9719 9721
[2173] 9726 9730 9734 9735 9737 9738 9743 9751 9757 9758 9761 9772
[2185] 9773 9778 9789 9792 9797 9800 9801 9808 9817 9821 9823 9831
[2197] 9838 9839 9850 9853 9857 9867 9878 9882 9889 9895 9901 9904
[2209] 9909 9925 9926 9931 9939 9943 9946 9949 9952 9956 9957 9960
[2221] 9966 9968 9969 9971 9973 9981 9991 9992 10004 10006 10008 10015
[2233] 10024 10025 10026 10030 10033 10037 10040 10048 10052 10077 10080 10081
[2245] 10104 10107 10109 10111 10121 10125 10142 10143 10148 10151 10153 10154
[2257] 10162 10170 10174 10175 10179 10187 10197 10206 10210 10211 10212 10215
[2269] 10223 10224 10225 10230 10240 10241 10245 10250 10251 10253 10255 10256
[2281] 10257 10261 10262 10264 10265 10269 10276 10283 10284 10286 10295 10299
[2293] 10300 10306 10318 10324 10333 10334 10347 10350 10351 10353 10386 10390
[2305] 10394 10398 10401 10415 10416 10423 10424 10426 10429 10431 10437 10439
[2317] 10442 10450 10458 10466 10470 10474 10475 10476 10477 10481 10482 10490
[2329] 10494 10499 10501 10506 10508 10509 10513 10517 10518 10519 10521 10524
[2341] 10537 10541 10550 10567 10568 10581 10589 10600 10601 10603 10609 10632
[2353] 10634 10636 10660 10663 10670 10675 10677 10679 10681 10682 10684 10689
[2365] 10691 10692 10695 10697 10701 10702 10703 10707 10708 10719 10720 10721
[2377] 10722 10729 10732 10733 10737 10740 10743 10744 10745 10748 10757 10762
[2389] 10763 10770 10781 10783 10788 10797 10798 10799 10801 10802 10803 10804
[2401] 10806 10812 10813 10814 10816 10819 10836 10842 10846 10852 10863 10866
[2413] 10869 10873 10874 10877 10878 10880 10891 10901 10912 10913 10916 10920
[2425] 10923 10928 10929 10931 10934 10936 10937 10939 10945 10948 10954 10956
[2437] 10958 10961 10962 10964 10966 10968 10971 10975 10976 10978 10981 10986
[2449] 10987 10993 10994 10995 11001 11005 11006 11007 11009 11011 11014 11018
[2461] 11022 11026 11028 11032 11033 11040 11044 11050 11056 11058 11059 11066
[2473] 11075 11076 11077 11083 11085 11088 11089 11099 11103 11104 11117 11119
[2485] 11120 11121 11122 11123 11131 11134 11137 11138 11142 11143 11144 11147
[2497] 11176 11177 11179 11184 11185 11194 11196 11198 11202 11208 11212 11213
[2509] 11214 11215 11218 11223 11229 11233 11235 11241 11243 11249 11256 11258
[2521] 11260 11262 11267 11278 11291 11307 11312 11316 11317 11318 11319 11321
[2533] 11327 11328 11331 11332 11333 11335 11337 11340 11343 11354 11356 11358
[2545] 11359 11360 11362 11366 11367 11382 11388 11392 11393 11394 11397 11398
[2557] 11399 11404 11406 11409 11413 11415 11424 11425 11429 11432 11435 11436
[2569] 11438 11439 11440 11445 11448 11450 11453 11455 11458 11463 11468 11470
[2581] 11478 11480 11482 11483 11489 11490 11497 11499 11510 11513 11518 11519
[2593] 11529 11530 11531 11534 11536 11544 11545 11546 11555 11567 11569 11571
[2605] 11572 11573 11575 11576 11579 11583 11584 11593 11601 11603 11605 11612
[2617] 11615 11620 11640 11643 11649 11656 11664 11670 11675 11681 11690 11691
[2629] 11692 11696 11716 11719 11723 11724 11728 11732 11736 11740 11741 11743
[2641] 11746 11747 11748 11754 11760 11762 11766 11768 11769 11771 11773 11778
[2653] 11782 11789 11797 11801 11804 11805 11807 11808 11809 11813 11817 11820
[2665] 11834 11835 11836 11845 11847 11855 11859 11864 11866 11879 11886 11895
[2677] 11897 11898 11900 11903 11913 11917 11922 11923 11924 11925 11926 11929
[2689] 11931 11939 11942 11944 11950 11951 11962 11972 11973 11974 11975 11976
[2701] 11980 11981 11982 11983 11987 11988 11989 11996 12001 12003 12011 12012
[2713] 12026 12031 12037 12041 12047 12056 12061 12065 12071 12078 12090 12091
[2725] 12094 12098 12099 12100 12111 12119 12124 12128 12140 12141 12149 12150
[2737] 12151 12153 12158 12168 12173 12176 12178 12179 12180 12181 12184 12185
[2749] 12187 12190 12193 12199 12201 12203 12209 12212 12213 12229 12230 12233
[2761] 12238 12247 12248 12249 12256 12264 12265 12267 12274 12275 12276 12283
[2773] 12290 12291 12296 12300 12301 12320 12326 12334 12335 12345 12347 12357
[2785] 12360 12361 12364 12370 12371 12372 12377 12382 12386 12395 12396 12400
[2797] 12401 12404 12408 12410 12418 12426 12427 12430 12442 12443 12446 12449
[2809] 12450 12451 12452 12456 12457 12460 12463 12479 12480 12482 12496 12502
[2821] 12513 12520 12540 12545 12549 12564 12567 12584 12587 12588 12596 12597
[2833] 12600 12604 12629 12633 12634 12635 12641 12646 12653 12655 12656 12658
[2845] 12659 12660 12676 12681 12682 12688 12693 12700 12701 12707 12727 12731
[2857] 12741 12745 12747 12751 12755 12764 12766 12782 12783 12785 12790 12794
[2869] 12806 12807 12812 12813 12815 12816 12818 12839 12842 12845 12847 12850
[2881] 12856 12857 12858 12862 12864 12865 12867 12875 12880 12887 12894 12912
[2893] 12916 12924 12931 12939 12945 12961 12964 12965 12967 12968 12969 12972
[2905] 12973 12975 12988 12991 13001 13003 13005 13007 13010 13014 13019 13020
[2917] 13023 13036 13037 13049 13055 13057 13061 13065 13066 13071 13075 13076
[2929] 13079 13092 13105 13107 13108 13117 13125 13129 13130 13132 13135 13142
[2941] 13144 13147 13156 13157 13158 13160 13166 13168 13181 13182 13183 13191
[2953] 13193 13196 13200 13201 13202 13205 13209 13212 13216 13218 13223 13242
[2965] 13246 13247 13250 13270 13273 13274 13278 13287 13292 13294 13297 13299
[2977] 13307 13310 13311 13321 13323 13325 13331 13334 13339 13344 13348 13349
[2989] 13351 13352 13353 13355 13359 13360 13369 13372 13373 13378 13381 13382
[3001] 13388 13392 13400 13403 13407 13413 13415 13416 13423 13425 13429 13431
[3013] 13432 13433 13436 13441 13442 13443 13448 13449 13465 13466 13470 13471
[3025] 13479 13496 13498 13499 13500 13502 13510 13513 13522 13526 13530 13534
[3037] 13546 13551 13558 13561 13563 13564 13569 13573 13574 13576 13581 13585
[3049] 13586 13598 13606 13617 13623 13627 13629 13632 13633 13641 13647 13649
[3061] 13650 13652 13656 13657 13661 13676 13677 13680 13688 13689 13690 13694
[3073] 13695 13699 13700 13704 13709 13710 13719 13720 13728 13729 13730 13731
[3085] 13736 13740 13752 13753 13754 13763 13764 13770 13785 13788 13796 13800
[3097] 13801 13803 13809 13813 13816 13823 13824 13832 13833 13834 13836 13838
[3109] 13842 13845 13846 13847 13848 13857 13865 13869 13870 13876 13878 13881
[3121] 13893 13897 13904 13914 13915 13918 13922 13932 13934 13937 13941 13947
[3133] 13960 13962 13967 13969 13970 13973 13974 13976 13978 13981 13982 13989
[3145] 13992 13999 14001 14002 14013 14016 14017 14019 14025 14030 14032 14042
[3157] 14044 14045 14052 14060 14064 14069 14078 14080 14081 14082 14083 14086
[3169] 14094 14100 14102 14107 14109 14117 14120 14125 14127 14131 14132 14134
[3181] 14137 14140 14141 14143 14144 14156 14157 14161 14163 14170 14176 14177
[3193] 14182 14186 14187 14192 14194 14197 14198 14204 14206 14210 14214 14216
[3205] 14218 14219 14220 14225 14226 14236 14238 14239 14244 14247 14264 14268
[3217] 14272 14274 14278 14279 14280 14287 14292 14309 14311 14318 14328 14330
[3229] 14331 14335 14336 14337 14354 14367 14369 14372 14378 14382 14385 14387
[3241] 14392 14393 14394 14399 14400 14404 14413 14419 14420 14422 14430 14436
[3253] 14445 14447 14450 14461 14464 14465 14473 14476 14481 14486 14493 14496
[3265] 14497 14499 14501 14505 14509 14512 14514 14519 14524 14525 14526 14528
[3277] 14529 14531 14533 14542 14544 14548 14549 14552 14556 14559 14570 14580
[3289] 14583 14587 14589 14591 14610 14613 14623 14625 14639 14646 14665 14668
[3301] 14672 14673 14674 14675 14678 14683 14685 14686 14689 14692 14698 14699
[3313] 14701 14705 14706 14718 14727 14731 14735 14736 14738 14740 14742 14744
[3325] 14745 14756 14761 14770 14773 14783 14787 14788 14791 14800 14801 14809
[3337] 14812 14815 14816 14823 14837 14842 14848 14853 14857 14860 14865 14866
[3349] 14867 14871 14879 14880 14882 14896 14900 14903 14924 14926 14928 14929
[3361] 14933 14936 14947 14954 14958 14959 14963 14976 14990 14992 14993 15000
[3373] 15007 15009 15010 15013 15026 15028 15030 15032 15038 15042 15050 15055
[3385] 15063 15067 15070 15072 15079 15081 15092 15093 15096 15100 15114 15115
[3397] 15116 15121 15126 15137 15141 15143 15151 15152 15159 15160 15163 15171
[3409] 15174 15176 15177 15184 15189 15191 15195 15202 15203 15205 15206 15210
[3421] 15213 15221 15229 15232 15233 15239 15240 15256 15262 15263 15265 15269
[3433] 15278 15282 15295 15311 15314 15315 15318 15322 15324 15326 15328 15332
[3445] 15335 15337 15345 15347 15348 15351 15363 15365 15369 15370 15386 15388
[3457] 15397 15408 15410 15419 15423 15427 15432 15436 15437 15452 15456 15463
[3469] 15467 15469 15475 15478 15487 15490 15491 15496 15505 15512 15513 15519
[3481] 15523 15526 15527 15530 15532 15535 15536 15544 15546 15547 15549 15554
[3493] 15560 15564 15565 15569 15570 15571 15572 15574 15576 15577 15584 15590
[3505] 15591 15596 15599 15600 15617 15619 15621 15622 15625 15632 15635 15643
[3517] 15648 15650 15651 15659 15660 15674 15676 15679 15693 15700 15701 15702
[3529] 15703 15707 15713 15727 15729 15739 15741 15748 15753 15758 15759 15761
[3541] 15762 15768 15769 15772 15774 15779 15784 15785 15788 15794 15798 15802
[3553] 15803 15806 15811 15812 15825 15826 15830 15835 15838 15840 15842 15843
[3565] 15846 15851 15852 15855 15860 15863 15879 15887 15892 15893 15906 15908
[3577] 15910 15911 15922 15933 15934 15935 15938 15939 15940 15943 15958 15959
[3589] 15963 15967 15968 15977 15980 15981 15984 15990 16006 16008 16021 16023
[3601] 16024 16029 16031 16034 16036 16037 16043 16045 16049 16053 16057 16058
[3613] 16060 16068 16070 16071 16076 16078 16081 16088 16094 16105 16115 16122
[3625] 16125 16126 16128 16130 16132 16133 16134 16135 16137 16140 16143 16145
[3637] 16146 16150 16158 16169 16172 16178 16181 16192 16199 16207 16211 16212
[3649] 16213 16218 16226 16235 16238 16239 16251 16253 16264 16274 16276 16278
[3661] 16287 16292 16295 16302 16308 16309 16311 16314 16315 16324 16328 16329
[3673] 16330 16331 16334 16338 16340 16341 16349 16350 16355 16359 16362 16363
[3685] 16366 16368 16369 16370 16386 16387 16388 16391 16392 16405 16414 16426
[3697] 16431 16438 16440 16441 16447 16451 16455 16457 16459 16460 16478 16480
[3709] 16485 16486 16490 16493 16501 16503 16506 16507 16509 16510 16516 16519
[3721] 16522 16526 16527 16528 16543 16546 16548 16551 16553 16554 16555 16557
[3733] 16560 16562 16563 16570 16571 16583 16586 16589 16591 16595 16602 16604
[3745] 16613 16617 16631 16633 16636 16637 16641 16649 16650 16652 16653 16656
[3757] 16660 16667 16674 16677 16681 16688 16694 16695 16699 16700 16701 16704
[3769] 16709 16714 16722 16723 16724 16731 16732 16739 16744 16752 16753 16762
[3781] 16765 16770 16772 16783 16784 16787 16789 16793 16795 16797 16806 16807
[3793] 16809 16813 16814 16815 16816 16821 16830 16832 16835 16838 16840 16841
[3805] 16848 16859 16869 16873 16885 16886 16893 16898 16902 16904 16907 16909
[3817] 16919 16921 16928 16934 16936 16946 16948 16967 16971 16973 16975 16976
[3829] 16980 16987 16995 16996 17000 17004 17009 17014 17019 17020 17021 17025
[3841] 17030 17031 17032 17037 17052 17055 17057 17059 17064 17067 17083 17094
[3853] 17098 17103 17106 17121 17123 17129 17140 17145 17146 17148 17149 17156
[3865] 17159 17160 17164 17165 17169 17170 17172 17174 17177 17182 17184 17189
[3877] 17191 17193 17194 17200 17201 17210 17219 17220 17226 17229 17233 17239
[3889] 17252 17258 17270 17280 17283 17293 17296 17306 17309 17311 17312 17324
[3901] 17325 17331 17341 17342 17352 17353 17356 17359 17363 17367 17370 17375
[3913] 17376 17377 17378 17379 17381 17390 17394 17401 17412 17415 17416 17421
[3925] 17423 17440 17445 17447 17448 17452 17462 17464 17466 17473 17478 17487
[3937] 17493 17496 17498 17499 17501 17504 17506 17515 17516 17518 17524 17526
[3949] 17528 17530 17535 17541 17542 17543 17545 17547 17548 17549 17554 17559
[3961] 17560 17570 17572 17575 17576 17579 17582 17586 17591 17597 17613 17614
[3973] 17616 17617 17618 17619 17620 17623 17625 17628 17635 17645 17650 17652
[3985] 17653 17662 17665 17666 17667 17668 17671 17674 17676 17680 17683 17685
[3997] 17692 17693 17695 17696 17701 17702 17708 17712 17714 17719 17725 17729
[4009] 17744 17754 17755 17756 17758 17761 17771 17772 17780 17783 17790 17792
[4021] 17793 17805 17806 17811 17823 17824 17825 17827 17838 17839 17842 17852
[4033] 17856 17869 17871 17877 17879 17881 17888 17899 17917 17924 17926 17932
[4045] 17934 17938 17939 17947 17950 17953 17955 17956 17957 17968 17969 17972
# combine the indexed loci with loci names
<- locNames(vac.pops)[index_178]
index_178names
# make a new genlight object with only those indexed loci
<- gl.keep.loc(vac.pops, loc.list = index_178names) vac.gl_178
Starting gl.keep.loc
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
List of loci to keep has been specified
Deleting all but the specified loci
Completed: gl.keep.loc
<- gl.compliance.check(vac.gl_178) vac.gl_178
Starting gl.compliance.check
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking coding of SNPs
SNP data scored NA, 0, 1 or 2 confirmed
Checking for population assignments
Population assignments confirmed
Checking locus metrics and flags
Recalculating locus metrics
Checking for monomorphic loci
Dataset contains monomorphic loci
Checking for loci with all missing data
Dataset contains loci with all missing dat
Checking whether individual names are unique.
Checking for individual metrics
Individual metrics confirmed
Spelling of coordinates checked and changed if necessary to
lat/lon
Completed: gl.compliance.check
vac.gl_178
********************
*** DARTR OBJECT ***
********************
** 34 genotypes, 4,056 SNPs , size: 91.2 Mb
missing data: 27488 (=19.93 %) scored as NA
** Genetic data
@gen: list of 34 SNPbin
@ploidy: ploidy of each individual (range: 2-2)
** Additional data
@ind.names: 34 individual labels
@loc.names: 4056 locus labels
@loc.all: 4056 allele labels
@position: integer storing positions of the SNPs [within 69 base sequence]
@pop: population of each individual (group size range: 5-22)
@other: a list containing: loc.metrics, ind.metrics, latlon, loc.metrics.flags, history, verbose
@other$ind.metrics: id, species, ind, pop, sex, lat, lon, e, n, karyo, PCR1, DAPI, FISH
@other$loc.metrics: AlleleID, CloneID, AlleleSequence, TrimmedSequence, Chrom_Komodo_dragon_NCBIv1, ChromPos_Komodo_dragon_NCBIv1, AlnCnt_Komodo_dragon_NCBIv1, AlnEvalue_Komodo_dragon_NCBIv1, Chrom_Varanus_JD_infer_z, ChromPos_Varanus_JD_infer_z, AlnCnt_Varanus_JD_infer_z, AlnEvalue_Varanus_JD_infer_z, Chrom_Varanus_JD_vac_aligned_chr6_7, ChromPos_Varanus_JD_vac_aligned_chr6_7, AlnCnt_Varanus_JD_vac_aligned_chr6_7, AlnEvalue_Varanus_JD_vac_aligned_chr6_7, Chrom_Varanus_JD_infer_w, ChromPos_Varanus_JD_infer_w, AlnCnt_Varanus_JD_infer_w, AlnEvalue_Varanus_JD_infer_w, SNP, SnpPosition, CallRate, OneRatioRef, OneRatioSnp, FreqHomRef, FreqHomSnp, FreqHets, PICRef, PICSnp, AvgPIC, AvgCountRef, AvgCountSnp, RepAvg, clone, uid, rdepth, maf, monomorphs, OneRatio, PIC
@other$latlon[g]: coordinates for all individuals are attached
<- gl.keep.pop(vac.gl_178, c("AA", "MA", "MM")) vac.gl_178
Starting gl.keep.pop
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Checking for presence of nominated populations
Retaining only populations AA, MA, MM
Warning: Resultant dataset may contain monomorphic loci
Locus metrics not recalculated
Completed: gl.keep.pop
table(vac.gl_178$pop)
AA MA MM
5 7 22
# Concantenate the SNPs for scaf_178 for phylogenetic inference-----------------
# The .FASTA file was used as input sequences for Genieous Prime
<-
fasta gl2fasta(
vac.gl_178,method = 3,
outfile = "./data/vac.gl_178_concat.fasta",
outpath = getwd(),
verbose = 3
)
Starting gl2fasta
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Warning: Dataset contains monomorphic loci which will be included in the output fasta file
Assigning ambiguity codes to heterozygote SNPs, concatenating SNPs
Removing loci for which SNP position is outside the length of the trimmed sequences
Generating haplotypes ... This may take some time
Completed: gl2fasta
# Supplementary Figure 4--------------------------------------------------------
<- gl.filter.callrate(vac.gl_178, threshold = 1) vac.pops_temp
Starting gl.filter.callrate
Processing genlight object with SNP data
Warning: data include loci that are scored NA across all individuals.
Consider filtering using gl <- gl.filter.allna(gl)
Warning: Data may include monomorphic loci in call rate
calculations for filtering
Recalculating Call Rate
Removing loci based on Call Rate, threshold = 1
Completed: gl.filter.callrate
gl.smearplot(
vac.pops_temp,ind.labels = T,
plot.colors = c('royalblue',
'yellow', 'red'),
label.size = 6,
group.pop = T
)
Processing genlight object with SNP data
Starting gl.smearplot
Completed: gl.smearplot
vac.pops_temp
********************
*** DARTR OBJECT ***
********************
** 34 genotypes, 874 SNPs , size: 89.8 Mb
missing data: 0 (=0 %) scored as NA
** Genetic data
@gen: list of 34 SNPbin
@ploidy: ploidy of each individual (range: 2-2)
** Additional data
@ind.names: 34 individual labels
@loc.names: 874 locus labels
@loc.all: 874 allele labels
@position: integer storing positions of the SNPs [within 69 base sequence]
@pop: population of each individual (group size range: 5-22)
@other: a list containing: loc.metrics, ind.metrics, latlon, loc.metrics.flags, history, verbose
@other$ind.metrics: id, species, ind, pop, sex, lat, lon, e, n, karyo, PCR1, DAPI, FISH
@other$loc.metrics: AlleleID, CloneID, AlleleSequence, TrimmedSequence, Chrom_Komodo_dragon_NCBIv1, ChromPos_Komodo_dragon_NCBIv1, AlnCnt_Komodo_dragon_NCBIv1, AlnEvalue_Komodo_dragon_NCBIv1, Chrom_Varanus_JD_infer_z, ChromPos_Varanus_JD_infer_z, AlnCnt_Varanus_JD_infer_z, AlnEvalue_Varanus_JD_infer_z, Chrom_Varanus_JD_vac_aligned_chr6_7, ChromPos_Varanus_JD_vac_aligned_chr6_7, AlnCnt_Varanus_JD_vac_aligned_chr6_7, AlnEvalue_Varanus_JD_vac_aligned_chr6_7, Chrom_Varanus_JD_infer_w, ChromPos_Varanus_JD_infer_w, AlnCnt_Varanus_JD_infer_w, AlnEvalue_Varanus_JD_infer_w, SNP, SnpPosition, CallRate, OneRatioRef, OneRatioSnp, FreqHomRef, FreqHomSnp, FreqHets, PICRef, PICSnp, AvgPIC, AvgCountRef, AvgCountSnp, RepAvg, clone, uid, rdepth, maf, monomorphs, OneRatio, PIC
@other$latlon[g]: coordinates for all individuals are attached
Further Study
Link to github repository of GENESPACE: https://github.com/jtlovell/GENESPACE/tree/master
Dobry, Zhu, et al. (2023)
Dobry, Wapstra, et al. (2023)