摘要：#!/bin/bash echo "par1 is gene,par2 is fst_file" #awk '{(if ($3==gene) print $0}' $2>gene_info grep -i $1 gene_info|awk '{print $2,$3,$4}'>${1}_info w 阅读全文

摘要：#!/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=$(( 阅读全文