# --window-pi-step,可选参数 L_donovani_all_samples_10kb.windowed.pi(不添加步长结果文件) 加入步长100b(--window-pi-step 100) 7,计算(Fst) vcftools --vcf L_donovani_all_samples.vcf --weir-fst-pop population1 --weir-fst-pop population2 --fst-window-size 10000 --out pop1_vs_pop2_F...
–window-pi, --window-pi-step计算窗口中的核酸多样性值 FST计算 –weir-fst-pop FILE:file must contain a list of individuals (one individual per line) from the VCF file that correspond to one population,可多次,生成weir.fst文件 –fst-window-size, --fst-window-step FSTc滑动窗口计算,重测序一...
--window-pi ：指定窗口大小 --window-pi-step ： 指定步长大小 out.windowed.pi 4. vcftools去除...
--window-pi 10000 --window-pi-step ##pi,核苷酸多样性,越大说明核苷酸多样性越高,越低说明两个座位DNA序列差异越小 #计算pi,滑动窗口为10kb,可选窗口滑动step参数 ## Fst --weir-fst-pop 20data/10number.txt --weir-fst-pop 20data/20number.txt --fst-window-size 500000 ## Fst, 分化系数,从0...
###计算pivcftools--vcf merged.vcf--window-pi500000--outpi###计算TajimaDvcftools--vcf marged.vcf--TajimaD500000--outTajimaD###计算HWvcftools--vcf merged.vcf--hardy--outHW###fst需要两个以上popvcftools--vcf SNP.vcf--weir-fst-pop1.txt--weir-fst-pop2.txt--out1_VS_2--fst-window...
1. 提取基因型信息：使用--extract-FORMAT-info选项，vcftools可以提取特定的基因型信息。2. 提取指定区域的变异信息：通过参数--vcf、--chr、--form-bp和--to-bp，用户可以指定一个区域，vcftools将提取该区域内的变异信息，输出到指定的文件中。3. 对vcf文件进行划窗处理：使用--window-pi和--...
## sliding window Fst值(滑动窗口,窗口10kb,步长5kb) vcftools --vcf input_data.vcf --window-pi 10000 --window-pi-step 5000 --out pop1_vs_pop2 & 1. 2. 3. 4. 5. 转换成Plink格式: vcftools --vcf input_data.vcf --plink --chr 1 --out output_in_plink ...
–window-pi –window-pi-step FST计算 –weir-fst-pop<filename></filename> –fst-window-size –fst-window-step 其它计算 –het –hardy –site-quality 主要用于提取VCF文件中每个位点的QUAL信``` --missing-indv --missing-site 计算每个位点的缺失率 ...
–window-pi –window-pi-step Measures the nucleotide diversity in windows, with the number provided as the window size. The output file has the suffix “.windowed.pi”. The latter is an optional argument used to specify the step size in between windows.滑动窗口计算sliding windows,参考A refere...
vcftools --vcf chr_4_imputed.vcf --chr chr4 --from-bp 25300000 --to-bp 25400000 --keep maize.csv --window-pi 1000 --out maize_1000.pi 3. R软件可视化该区域 zemaz = read.table(file='~/maize_hapmap/maize_1000.pi',header = T) ...