23900000 SAM alignment record pairs processed.Warning: Mate records missing for 270457 records; first such record: <SAM_Alignment object: Paired-end read 'A00224:87:H5WHYDSXX:1:2474:9118:3098' aligned to utg9_pilon_pilon_pilon:[3553256,3553424)/+>.24000000 SAM alignment record pairs processed....
R语言代码有兴趣可以了解一下:
###########genotype#############################genotype<-function(){count<-0 if (population_structure=="RIL"){x<-runif(1) if (x<=0.5){count<-1}else{count<-0} }else{for(i in 1:2){x<-runif(1) if (x<=0.5){number<-0.5}else{n...