library(Matrix)

count <- 0
posA <- c()
refA <- c()
altA <- c()
GalleleI <- Matrix(0,ncol=0,nrow=10000,sparse=TRUE)
l0 <- 0
t1 <- 0

#########START

setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data30")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)

setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data31")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data32")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data33")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data34")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data35")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data36")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data37")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data38")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data39")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data310")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data311")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data312")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data313")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data314")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data315")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data316")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data317")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data318")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data319")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data320")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data321")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data322")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data323")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data324")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data325")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data326")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data327")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data328")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data329")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data330")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data331")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



#setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data332")
#load("res1.Rda")
#l1 <- length(alt)
#l0 <- l0+l1
#altA <- c(altA,alt)
#refA <- c(refA,ref)
#post <- pos1+t1
#posA <- c(posA,post)
#t1 <- posA[l0]
#GalleleI <- cbind2(GalleleI,hap)
##rm(alt)
#rm(ref)
#rm(pos1)
#rm(hap)
#rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data333")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data334")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data335")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



#setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data336")
#load("res1.Rda")
#l1 <- length(alt)
#l0 <- l0+l1
#altA <- c(altA,alt)
#refA <- c(refA,ref)
#post <- pos1+t1
#posA <- c(posA,post)
#t1 <- posA[l0]
#GalleleI <- cbind2(GalleleI,hap)
##rm(alt)
#rm(ref)
#rm(pos1)
#rm(hap)
#rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data337")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data338")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data339")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data340")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data341")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data342")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



#setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data343")
#load("res1.Rda")
#l1 <- length(alt)
#l0 <- l0+l1
#altA <- c(altA,alt)
#refA <- c(refA,ref)
#post <- pos1+t1
#posA <- c(posA,post)
#t1 <- posA[l0]
#GalleleI <- cbind2(GalleleI,hap)
##rm(alt)
#rm(ref)
#rm(pos1)
#rm(hap)
#rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data344")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data345")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data346")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data347")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data348")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data349")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data350")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data351")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data352")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data353")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data354")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data355")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data356")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data357")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data358")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data359")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data360")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data361")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data362")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data363")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data364")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data365")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data366")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data367")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data368")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data369")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data370")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data371")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data372")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data373")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data374")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data375")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data376")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data377")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data378")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data379")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data380")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data381")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data382")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data383")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data384")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data385")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data386")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data387")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data388")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data389")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data390")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data391")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data392")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data393")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data394")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data395")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data396")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data397")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data398")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data399")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3100")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3101")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3102")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3103")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3104")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3105")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3106")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3107")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3108")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3109")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3110")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3111")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3112")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3113")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3114")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3115")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3116")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3117")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3118")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3119")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3120")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3121")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3122")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3123")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3124")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



#setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3125")
#load("res1.Rda")
#l1 <- length(alt)
#l0 <- l0+l1
#altA <- c(altA,alt)
#refA <- c(refA,ref)
#post <- pos1+t1
#posA <- c(posA,post)
#t1 <- posA[l0]
#GalleleI <- cbind2(GalleleI,hap)
##rm(alt)
#rm(ref)
#rm(pos1)
#rm(hap)
#rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3126")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3127")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3128")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3129")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3130")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3131")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3132")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3133")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3134")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3135")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3136")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3137")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3138")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3139")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3140")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3141")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3142")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3143")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3144")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3145")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3146")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3147")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3148")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3149")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)



setwd("/seppdata/sepp/linkage/software/simulation/sfscode/data3150")
load("res1.Rda")
l1 <- length(alt)
l0 <- l0+l1
altA <- c(altA,alt)
refA <- c(refA,ref)
post <- pos1+t1
posA <- c(posA,post)
t1 <- posA[l0]
GalleleI <- cbind2(GalleleI,hap)
rm(alt)
rm(ref)
rm(pos1)
rm(hap)
rm(dataG)
count <- count+1
print(t1)
print(l0)
print(count)


setwd("/seppdata/sepp/linkage/software/simulation/sfscode/")

save(l0,refA,altA,posA,t1,GalleleI,file="allData.Rda")


freq <- colSums(GalleleI)/10000

> summary(freq)
    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
0.000000 0.000100 0.000100 0.007482 0.000300 0.500000 

> length(which(freq>0.0001))
[1] 450288
> length(which(freq>0.0002))
[1] 290749
> length(which(freq>0.0003))
[1] 223352
> length(which(freq>0.0005))
[1] 163389
> length(which(freq>0.001))
[1] 113957
> length(which(freq>0.002))
[1] 84912
> length(which(freq>0.005))
[1] 61426
> length(which(freq>0.01))
[1] 49613
> length(which(freq>0.02))
[1] 41098
> length(which(freq>0.05))
[1] 32363