摘要： #!/bin/bash #set -Eeuo pipefail if [ $# -lt 2 ] then echo "par1 is vcf,par2 is pvalue,sample_size=297" exit fi #vcftools --vcf $1 --freq --out ${1}maf 阅读全文

摘要： #!/bin/bash #LD PLOT #author lee echo"par1 is vcf,par2 is top_chr" Genome=/public/home/caisl/lee/genome/rice/msu.gff3 while read chr pos do pos_up=$(( 阅读全文